Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,948,180 T→A 6.5% *151C (TGA→TGT nudB ← dihydroneopterin triphosphate pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,948,1800TA 6.5% 298.1 / 6.4 105*151C (TGA→TGTnudBdihydroneopterin triphosphate pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base T (42/56);  new base A (3/4);  total (45/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 106 total aligned reads displayed.

ACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACC  >  NC_000913/1948032‑1948328
                                                                                                                                                    |                                                                                                                                                    
aCCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTc                                                                                                                                                     <  1:1235082/150‑1 (MQ=255)
  cGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAg                                                                                                                                                   <  2:891525/150‑1 (MQ=255)
   gCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAgg                                                                                                                                                  >  2:552060/1‑150 (MQ=255)
    cccAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGc                                                                                                                                                 <  2:723506/150‑1 (MQ=255)
      cAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAg                                                                                                                                               >  1:1788595/1‑150 (MQ=255)
       aGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGc                                                                                                                                              >  2:1662609/1‑150 (MQ=255)
        tcTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGATTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAATATATCTCCAGATAGCCCTGCCTGAACAGGCAggg                                                                                                                                             <  2:1783530/149‑3 (MQ=255)
        gCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCg                                                                                                                                             >  1:1532783/1‑150 (MQ=255)
         cTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGt                                                                                                                                            >  1:1517383/1‑150 (MQ=255)
         cTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGAACAGGCAgggc                                                                                                                                            <  2:1509179/150‑4 (MQ=255)
         cTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGAACAGGCAgggc                                                                                                                                            <  2:1325091/150‑4 (MQ=255)
          ttAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGtt                                                                                                                                           <  2:1884129/150‑1 (MQ=255)
          ttAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGtt                                                                                                                                           <  2:1759659/150‑1 (MQ=255)
           tAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTa                                                                                                                                          <  1:1428499/150‑1 (MQ=255)
           tAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTa                                                                                                                                          <  2:1431086/150‑1 (MQ=255)
             gCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAAt                                                                                                                                        <  2:1530364/150‑1 (MQ=255)
             gCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGAACAGGCAgggctatc                                                                                                                                        <  1:2036695/150‑8 (MQ=255)
              ccGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAAtt                                                                                                                                       <  2:1228716/150‑1 (MQ=255)
              ccGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAAtt                                                                                                                                       <  2:986534/150‑1 (MQ=255)
                gcgGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTAc                                                                                                                                     <  2:1147725/150‑1 (MQ=255)
                  ggTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACaa                                                                                                                                   <  1:1816460/150‑1 (MQ=255)
                       ccAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGt                                                                                                                              >  2:554083/1‑150 (MQ=255)
                         aGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTc                                                                                                                            <  2:1603814/150‑1 (MQ=255)
                          gCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCa                                                                                                                           <  1:259675/150‑1 (MQ=255)
                            tCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAAt                                                                                                                         <  1:252788/150‑1 (MQ=255)
                             cACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATc                                                                                                                        >  1:994660/1‑150 (MQ=255)
                             cACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATc                                                                                                                        >  2:2023215/1‑150 (MQ=255)
                             cACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATc                                                                                                                        >  2:1176823/1‑150 (MQ=255)
                                     aTTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCg                                                                                                                >  2:1114973/1‑150 (MQ=255)
                                          aGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGc                                                                                                           >  1:318928/1‑150 (MQ=255)
                                          aGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGc                                                                                                           <  1:685805/150‑1 (MQ=255)
                                            tGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGCCCTGCCATAAAAATATATCCAGATAGCCCTGCCTGTTCCGGCAGCGTTAATTACAAACTGTAAATACCCCTGCCGTTTGCTc                                                                                                         >  2:1649397/1‑150 (MQ=255)
                                             gAAGATTTTACCGCGCTTAGCATCATGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCAGTTCAGGCAGCGTTAATTACAAACTGATCAATCGCCTGCCGGTTGCTcc                                                                                                        >  2:1273168/1‑150 (MQ=255)
                                              aaGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCa                                                                                                       <  1:1989924/150‑1 (MQ=255)
                                              aaGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGCTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCa                                                                                                       <  2:735449/150‑1 (MQ=255)
                                                  ttttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGCCCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTTATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGAc                                                                                                   >  2:1849170/1‑150 (MQ=255)
                                                  ttttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAggcc                                                                                                   >  1:773663/1‑148 (MQ=255)
                                                  ttttACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGAc                                                                                                   <  1:1435868/150‑1 (MQ=255)
                                                     tACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTa                                                                                                <  2:2048822/150‑1 (MQ=255)
                                                      aCCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAATTATCTCCTGATAGCCGTGCCTGTTCAGGCAGCGGTAATTACAATCTGTTGTATCGCCTGCCGGTTGCTCCAGGACTTAg                                                                                               <  1:259419/150‑1 (MQ=255)
                                                      aCCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAg                                                                                               <  1:75417/150‑1 (MQ=255)
                                                       ccGCGCTTAGCATCCTGCGCAGCTTTACGATGGCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGt                                                                                              <  2:982186/150‑1 (MQ=255)
                                                        cgcgCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTg                                                                                             <  1:1986814/150‑1 (MQ=255)
                                                         gcgcTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGa                                                                                            >  2:1142449/1‑150 (MQ=255)
                                                           gcTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAgc                                                                                          >  1:710400/1‑150 (MQ=255)
                                                                     cTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCgcag                                                                                >  1:332209/1‑150 (MQ=255)
                                                                      tGCGCAGCTTTACGCTGTCTGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAATCTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCgcagca                                                                              <  1:1342288/150‑1 (MQ=255)
                                                                      tGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCgcagc                                                                               >  1:1754962/1‑150 (MQ=255)
                                                                       gcgcAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCgcagca                                                                              <  1:780646/150‑1 (MQ=255)
                                                                        cgcAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGcagcag                                                                             >  1:803324/1‑150 (MQ=255)
                                                                             cTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCa                                                                        <  2:804583/150‑1 (MQ=255)
                                                                                tACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCa                                                                     >  2:1957951/1‑150 (MQ=255)
                                                                                  cGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAg                                                                   <  2:896772/150‑1 (MQ=255)
                                                                                     tGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGccc                                                                >  1:174506/1‑149 (MQ=255)
                                                                                      gTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCAc                                                               <  2:893029/150‑1 (MQ=255)
                                                                                      gTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCAc                                                               <  1:1294785/150‑1 (MQ=255)
                                                                                       tCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACt                                                              <  2:275327/150‑1 (MQ=255)
                                                                                          ggTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGTTAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGt                                                           <  2:966360/150‑1 (MQ=255)
                                                                                           gtgtTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTa                                                          <  1:826610/150‑1 (MQ=255)
                                                                                             gtTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCCCCGCAGCAGGCGCATCAAGCCACTTGTAAg                                                        <  1:855284/150‑1 (MQ=255)
                                                                                               tGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCCCCGCAGCAGGCGCATCAAGCCACTTGTAAGcc                                                      <  2:53259/150‑1 (MQ=255)
                                                                                                ggCCCATTTACAATGACCTGCCATAAAAATATCTCCAGATACCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGcct                                                     <  2:50174/150‑2 (MQ=255)
                                                                                                  cccATTTACTATGACCTGTCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCATGCCACTTGTAAGCCAgt                                                   >  2:199922/1‑149 (MQ=255)
                                                                                                  cccATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGa                                                   <  1:1176824/150‑1 (MQ=255)
                                                                                                  cccATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGa                                                   >  1:1181580/1‑150 (MQ=255)
                                                                                                  cccATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATAACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGa                                                   <  2:258024/150‑1 (MQ=255)
                                                                                                    cATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATg                                                 <  2:988537/150‑1 (MQ=255)
                                                                                                    cATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATg                                                 <  2:1535714/150‑1 (MQ=255)
                                                                                                       ttACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTc                                              <  2:512152/150‑1 (MQ=255)
                                                                                                        tACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCATTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCa                                             >  2:802510/1‑150 (MQ=255)
                                                                                                          cTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGt                                           <  1:1482987/150‑1 (MQ=255)
                                                                                                          cTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGt                                           >  2:294694/1‑150 (MQ=255)
                                                                                                                ccTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCACGCCCCTTGTAAGCCAGATGTTCAGTGAAAAc                                     >  1:1035887/1‑150 (MQ=255)
                                                                                                                  tGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGa                                   >  2:1669720/1‑150 (MQ=255)
                                                                                                                   gCCATAAAATTATCTCCAGATAGCCCTGCGTGTTCCGGCAGCGTTAATTAAAAACTGTTCAATCGCCTGCCGGTTGCTCCCGCACTTAGTGAGCGCAGCCGCAGCAGGCGCATCAAGCCCATTCTAAGCCAGATGTGCAGTCACAACGag                                  >  2:190343/1‑149 (MQ=255)
                                                                                                                   gCCATAAAAATATCTCCAGATAGCCCTGCATGTTAAGTCAGCGTTAATTACAAACTATTCAATCGCCTGCACGATGCTCCAGGACTTAGTGAGCGCCGGCGCAGCACCCGGATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGAt                                  >  2:1909093/1‑150 (MQ=255)
                                                                                                                    ccATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATc                                 <  1:914101/150‑1 (MQ=255)
                                                                                                                    ccATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATc                                 <  1:1358799/150‑1 (MQ=255)
                                                                                                                     cATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCt                                <  1:12500/150‑1 (MQ=255)
                                                                                                                     cATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCt                                <  1:1431186/150‑1 (MQ=255)
                                                                                                                       tAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCCGCGTTAATTACAAACTGGTCAATCGCCTGCCGGTTGCGCCAGGACTTAGTGAGCGCCGCCGCAGAAGGCGCATCACGCCGCTTGTAAGCCAGAGGTGCAGTGAAAACGATCTGc                              >  2:367354/1‑150 (MQ=255)
                                                                                                                        aaaaaTATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGcc                             <  1:1295850/150‑1 (MQ=255)
                                                                                                                        aaaaaTATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGcc                             <  2:273205/150‑1 (MQ=255)
                                                                                                                         aaaaTATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCg                            >  1:1428686/1‑150 (MQ=255)
                                                                                                                                ctcCAGATAGCCCTGCCTGTTCAGGCCGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGc                     >  1:1056333/1‑150 (MQ=255)
                                                                                                                                  ccAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCAACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCgg                   >  2:1264300/1‑150 (MQ=255)
                                                                                                                                       tAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAgcg              <  1:1924551/150‑1 (MQ=255)
                                                                                                                                        aacgctgcctgAACCGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCCGGCGCATCAAGCCACTTGTACGCCAGATGTTCAGTGAAAACGCTCTGCCGCTCGTGCGGCAgcgc             >  1:1787550/5‑150 (MQ=255)
                                                                                                                                        aacgctgcctgAACAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAgcgc             >  1:142035/5‑150 (MQ=255)
                                                                                                                                        aacgctgcctgAACAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCAATTGTAAGCCAGATGGTCAGTGAAAACGATCTGCCGCTCGTGCGGAAgagc             >  1:2014249/5‑150 (MQ=255)
                                                                                                                                         gCCCTGCCTGTTCAGGCAGTGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCa            >  1:1263582/1‑150 (MQ=255)
                                                                                                                                         gCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCa            <  2:827816/150‑1 (MQ=255)
                                                                                                                                          cccTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCGAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCaa           <  1:636182/150‑1 (MQ=255)
                                                                                                                                                ctgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAgc     >  1:967004/1‑149 (MQ=255)
                                                                                                                                                 tgTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGCTCTGCCGCTCGTGCGGAAGCGCAAGCCAGac    >  1:875985/1‑149 (MQ=255)
                                                                                                                                                  gTTCAGGCAGCGTTAATTACACACTGTTCCATCCCCTGCCGGTTGCGCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAAc   <  1:381644/150‑1 (MQ=255)
                                                                                                                                                  gTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGTCGCATCAAGCCACTAGTAAGCCAGAGGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAAc   >  1:1017457/1‑150 (MQ=255)
                                                                                                                                                  gTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAAc   >  1:155319/1‑150 (MQ=255)
                                                                                                                                                  gTTAAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGGCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAAc   <  2:364132/150‑1 (MQ=255)
                                                                                                                                                   ttCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAAcc  <  2:1517382/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACCGCCCAGCTTAGCCGCGGTTACCAGCTCACGAATGATTTTAGTGAAGATTTTACCGCGCTTAGCATCCTGCGCAGCTTTACGATGTCTGGTGTTGGCCCATTTACTATGACCTGCCATAAAAATATCTCCAGATAGCCCTGCCTGTTCAGGCAGCGTTAATTACAAACTGTTCAATCGCCTGCCGGTTGCTCCAGGACTTAGTGAGCGCCGCCGCAGCAGGCGCATCAAGCCACTTGTAAGCCAGATGTTCAGTGAAAACGATCTGCCGCTCGTGCGGAAGCGCAAGACAGAACC  >  NC_000913/1948032‑1948328

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: