Predicted mutation
evidence seq id position mutation annotation gene description
RA 3 30,670 A→G intergenic (+130/+189) unknown → / ← unknown S‑ribosylhomocysteine lyase (EC 4.4.1.21) @ Autoinducer‑2 production protein LuxS/LSU rRNA pseudouridine(746) synthase (EC 5.4.99.29) @ tRNA pseudouridine(32) synthase (EC 5.4.99.28)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*330,6700AG98.7% 219.9 / ‑3.0 78intergenic (+130/+189)unknown/unknownS‑ribosylhomocysteine lyase (EC 4.4.1.21) @ Autoinducer‑2 production protein LuxS/LSU rRNA pseudouridine(746) synthase (EC 5.4.99.29) @ tRNA pseudouridine(32) synthase (EC 5.4.99.28)
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (32/45);  minor base . (1/0);  total (33/45)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.18e-01

TACACCGAGCATTCTTTGGAAGAAGCACATCAAATCGCACAAAATGTGTTGGATCGCGGTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGAAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAACCCGACCAGTGTGGTGAATTCCATGGTTTCGCCGGTGATCGGGTGCCGGA  >  3/30424‑30914
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
tACACCGAGCATTCTTTGGAAGAAGCACATCAAATCGCACAAAATGTGTTGGATCGCGGTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAgg                                                                                                                                                                                                                                                  <  1:307681/251‑1 (MQ=255)
 acacCGAGCATTCTTTGGAAGAAGCACATCAAATCGCACAAAATGTGTTGGATCGCGGTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGc                                                                                                                                                                                                                                                 <  1:639714/251‑1 (MQ=255)
              tttGGAAGAAGCACATCAAATCGCACAAAATGTGTTGGATCGCGGTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGa                                                                                                                                                                                                                                    <  1:276421/251‑1 (MQ=255)
                          acaTCAAATCGCACAAAATGTGTTGGATCGCGGTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATaa                                                                                                                                                                                                                        <  1:29915/251‑1 (MQ=255)
                              cAAATCGCACAAAATGTGTTGGATCGCGGTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGAc                                                                                                                                                                                                                    <  1:387230/251‑1 (MQ=255)
                              cAAATCGCACAAAATGTGTTGGATCGCGGTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGAc                                                                                                                                                                                                                    <  1:365418/251‑1 (MQ=255)
                                   cGCACAAAATGTGTTGGATCGCGGTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACaataa                                                                                                                                                                                                               <  1:118625/251‑1 (MQ=255)
                                   cGCACAAAATGTGTTGGATCGCGGTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACaataa                                                                                                                                                                                                               <  1:44672/251‑1 (MQ=255)
                                       caAAATGTGTTGGATCGCGGTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTg                                                                                                                                                                                                            >  1:190754/1‑250 (MQ=255)
                                             gtgtTGGATCGCGGTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTgaaaaaga                                                                                                                                                                                                     >  1:457669/1‑251 (MQ=255)
                                                  ggATCGCGGTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGAc                                                                                                                                                                                                                    <  1:358736/231‑1 (MQ=255)
                                                        cgGTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACaataa                                                                                                                                                                                                               >  1:703765/1‑230 (MQ=255)
                                                          gTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGca                                                                                                                                                                                        >  1:555348/1‑251 (MQ=255)
                                                          gTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGc                                                                                                                                                                                         <  1:854217/250‑1 (MQ=255)
                                                          gTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTAAAAAAGAAAAATATTTGGca                                                                                                                                                                                        >  1:490262/1‑251 (MQ=255)
                                                              cGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGcaatac                                                                                                                                                                                    <  1:47730/251‑1 (MQ=255)
                                                                      ataaaAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGa                                                                                                                                                                            <  1:820506/251‑1 (MQ=255)
                                                                       taaaAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGAc                                                                                                                                                                           <  1:159047/251‑1 (MQ=255)
                                                                                  aTTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGc                                                                                                                                                                  <  1:157549/249‑1 (MQ=255)
                                                                                       aCCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATaa                                                                                                                                                                               >  1:60140/1‑231 (MQ=255)
                                                                                       aCCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGa                                                                                                                                                                            >  1:328023/1‑234 (MQ=255)
                                                                                                  aaGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTg                                                                                                                                                <  1:538209/251‑1 (MQ=255)
                                                                                                          taAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGAc                                                                                                                                                                                                                    <  1:112354/175‑1 (MQ=255)
                                                                                                              aCAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAgg                                                                                                                                                                                                                      >  1:514759/1‑169 (MQ=255)
                                                                                                              aCAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGc                                                                                                                                     >  1:204968/1‑250 (MQ=255)
                                                                                                               cAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAg                                                                                                                                                                                                                                     <  1:533982/153‑1 (MQ=255)
                                                                                                               cAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGc                                                                                                                                     >  1:750635/1‑249 (MQ=255)
                                                                                                                 aTAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCg                                                                                                                                 >  1:151834/1‑251 (MQ=255)
                                                                                                                      gagaCGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGgaaa                                                                                                                             <  1:554495/250‑1 (MQ=255)
                                                                                                                         aCGTTTATGACCGCCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAAGAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGgaaaaga                                                                                                                          <  1:741267/250‑1 (MQ=255)
                                                                                                                          cGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAAtctc                                                                                                                                                                                                                                          >  1:803805/1‑137 (MQ=255)
                                                                                                                           gTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGcaatac                                                                                                                                                                                    <  1:309779/190‑1 (MQ=255)
                                                                                                                           gTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGgaaaagaaaa                                                                                                                       >  1:82567/1‑251 (MQ=255)
                                                                                                                                    caccaTCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGaaaagaaaagag                                                                                                                    >  1:680910/1‑243 (MQ=255)
                                                                                                                                      ccaTCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGa                                                                                                            >  1:700694/1‑251 (MQ=255)
                                                                                                                                         tCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGc                                                                                                                                     >  1:288054/1‑223 (MQ=255)
                                                                                                                                         tCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGgaaa                                                                                                                             <  1:241527/231‑1 (MQ=255)
                                                                                                                                         tCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAAtt                                                                                                         >  1:467745/1‑251 (MQ=255)
                                                                                                                                            gTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGc                                                                                                                                                                                                                                                 >  1:814281/1‑112 (MQ=255)
                                                                                                                                            gTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGc                                                                                                                                                                                         <  1:254634/168‑1 (MQ=255)
                                                                                                                                            gTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCaa                                                                                                                                                                     <  1:353565/188‑1 (MQ=255)
                                                                                                                                            gTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTAt                                                                                                                                        >  1:299906/1‑217 (MQ=255)
                                                                                                                                            gTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGgaaaagaaa                                                                                                                        <  1:460586/233‑1 (MQ=255)
                                                                                                                                              aCATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAgg                                                                                                    <  1:305717/251‑1 (MQ=255)
                                                                                                                                                   aCCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAAc                                                                                                                                                                        <  1:606427/178‑1 (MQ=255)
                                                                                                                                                    ccTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGttt                                                                                               >  1:291262/1‑250 (MQ=255)
                                                                                                                                                        aaaGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACAt                                                                                          <  1:853607/251‑1 (MQ=255)
                                                                                                                                                                    aTTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGc                                                                                                                                                                  >  1:153940/1‑167 (MQ=255)
                                                                                                                                                                              ggCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATc                                                                     >  1:404397/1‑250 (MQ=255)
                                                                                                                                                                                     ggATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATG‑GAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGgaaaagaaa                                                                                                                        >  1:635841/1‑191 (MQ=255)
                                                                                                                                                                                            gATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTaa                                                      >  1:56609/1‑251 (MQ=255)
                                                                                                                                                                                               tGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTACAAc                                                   <  1:876937/251‑1 (MQ=255)
                                                                                                                                                                                               tGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAAc                                                   <  1:852660/251‑1 (MQ=255)
                                                                                                                                                                                               tGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAAc                                                   <  1:438419/251‑1 (MQ=255)
                                                                                                                                                                                                 ccGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAAcc                                                  <  1:27250/250‑1 (MQ=255)
                                                                                                                                                                                                 ccGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAAcc                                                  <  1:453707/250‑1 (MQ=255)
                                                                                                                                                                                                    gTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGc                                                                                                                                                                  >  1:607140/1‑135 (MQ=255)
                                                                                                                                                                                                    gTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGaaacgggttc                                                                                                                                                    <  1:186059/149‑9 (MQ=255)
                                                                                                                                                                                                    gTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAACCCGa                                               >  1:300597/1‑250 (MQ=255)
                                                                                                                                                                                                              ccAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAACCCGACCAGTgtggtg                                    <  1:656232/251‑1 (MQ=255)
                                                                                                                                                                                                              ccAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAACCCGACCAGTgtggtg                                    <  1:342408/251‑1 (MQ=255)
                                                                                                                                                                                                               cAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAACCCGACCAGTGTGGTGa                                   <  1:870757/251‑1 (MQ=255)
                                                                                                                                                                                                               cAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAACCCGACCAGTGTGGTGa                                   <  1:255951/251‑1 (MQ=255)
                                                                                                                                                                                                               cAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAACCCGACCAGTGTGGTGa                                   <  1:767550/251‑1 (MQ=255)
                                                                                                                                                                                                                  cTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGa                                                                                                                                                                            >  1:681653/1‑111 (MQ=255)
                                                                                                                                                                                                                         aGCCTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGa                                                                                                                                                                            >  1:528239/1‑104 (MQ=255)
                                                                                                                                                                                                                           ccTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAACCCGACCAGTGTGGTGAATTCCATgg                          <  1:423401/248‑1 (MQ=255)
                                                                                                                                                                                                                           ccTGTCCCCTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAACCCGACCAGTGTGGTGAATTCCATGGt                         <  1:503497/249‑1 (MQ=255)
                                                                                                                                                                                                                                   cTTCCGGCCCGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAACCCGACCAGTGTGGTGa                                   <  1:176201/231‑1 (MQ=255)
                                                                                                                                                                                                                                          cccGCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAACCCGACCAGTGTGGTGAATTCCATGGTTTCGCCGGTGATCggg         <  1:557252/250‑1 (MQ=255)
                                                                                                                                                                                                                                             gCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATc                                                                     >  1:777999/1‑187 (MQ=255)
                                                                                                                                                                                                                                             gCCATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAACCCGACCAGTGTGGTGAAt                                 <  1:819932/223‑1 (MQ=255)
                                                                                                                                                                                                                                              ccATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAAc                                                   <  1:180776/204‑1 (MQ=255)
                                                                                                                                                                                                                                              ccATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAACCCGACCAGTGTGGTGa                                   <  1:687458/220‑1 (MQ=255)
                                                                                                                                                                                                                                              ccATCATGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAACCCGACCAGTGTGGTGAATTCCATGGTTTCGCCGGTGATCGGGTGCCg    >  1:364365/1‑251 (MQ=255)
                                                                                                                                                                                                                                                atcatGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTc                                                                                              >  1:598908/1‑159 (MQ=255)
                                                                                                                                                                                                                                                atcatGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAACCCGACCAGTGTGGTGAATTCCATGGTTTCGCCGGTGATCGGGTGCCGGa  >  1:639523/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 tcatGGAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAacac                                                                               >  1:143943/1‑173 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    
TACACCGAGCATTCTTTGGAAGAAGCACATCAAATCGCACAAAATGTGTTGGATCGCGGTATCGGCGTAAATAAAAATGAAGATTTAACCTTAGACGAAAGCTTATTAAAACAATAAGGAGACGTTTATGACCACCATCGGTACATAACCTGAAAGCATTCGTCATTTTTTAGTGGCAGAAGGATTCGGATTGCCGGTTCTAATCGCCAACTTCATCAGCCTGTCCCCTTCCGGCCCGCCATCATGAAAGGCAAATCTCAAAAGATCGGAATGATAAGGACAATAATTGAAAAAGAAAAATATTTGGCAATACAATAAAGACAACCAAAGCAATGAACCCGTTTCATTGCTTTTTATCGCCTCGGAAAAGAAAAGTGCGGTCGGAATTTTAGGTGTTTCACATCGAAAAAACACAGGAAAAATCCACCGCACTTTTTAAAACCCGACCAGTGTGGTGAATTCCATGGTTTCGCCGGTGATCGGGTGCCGGA  >  3/30424‑30914

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑