Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 208,552 T→A 6.7% L1143I (TTA→ATA)  dnaE → DNA polymerase III subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913208,5520TA 6.7% 282.7 / 8.1 90L1143I (TTA→ATA) dnaEDNA polymerase III subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (43/41);  new base A (4/2);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

TTGATGACCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGTTAGCC  >  NC_000913/208403‑208699
                                                                                                                                                     |                                                                                                                                                   
ttGATGACCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTttat                                                                                                                                                     <  1:1005601/150‑1 (MQ=255)
 tgatgaCCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTttatt                                                                                                                                                    >  2:1343298/1‑150 (MQ=255)
 tgatgaCCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTttatt                                                                                                                                                    >  1:752779/1‑150 (MQ=255)
  gatgaCCAGCTTTTAAACCGACTCGGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTttatta                                                                                                                                                   >  2:649137/1‑150 (MQ=255)
  gatgaCCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTttatta                                                                                                                                                   <  2:175539/150‑1 (MQ=255)
  gatgaCCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTaa                                                                                                                                                  >  2:854736/1‑150 (MQ=255)
    tgaCCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTaaa                                                                                                                                                 >  2:1188147/1‑150 (MQ=255)
     gaCCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAAc                                                                                                                                                <  1:1029172/150‑1 (MQ=255)
       ccAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGa                                                                                                                                              >  1:783162/1‑150 (MQ=255)
         aGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGAtc                                                                                                                                            >  1:138406/1‑150 (MQ=255)
          gCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGAtct                                                                                                                                           >  1:438793/1‑150 (MQ=255)
          gCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTAATAAACGAtcg                                                                                                                                           <  1:365917/150‑2 (MQ=255)
          ccTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGGTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTACACGAtct                                                                                                                                           <  2:505809/149‑1 (MQ=255)
             tttAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCg                                                                                                                                        >  1:883661/1‑150 (MQ=255)
             tttAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCg                                                                                                                                        >  2:33659/1‑150 (MQ=255)
               tAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTg                                                                                                                                      >  1:303328/1‑150 (MQ=255)
                aaaCCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTgg                                                                                                                                     >  1:589188/1‑150 (MQ=255)
                aaaCCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTgg                                                                                                                                     <  1:649137/150‑1 (MQ=255)
                aaaCCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATCCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTgg                                                                                                                                     <  1:112287/150‑1 (MQ=255)
                aaaCCGACTCCGTCAGTCTCTGGAACCCCACCGCTATGGGACAATTCCAGTACATCACTACAATCAGAGCCCGGATGCACGCGCGCCCTTGCGTTTTGGCGCGACGTGCCGTGTCTCTACGAGCGATCGTTTATTTAACGATATCCGTgg                                                                                                                                     >  2:924552/1‑150 (MQ=255)
                  aCCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTgg                                                                                                                                     >  2:676843/1‑148 (MQ=255)
                  aCCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTgg                                                                                                                                     <  1:676843/148‑1 (MQ=255)
                         ccGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTgg                                                                                                                            >  2:1056314/1‑150 (MQ=255)
                               gTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGa                                                                                                                      >  1:933114/1‑150 (MQ=255)
                                tctcTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAg                                                                                                                     <  1:1035079/150‑1 (MQ=255)
                                 ctctGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGc                                                                                                                    <  2:1156633/150‑1 (MQ=255)
                                    tGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAgg                                                                                                                 >  1:832700/1‑150 (MQ=255)
                                     ggAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGt                                                                                                                >  2:1299136/1‑150 (MQ=255)
                                      gAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTg                                                                                                               >  2:558651/1‑150 (MQ=255)
                                            cACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTg                                                                                                         <  2:1378807/150‑1 (MQ=255)
                                             aCCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTgg                                                                                                        <  2:589188/150‑1 (MQ=255)
                                              ccGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGa                                                                                                                      <  2:1003419/135‑1 (MQ=255)
                                              ccGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGa                                                                                                                      >  1:1003419/1‑135 (MQ=255)
                                               cGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAg                                                                                                      <  1:1191594/150‑1 (MQ=255)
                                               cGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAg                                                                                                      <  2:359971/150‑1 (MQ=255)
                                                gCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGt                                                                                                     >  1:434933/1‑150 (MQ=255)
                                                   ctGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTg                                                                                                  <  1:1188147/150‑1 (MQ=255)
                                                         cAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTTGGAGCAGGTGGAACTGGAGTTTGACTAAt                                                                                            >  1:67365/1‑150 (MQ=255)
                                                           aTTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTa                                                                                              <  1:1223335/146‑1 (MQ=255)
                                                           aTTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTa                                                                                              >  2:1223335/1‑146 (MQ=255)
                                                             tCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAg                                                                                        <  1:925877/150‑1 (MQ=255)
                                                                    cATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACt                                                                                 >  2:1253383/1‑150 (MQ=255)
                                                                      tctcTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTAt                                                                               <  2:938526/150‑1 (MQ=255)
                                                                       ctctACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTagg                                                                              >  2:1205314/1‑148 (MQ=255)
                                                                       ctctACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATg                                                                              <  2:425106/150‑1 (MQ=255)
                                                                         ctactaTCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAg                                                                            >  1:338249/1‑150 (MQ=255)
                                                                         ctactaTCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAg                                                                            <  1:33659/150‑1 (MQ=255)
                                                                           actaTCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTc                                                                          >  2:1317982/1‑150 (MQ=255)
                                                                           actaTCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTc                                                                          >  1:414381/1‑150 (MQ=255)
                                                                                cAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAAt                                                                     >  1:440622/1‑150 (MQ=255)
                                                                                    gggCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTcc                                                                 >  2:386645/1‑150 (MQ=255)
                                                                                       cGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTg                                                              <  1:310580/150‑1 (MQ=255)
                                                                                       cGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTg                                                              <  2:379699/150‑1 (MQ=255)
                                                                                          aTGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGAtt                                                           <  2:550949/150‑1 (MQ=255)
                                                                                           tGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGAttt                                                          >  1:111126/1‑150 (MQ=255)
                                                                                           tGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGAttt                                                          >  1:492667/1‑150 (MQ=255)
                                                                                            gCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGAtttt                                                         <  1:1343298/150‑1 (MQ=255)
                                                                                             cACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTg                                                        <  2:1055072/150‑1 (MQ=255)
                                                                                             cACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTg                                                        <  2:4631/150‑1 (MQ=255)
                                                                                               cgcgcgcgGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGaa                                                      >  1:103620/1‑150 (MQ=255)
                                                                                                gcgcgcgGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAAc                                                     <  1:386645/150‑1 (MQ=255)
                                                                                                    gcgGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGcc                                                 <  1:558651/150‑1 (MQ=255)
                                                                                                          gCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGc                                           <  2:440622/150‑1 (MQ=255)
                                                                                                          gCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGc                                           <  2:783162/150‑1 (MQ=255)
                                                                                                          gCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGc                                           <  1:1180104/150‑1 (MQ=255)
                                                                                                            gTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTTGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCag                                         <  2:67365/150‑1 (MQ=255)
                                                                                                              tttGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCagag                                       <  1:1253383/150‑1 (MQ=255)
                                                                                                               ttGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGc                                      <  2:1473682/150‑1 (MQ=255)
                                                                                                                tGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCt                                     <  1:179882/150‑1 (MQ=255)
                                                                                                                  gcgcGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAATCGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTgg                                   <  2:214691/150‑1 (MQ=255)
                                                                                                                    gcgACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGGGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGaa                                 <  1:1315815/150‑1 (MQ=255)
                                                                                                                     cgACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAg                                <  2:745379/150‑1 (MQ=255)
                                                                                                                      gACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGc                               >  1:958727/1‑150 (MQ=255)
                                                                                                                          tGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGaaa                           >  1:199027/1‑150 (MQ=255)
                                                                                                                           ggCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGCGCTGGAAGCGaaaa                          >  1:1336093/1‑150 (MQ=255)
                                                                                                                             cGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATc                        >  1:965385/1‑150 (MQ=255)
                                                                                                                                 tctctcCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGAtt                    >  1:544162/1‑150 (MQ=255)
                                                                                                                                  ctctcCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATtc                   >  1:1252924/1‑150 (MQ=255)
                                                                                                                                   tctcCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATtct                  >  1:1459385/1‑150 (MQ=255)
                                                                                                                                   tctcCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATtct                  <  2:1354645/150‑1 (MQ=255)
                                                                                                                                    ctcCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATtctc                 >  2:1096673/1‑150 (MQ=255)
                                                                                                                                    ctcCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATtctc                 >  2:29678/1‑150 (MQ=255)
                                                                                                                                     tcCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTctct                <  2:526357/150‑1 (MQ=255)
                                                                                                                                     ccacggagATCGTTTAATAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTctct                >  1:1127193/8‑150 (MQ=255)
                                                                                                                                      cacggagATCGTTTAATAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGc                               >  1:641867/7‑134 (MQ=255)
                                                                                                                                      cacggagATCGTTTAATAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGc                               <  2:641867/128‑1 (MQ=255)
                                                                                                                                      cacggagATCGTTTAATAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTg               >  2:276310/7‑150 (MQ=255)
                                                                                                                                        gAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGAc             <  1:1317982/150‑1 (MQ=255)
                                                                                                                                         ggCGATCGTTTAATAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACt            >  1:1138788/2‑150 (MQ=255)
                                                                                                                                          gCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTg           >  1:315458/1‑150 (MQ=255)
                                                                                                                                             aTCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCgg        >  2:1208495/1‑150 (MQ=255)
                                                                                                                                              tCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGt       <  2:111126/150‑1 (MQ=255)
                                                                                                                                              tCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGGCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGt       <  2:883661/150‑1 (MQ=255)
                                                                                                                                               cGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGtt      >  1:1408794/1‑150 (MQ=255)
                                                                                                                                               cGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGtt      >  1:762618/1‑150 (MQ=255)
                                                                                                                                                   tattaAACGATCTCCGTGGCCTCATTGGTTCGGCGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGCCTGCGGTTAGcc  >  2:529858/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTGATGACCAGCTTTTAAACCGACTCCGTCAGTCTCTGGAACCCCACCGCTCTGGGACAATTCCAGTACATCTCTACTATCAGAGGGCGGATGCACGCGCGCGGTTGCGTTTTGGCGCGACGTGGCGTGTCTCTCCGAGCGATCGTTTATTAAACGATCTCCGTGGCCTCATTGGTTCGGAGCAGGTGGAACTGGAGTTTGACTAATACAGGAATACTATGAGTCTGAATTTCCTTGATTTTGAACAGCCGATTGCAGAGCTGGAAGCGAAAATCGATTCTCTGACTGCGGTTAGCC  >  NC_000913/208403‑208699

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