Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,374,899 T→A 6.9% L11I (TTA→ATA) ‡ dcuD → putative transporter DcuD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,374,8990TA 6.9% 365.7 / 10.3 115L11I (TTA→ATA) ‡dcuDputative transporter DcuD
Reads supporting (aligned to +/- strand):  ref base T (45/62);  new base A (2/6);  total (47/68)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

Only 100 of 118 total aligned reads displayed.

CTACAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCG  >  NC_000913/3374753‑3375046
                                                                                                                                                  |                                                                                                                                                   
ctACAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATTGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTattaa                                                                                                                                                  <  1:1597781/150‑1 (MQ=255)
   cAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTattaatta                                                                                                                                               >  2:1316427/1‑150 (MQ=255)
    aaTTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTAc                                                                                                                                              <  1:1128877/150‑1 (MQ=255)
        aaaaaCACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATg                                                                                                                                          <  2:825068/150‑1 (MQ=255)
        aaaaaCACGACATCCACAACATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCTTTTATGGGATTGACCGTTTTCTTTTGACACGTAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATg                                                                                                                                          <  2:1457517/150‑1 (MQ=255)
         aaaaCACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATgg                                                                                                                                         <  1:156571/150‑1 (MQ=255)
         aaaaCACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATgg                                                                                                                                         <  2:506675/150‑1 (MQ=255)
          aaaCACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTAATCGTATTAATTACGATggg                                                                                                                                        >  1:152807/1‑150 (MQ=255)
           aaCACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGc                                                                                                                                       <  1:1132295/150‑1 (MQ=255)
           aaCACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTAATTAATACGATgaca                                                                                                                                       <  1:1078647/150‑4 (MQ=255)
           aaCACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTAATTAATACGATgaca                                                                                                                                       <  2:865075/150‑4 (MQ=255)
            acacGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCt                                                                                                                                      <  1:1523249/150‑1 (MQ=255)
            acacGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTAATTAATACGATgacag                                                                                                                                      <  2:1189776/150‑5 (MQ=255)
            acacGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTAATTAATACGATgacag                                                                                                                                      <  1:1158374/150‑5 (MQ=255)
             cacGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCTACAATGTTCGGCATAATTATATCTGTGATCGTAATTAATACGATgactgt                                                                                                                                     <  1:1229980/150‑6 (MQ=255)
             cacGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTAATTAATACGATgacag                                                                                                                                      <  2:328896/149‑5 (MQ=255)
             cacGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTAATTAATACGATgacag                                                                                                                                      >  1:328896/1‑145 (MQ=255)
                gACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTAttt                                                                                                                                  <  1:1264037/150‑1 (MQ=255)
                 aCATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTg                                                                                                                                 >  2:1056361/1‑150 (MQ=255)
                 aCATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTg                                                                                                                                 <  1:730718/150‑1 (MQ=255)
                     ccGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATcc                                                                                                                             <  2:124410/150‑1 (MQ=255)
                     ccGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATcc                                                                                                                             <  2:1571022/150‑1 (MQ=255)
                     ccGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATcc                                                                                                                             <  2:1534618/150‑1 (MQ=255)
                     ccGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATcc                                                                                                                             >  2:303193/1‑150 (MQ=255)
                      cGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCt                                                                                                                            >  2:270076/1‑150 (MQ=255)
                       gCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTg                                                                                                                           <  2:508253/150‑1 (MQ=255)
                        cACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGa                                                                                                                          <  2:306816/150‑1 (MQ=255)
                         aCCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGaa                                                                                                                         >  1:1173092/1‑150 (MQ=255)
                              aaaTAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCTGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACt                                                                                                                    <  1:1113967/150‑1 (MQ=255)
                              aaaTAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACt                                                                                                                    <  1:41855/150‑1 (MQ=255)
                              aaaTACCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCTTTTATTGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACt                                                                                                                    <  1:373887/150‑1 (MQ=255)
                               aaTATCCTTTAAAAAAATAAAACATAGATGTTTTAAATATCCAGTATATGGCATTGACAGTTTTCTATTCAAACGGAGTTCAACAATGTTCTGCATAATTATATCTTTAATACTATTAATTACTATGTGATATTTGATCCTGAAAAAATa                                                                                                                   >  2:561449/1‑150 (MQ=255)
                               aaTAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTa                                                                                                                   >  1:412811/1‑150 (MQ=255)
                               aaTAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGATCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTa                                                                                                                   >  1:910193/1‑150 (MQ=255)
                                aTAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTAc                                                                                                                  <  2:1623400/150‑1 (MQ=255)
                                aTAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTAc                                                                                                                  <  2:10611/150‑1 (MQ=255)
                                    ccTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAAc                                                                                                              >  1:1497169/1‑150 (MQ=255)
                                     cTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAAcc                                                                                                             <  1:896072/150‑1 (MQ=255)
                                        gCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCa                                                                                                          >  1:131074/1‑150 (MQ=255)
                                          aaaaaaTATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAgg                                                                                                        >  2:426090/1‑150 (MQ=255)
                                          aaaaaaTATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCATGAAAAACTACAAACCTCAgg                                                                                                        >  1:715130/1‑150 (MQ=255)
                                           aaaaaTATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAggt                                                                                                       >  1:1094678/1‑150 (MQ=255)
                                            aaaaTATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAAATACAAACCTCAggtg                                                                                                      <  2:150433/150‑1 (MQ=255)
                                              aaTATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAggtggt                                                                                                    <  1:934396/150‑1 (MQ=255)
                                              aaTATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAggtggt                                                                                                    <  1:940183/150‑1 (MQ=255)
                                                 ataACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgct                                                                                                 >  2:1326715/1‑150 (MQ=255)
                                                  taACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctg                                                                                                >  1:672779/1‑150 (MQ=255)
                                                    aCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctggc                                                                                              <  2:1559796/150‑1 (MQ=255)
                                                    aCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctggc                                                                                              <  1:1498152/150‑1 (MQ=255)
                                                    aCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctggc                                                                                              <  1:503689/150‑1 (MQ=255)
                                                    aCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctggc                                                                                              <  2:1285188/150‑1 (MQ=255)
                                                        cGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGcc                                                                                          <  2:1645436/150‑1 (MQ=255)
                                                         gttgttTTCAATCTGCCGTTTATGGGATTGCCCGTTTTATTTTGACACGGAGTTCAACAATGTTAGGCATAATTATATCTGTCATCGTATTAATTAGGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCg                                                                                         >  1:1515706/1‑150 (MQ=255)
                                                          ttggtTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGc                                                                                        <  1:1491582/150‑1 (MQ=255)
                                                           tgttTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCa                                                                                       <  2:317405/150‑1 (MQ=255)
                                                             ttttCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTGAACAATGTTCGGCATAATTATATCTGTCATCGTTTTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAgg                                                                                     <  1:1411501/150‑1 (MQ=255)
                                                             ttttCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTTAGGTGGTGCTGGCTGCCGCAgg                                                                                     <  1:1378501/150‑1 (MQ=255)
                                                               ttCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTa                                                                                   <  2:1497169/150‑1 (MQ=255)
                                                                tCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTAt                                                                                  >  2:859343/1‑150 (MQ=255)
                                                                tCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTAt                                                                                  <  1:270076/150‑1 (MQ=255)
                                                                       gCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTTTTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCtg                                                                           <  2:566441/150‑1 (MQ=255)
                                                                         cGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCACAGGTGGTGCTGGCTGACGCAGGTATCTTCCtgat                                                                         >  1:586215/1‑150 (MQ=255)
                                                                                 ggATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCg                                                                 >  1:1001114/1‑150 (MQ=255)
                                                                                  taTTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCgg                                                                <  1:674757/149‑1 (MQ=255)
                                                                                    ttGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGgtg                                                              >  2:1221585/1‑150 (MQ=255)
                                                                                        ccGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTg                                                          >  2:1187635/1‑150 (MQ=255)
                                                                                           ttttCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGtt                                                       >  2:224153/1‑150 (MQ=255)
                                                                                             ttCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAg                                                     <  1:1221585/150‑1 (MQ=255)
                                                                                             ttCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAg                                                     <  1:246946/150‑1 (MQ=255)
                                                                                               cTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAggg                                                   >  2:550381/1‑150 (MQ=255)
                                                                                                ttttGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGt                                                  >  1:526943/1‑150 (MQ=255)
                                                                                                 tttGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGtt                                                 >  2:1074708/1‑150 (MQ=255)
                                                                                                     acacGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCggt                                             >  2:1477406/1‑150 (MQ=255)
                                                                                                     acacGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCggt                                             >  1:841679/1‑150 (MQ=255)
                                                                                                         ggAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGgtggtg                                         <  2:1094678/150‑1 (MQ=255)
                                                                                                         ggAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGgtggtg                                         <  2:1629226/150‑1 (MQ=255)
                                                                                                           aGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTa                                       >  2:1385300/1‑150 (MQ=255)
                                                                                                           aGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTa                                       <  2:47624/150‑1 (MQ=255)
                                                                                                            tttCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTAc                                      <  1:1253046/149‑1 (MQ=255)
                                                                                                            gTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTAc                                      <  1:550381/150‑1 (MQ=255)
                                                                                                             ttCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACt                                     <  1:1029431/150‑1 (MQ=255)
                                                                                                               caacaaTGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCg                                   >  2:1121358/1‑150 (MQ=255)
                                                                                                               caacaaTGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCg                                   >  2:452334/1‑150 (MQ=255)
                                                                                                               caacaaTGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCg                                   <  1:334580/150‑1 (MQ=255)
                                                                                                                  caaTGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATc                                >  1:205391/1‑150 (MQ=255)
                                                                                                                   aaTGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATcc                               <  2:261761/150‑1 (MQ=255)
                                                                                                                    aTGTTCGGCATAATTATATCTGTCATCGTATTAATTAGGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATccc                              <  1:621783/150‑1 (MQ=255)
                                                                                                                     tGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCcca                             >  1:903186/1‑150 (MQ=255)
                                                                                                                            cATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAag                      >  2:1468023/1‑150 (MQ=255)
                                                                                                                              tAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAagca                    >  2:775077/1‑150 (MQ=255)
                                                                                                                                aTTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAagcagc                  >  2:1498417/1‑150 (MQ=255)
                                                                                                                                aTTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTAAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAagcagc                  >  2:1117525/1‑150 (MQ=255)
                                                                                                                                  aatatCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGGGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGGGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCgg                <  2:1095980/149‑1 (MQ=255)
                                                                                                                                       cTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGATAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACt           >  2:944489/1‑150 (MQ=255)
                                                                                                                                        gcccATCGTAATTAATACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACtt          >  2:1526091/4‑150 (MQ=255)
                                                                                                                                         gTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATccc                              <  1:380050/129‑1 (MQ=255)
                                                                                                                                         gTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATccc                              >  2:380050/1‑129 (MQ=255)
                                                                                                                                               gTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGAtcgtc   <  2:672779/150‑1 (MQ=255)
                                                                                                                                               gTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGAtcgtc   >  1:802739/1‑150 (MQ=255)
                                                                                                                                                tATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGAtcgtcg  <  1:1098670/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTACAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCG  >  NC_000913/3374753‑3375046

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: