Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 781,670 A→T 8.2% intergenic (+18/‑415) lysQ → / → nadA tRNA‑Lys/quinolinate synthase, subunit A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913781,6700AT 8.2% 238.2 / 13.6 95intergenic (+18/‑415)lysQ/nadAtRNA‑Lys/quinolinate synthase, subunit A
Reads supporting (aligned to +/- strand):  ref base A (50/36);  new base T (6/2);  total (57/38)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAA  >  NC_000913/781523‑781816
                                                                                                                                                   |                                                                                                                                                  
gTCATCCTGCCCGACCCTCCATCCTGACTGCTTTAGGCAGGATAAGCCCGCCAGGGGTCGTTAGCTGTGGTGGTAGTGCTGTGGTCTTTTATTCTGGGGGTCGCAGGTTGGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTaaa                                                                                                                                                  <  2:1224849/150‑1 (MQ=12)
       tGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtt                                                                                                                                           >  2:866229/1‑150 (MQ=255)
       tGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtt                                                                                                                                           >  2:1441121/1‑150 (MQ=255)
       tGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtt                                                                                                                                           >  2:1086440/1‑150 (MQ=255)
         ctcGACCCACCATCCTTAATGATTAAGTCATCATAATCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttt                                                                                                                                         <  2:1361570/148‑1 (MQ=255)
         cACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGGAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttt                                                                                                                                         <  2:763534/150‑1 (MQ=255)
         cACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCATGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttt                                                                                                                                         <  2:1361728/150‑1 (MQ=255)
          ccgcccccccATCCTGAATGATTAAGGCAGCATACGCCCGCATGGGGTCGTTAGCTCCGTTGGTAGTGCAGTTGACTTTTAATCAATTGGTCGCTGGTGCGAATCCTGCACGACCCACCACTGTAAAAAAGCGCCTTTAGGGCGCttttt                                                                                                                                        <  1:402252/146‑1 (MQ=17)
           cGACCCACCATCCTGTATGATTAAGTCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCtttttt                                                                                                                                       <  1:2136437/150‑1 (MQ=255)
            gACCCACCATCCTGAATTATTAATGCAGCATAATCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTg                                                                                                                                      <  1:1908776/150‑1 (MQ=255)
            gACCCACCATCCTGAATTATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCTTTAGGGCGCtttttta                                                                                                                                      <  1:1666769/150‑2 (MQ=255)
                caccaTCCTGACTGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                  <  2:379211/150‑1 (MQ=255)
                caccaTCCTGAATGATTAAGGCAGCATAATCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                  <  2:1356015/150‑1 (MQ=255)
                caccaTCCTGAATGACTACGGCAGCATAATCCCGCTAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTAt                                                                                                                                  <  2:1356001/150‑1 (MQ=255)
                        tGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCATTTGGTAGAGCAGTTGACTTTTAAAAAATTGGTCGCAGGTTAGAATCCTTCACGACCCACCAAAGTAAAAAATCGCCCTAAATGCCCATTTTTGCTATCTGCGATa                                                                                                                          >  1:1374221/1‑150 (MQ=17)
                         gAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGAAAc                                                                                                                         >  2:2269421/1‑150 (MQ=255)
                             gATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTTCGAAGCTCaa                                                                                                                     >  2:1442877/1‑150 (MQ=255)
                             gATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATATGCGATACTCaa                                                                                                                     >  2:2229477/1‑150 (MQ=255)
                              aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCGATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCcaa                                                                                                                    >  2:1973638/1‑150 (MQ=255)
                              aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                    <  1:1599681/150‑1 (MQ=255)
                              aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                    >  2:2454536/1‑150 (MQ=255)
                              aTTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCaaa                                                                                                                    >  2:707560/1‑150 (MQ=255)
                               ttAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTCAAGGCGCTTTTTTGCTATCTGCGAGACTCAAAg                                                                                                                   >  1:1998697/1‑150 (MQ=255)
                               ttAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAg                                                                                                                   >  1:1998686/1‑150 (MQ=255)
                                  aGGCAGCCTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGAtt                                                                                                                <  1:49970/150‑1 (MQ=255)
                                  aGGCAGCATAATCCCGCCAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGAtt                                                                                                                <  1:191127/150‑1 (MQ=255)
                                  aGGCACCATAATCCCGCTTGGGGTCGTTCTCTCATTTGGTAGAGCATTCGACTTTAAATCCTTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGAtt                                                                                                                <  1:488027/150‑1 (MQ=14)
                                      agcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGaa                                                                                                            >  2:75919/1‑150 (MQ=255)
                                      agcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGaa                                                                                                            <  2:591721/150‑1 (MQ=255)
                                      agcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTAAAAGATTCGaa                                                                                                            >  2:75787/1‑150 (MQ=255)
                                      agcaTAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTGGAATCAAGGGGTCGCCGGGTCGACTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCGGGTGTGCTAGCGGCGATACGCCAAGaggcgac                                                                                                            >  2:292828/1‑144 (MQ=14)
                                             taccggtaGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                     <  2:1951831/143‑1 (MQ=18)
                                             tCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgcag                                                                                                     <  1:1693428/150‑1 (MQ=39)
                                               ccGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTgaagca                                                                                                   >  2:1668521/1‑150 (MQ=38)
                                               ccGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCACCGATACTAAAATATGCGAAACTgcacca                                                                                                   >  2:1107295/1‑150 (MQ=11)
                                                cGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAAATTGTAGCAGGTTCGAATCCTGCACGACCCACCAAAGTAAAAAATCGCCCGAAAGGCGCTTTTTTTCTATCTGCGATACTCAAAGATTCGAACcagcagcag                                                                                                  >  1:643601/1‑150 (MQ=11)
                                                  actcgggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCAAAAGGCGCTTTTTTGCTAAAAGCGATACTCAAAGATTCGAACCTGCAGCAggg                                                                                                >  2:1755494/5‑149 (MQ=11)
                                                      gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                            <  1:1755494/150‑1 (MQ=32)
                                                      gggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttga                                                                                            >  2:107379/1‑150 (MQ=32)
                                                       ggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                                                           <  1:1866242/150‑1 (MQ=32)
                                                       ggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                                                           <  1:787930/150‑1 (MQ=32)
                                                       ggTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                                                           <  1:12895/150‑1 (MQ=32)
                                                       ggTCGTTAGCTCAGTTGGTACAGCAGTTGACTTTTAATCAATTGGTCGCAGGTACGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgag                                                                                           <  1:1906065/150‑1 (MQ=17)
                                                         tCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATAAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATATGCGATACTCAAAGATTCGAAACTGCAGCAGGTttgagtt                                                                                         >  1:765114/1‑150 (MQ=12)
                                                           gTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                       <  1:264588/150‑1 (MQ=17)
                                                           gTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                       <  1:264585/150‑1 (MQ=17)
                                                           gTTAGCTCCGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTttgagttga                                                                                       <  1:1006442/150‑1 (MQ=17)
                                                                cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                  >  1:2432370/1‑150 (MQ=32)
                                                                cTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCa                                                                                  >  1:1743595/1‑150 (MQ=32)
                                                                   aGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCg                                                                               <  1:360121/150‑1 (MQ=32)
                                                                    gTTGGTAGCGCAGTTGACTTTTCATCACTTGGGCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAACAAGCGCCCTACAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGa                                                                              <  2:1191550/150‑1 (MQ=11)
                                                                          agagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGATGAGCGCAGCGAaacacc                                                                        >  1:1263289/1‑148 (MQ=17)
                                                                           gagCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTCTCAGCGATACTCAAAGATTCGAACCTGCCGCAGGTGTGAGTTGAGCGCAGCGAAACAACg                                                                       >  2:2425811/1‑150 (MQ=255)
                                                                            agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACgg                                                                      >  2:410122/1‑150 (MQ=32)
                                                                            agCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATACGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACgg                                                                      >  2:281553/1‑150 (MQ=17)
                                                                                gTTGTCTTTTAATCATTTGGTCGCAGGTTCGAATCCTGCACGACCCACAAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAATTTTCGAACCTGCAGCAGGTTTTAGTTGAGCGCAGCTAAACAACGGAGcc                                                                  <  2:1553320/150‑1 (MQ=11)
                                                                                gTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACTAACCACCAATGTAAAAAAGCGACCTAAAGTCGCTTTTTTGCTATATGCTATACTCAAAGATTCGAACCTGCAGCAGTTTTGAGTTGATCTCAGCGATACAACGGAgac                                                                  >  1:972322/1‑148 (MQ=11)
                                                                                   gACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCt                                                               >  1:1328716/1‑150 (MQ=32)
                                                                                       tttAATCATTTGGTCGCAGGTTCTAATCCTGCACGACCCACCAATGTAAAAATGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAATTATTCGAACCTGCAGCAGGTTTTTGTTGAGCGCAGCGAAACAACGGAGCCGCTcgcg                                                           <  2:885541/150‑1 (MQ=11)
                                                                                           aTCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACGACGGAGCCGCTCgcggcgc                                                       >  2:1442594/1‑149 (MQ=17)
                                                                                              acttGGTCGCAGGTTCGACTCCTGCACGACCCCCCACTGTAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                    <  2:1743983/148‑1 (MQ=255)
                                                                                              aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCACTGTAAAAAAGCGCCCTACAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                    <  2:1619191/150‑1 (MQ=17)
                                                                                              aaTTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACgg                                                    <  2:1619220/150‑1 (MQ=32)
                                                                                                      gCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAACGGCGCTTTTTTGCTATCTGCGATACTCAAATATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAgg                                            <  2:107018/150‑1 (MQ=12)
                                                                                                      gCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGTTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAgg                                            <  2:107040/150‑1 (MQ=17)
                                                                                                           ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCGAAGATTCGACCCTGCCGGAGCGGTGAGTTGAGCGCCGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGATGGgcgag                                       >  1:854546/1‑150 (MQ=11)
                                                                                                           ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGACGGgcgag                                       >  2:328171/1‑150 (MQ=17)
                                                                                                           ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                       >  2:505839/1‑150 (MQ=32)
                                                                                                           ttCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgag                                       >  1:864806/1‑150 (MQ=32)
                                                                                                              gAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                    >  2:1521121/1‑150 (MQ=32)
                                                                                                              gAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAAAAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                    >  1:1498927/1‑150 (MQ=17)
                                                                                                              gAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGAAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGgcgagcga                                    >  1:1049051/1‑150 (MQ=17)
                                                                                                               taTCCTGCACGACCCACCACTGTAAACATGCGCCGTAAAGGCGCTTTTTTGCTATCTGCGTTACTCAACGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGCAGGgcgagcgac                                   <  2:2534887/149‑2 (MQ=11)
                                                                                                               aaTCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCGTTTTTGCTATCGGCGATACTCAAAGATTAGAACCTGCAGCAGGTTTGAGTTGTGCGCCGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGCAGGGCGagcgaa                                   >  1:845174/1‑150 (MQ=11)
                                                                                                                 tCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGAAACAGCAATATGAGAAGCTGCAGCAGGTTAGAGGGTAGCGCATCGAACCGACGGAGCCGCTCGCGCCGACGGCCCGGGGGGCGagagacgc                                 >  2:2076767/1‑150 (MQ=2)
                                                                                                                  ccTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGagcgaagcg                                >  2:1884215/1‑150 (MQ=32)
                                                                                                                  ccTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCGGCGATACTCAAAGATGCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACACCGGAGCCGCTCGCGGCGACGGCCCGGGGGGCGagcgaagcg                                >  2:1884183/1‑150 (MQ=11)
                                                                                                                              accaATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCa                    <  2:643601/150‑1 (MQ=32)
                                                                                                                                   tGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGAcc               >  1:638243/1‑150 (MQ=32)
                                                                                                                                     tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGCGTCATCCTGCACGACcca             >  1:1790327/1‑150 (MQ=17)
                                                                                                                                     tAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca             >  1:320791/1‑150 (MQ=32)
                                                                                                                                     caaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACcca             >  2:1427141/2‑150 (MQ=12)
                                                                                                                                      aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac            >  2:1690145/1‑150 (MQ=12)
                                                                                                                                      aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac            >  2:2139923/1‑150 (MQ=12)
                                                                                                                                      aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac            >  2:1689999/1‑150 (MQ=12)
                                                                                                                                      aaaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTAGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac            >  2:303745/1‑150 (MQ=255)
                                                                                                                                      aaaaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccac            >  1:1831401/1‑150 (MQ=32)
                                                                                                                                       aaaaaGCGCCTTTAGGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacc           >  1:2307365/1‑150 (MQ=12)
                                                                                                                                        aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca          >  1:847538/1‑150 (MQ=32)
                                                                                                                                        aaaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACccacca          >  1:1468114/1‑150 (MQ=32)
                                                                                                                                         aaaGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTGTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCaa         >  2:1013508/1‑150 (MQ=17)
                                                                                                                                              gcCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCTGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaa    <  1:1441121/150‑1 (MQ=17)
                                                                                                                                              gcCCTAAAGGCGATTTTTTGATATTTGCTATACTAAAATTTTCTAACCTGCTGCTGTTTTGTTTTGCGCGCCGCTAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaa    <  1:1086440/150‑1 (MQ=11)
                                                                                                                                                ccTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGCCCCACCAATGTaaaa  >  1:1653709/1‑150 (MQ=17)
                                                                                                                                                ccTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa  >  2:264746/1‑150 (MQ=32)
                                                                                                                                                ccTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTaaaa  >  1:1423433/1‑150 (MQ=32)
                                                                                                                                                   |                                                                                                                                                  
GTCATCCTGCACGACCCACCATCCTGAATGATTAAGGCAGCATAATCCCGCAAGGGGTCGTTAGCTCAGTTGGTAGAGCAGTTGACTTTTAATCAATTGGTCGCAGGTTCGAATCCTGCACGACCCACCAATGTAAAAAAGCGCCCTAAAGGCGCTTTTTTGCTATCTGCGATACTCAAAGATTCGAACCTGCAGCAGGTTTGAGTTGAGCGCAGCGAAACAACGGAGCCGCTCGCGGCGACGGCCCGAAGGGCGAGCGAAGCGAGTCATCCTGCACGACCCACCAATGTAAAA  >  NC_000913/781523‑781816

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