Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,374,901 A→T 7.0% L11F (TTA→TTT) ‡ dcuD → putative transporter DcuD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,374,9010AT 7.0% 330.5 / 12.4 115L11F (TTA→TTT) ‡dcuDputative transporter DcuD
Reads supporting (aligned to +/- strand):  ref base A (45/62);  new base T (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 = 6.87e-01

Only 100 of 122 total aligned reads displayed.

CTACAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATAT  >  NC_000913/3374753‑3375050
                                                                                                                                                    |                                                                                                                                                     
ctACAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATTGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTattaa                                                                                                                                                      <  1:1597781/150‑1 (MQ=255)
  aCAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTattaatt                                                                                                                                                    >  1:47624/1‑150 (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                                                                                                                                          <  1:1158374/150‑5 (MQ=255)
            acacGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTAATTAATACGATgacag                                                                                                                                          <  2:1189776/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)
                 aCATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTg                                                                                                                                     <  1:730718/150‑1 (MQ=255)
                     ccGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATcc                                                                                                                                 <  2:1571022/150‑1 (MQ=255)
                     ccGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATcc                                                                                                                                 <  2:1534618/150‑1 (MQ=255)
                      cGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCt                                                                                                                                >  2:270076/1‑150 (MQ=255)
                        cACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGa                                                                                                                              <  2:306816/150‑1 (MQ=255)
                         aCCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGaa                                                                                                                             >  1:1173092/1‑150 (MQ=255)
                              aaaTAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACt                                                                                                                        <  1:41855/150‑1 (MQ=255)
                              aaaTACCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCTTTTATTGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACt                                                                                                                        <  1:373887/150‑1 (MQ=255)
                               aaTATCCTTTAAAAAAATAAAACATAGATGTTTTAAATATCCAGTATATGGCATTGACAGTTTTCTATTCAAACGGAGTTCAACAATGTTCTGCATAATTATATCTTTAATACTATTAATTACTATGTGATATTTGATCCTGAAAAAATa                                                                                                                       >  2:561449/1‑150 (MQ=255)
                               aaTAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGATCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTa                                                                                                                       >  1:910193/1‑150 (MQ=255)
                                aTAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTAc                                                                                                                      <  2:10611/150‑1 (MQ=255)
                                aTAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTAc                                                                                                                      <  2:1623400/150‑1 (MQ=255)
                                aTAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTAc                                                                                                                      <  1:636983/150‑1 (MQ=255)
                                     cTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAAcc                                                                                                                 <  1:896072/150‑1 (MQ=255)
                                       tGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTc                                                                                                               <  2:1473497/150‑1 (MQ=255)
                                        gCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCa                                                                                                              >  1:131074/1‑150 (MQ=255)
                                         cAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAg                                                                                                             >  2:1029431/1‑150 (MQ=255)
                                           aaaaaTATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAggt                                                                                                           >  1:1094678/1‑150 (MQ=255)
                                            aaaaTATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAAATACAAACCTCAggtg                                                                                                          <  2:150433/150‑1 (MQ=255)
                                              aaTATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAggtggt                                                                                                        <  1:940183/150‑1 (MQ=255)
                                              aaTATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAggtggt                                                                                                        <  1:934396/150‑1 (MQ=255)
                                                 ataACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgct                                                                                                     >  2:1326715/1‑150 (MQ=255)
                                                  taACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctg                                                                                                    >  1:672779/1‑150 (MQ=255)
                                                    aCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctggc                                                                                                  <  1:1498152/150‑1 (MQ=255)
                                                    aCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctggc                                                                                                  <  2:1285188/150‑1 (MQ=255)
                                                    aCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctggc                                                                                                  <  2:1559796/150‑1 (MQ=255)
                                                    aCATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctggc                                                                                                  <  1:503689/150‑1 (MQ=255)
                                                      aTCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTgctggctg                                                                                                <  1:87933/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)
                                                                                 ggATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCg                                                                     >  1:1001114/1‑150 (MQ=255)
                                                                                  taTTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCgg                                                                    <  1:674757/149‑1 (MQ=255)
                                                                                   aTTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGgt                                                                   <  2:131074/150‑1 (MQ=255)
                                                                                    ttGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGgtg                                                                  >  2:1221585/1‑150 (MQ=255)
                                                                                      gACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTc                                                                >  1:1250089/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:325026/1‑150 (MQ=255)
                                                                                               cTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAggg                                                       >  2:550381/1‑150 (MQ=255)
                                                                                               cTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAggg                                                       <  1:1412552/150‑1 (MQ=255)
                                                                                                ttttGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGt                                                      >  1:526943/1‑150 (MQ=255)
                                                                                                 tttGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGtt                                                     >  2:1074708/1‑150 (MQ=255)
                                                                                                    gACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCgg                                                  >  2:481795/1‑150 (MQ=255)
                                                                                                     acacGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCggt                                                 >  2:1477406/1‑150 (MQ=255)
                                                                                                     acacGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCggt                                                 <  2:715130/150‑1 (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: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                                       <  1:334580/150‑1 (MQ=255)
                                                                                                               caacaaTGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCg                                       >  2:452334/1‑150 (MQ=255)
                                                                                                               caacaaTGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCg                                       >  2:1121358/1‑150 (MQ=255)
                                                                                                                  caaTGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATc                                    >  1:205391/1‑150 (MQ=255)
                                                                                                                   aaTGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATcc                                   <  2:261761/150‑1 (MQ=255)
                                                                                                                    aTGTTCGGCATAATTATATCTGTCATCGTATTAATTAGGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATccc                                  <  1:621783/150‑1 (MQ=255)
                                                                                                                            cATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAag                          >  2:1468023/1‑150 (MQ=255)
                                                                                                                                aTTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTAAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAagcagc                      >  2:1117525/1‑150 (MQ=255)
                                                                                                                                  aatatCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGGGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGGGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCgg                    <  2:1095980/149‑1 (MQ=255)
                                                                                                                                      tCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTAc                <  1:1121358/150‑1 (MQ=255)
                                                                                                                                       cTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGATAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACt               >  2:944489/1‑150 (MQ=255)
                                                                                                                                        gcccATCGTAATTAATACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACtt              >  2:1526091/4‑150 (MQ=255)
                                                                                                                                         gTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATccc                                  >  2:380050/1‑129 (MQ=255)
                                                                                                                                         gTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATccc                                  <  1:380050/129‑1 (MQ=255)
                                                                                                                                               gTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGAtcgtc       >  1:802739/1‑150 (MQ=255)
                                                                                                                                               gTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGAtcgtc       <  2:672779/150‑1 (MQ=255)
                                                                                                                                                tATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGAtcgtcg      <  1:1098670/150‑1 (MQ=255)
                                                                                                                                                   taattaCGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGa     <  2:122105/148‑1 (MQ=255)
                                                                                                                                                    aattaCGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGatat  >  2:281637/1‑150 (MQ=255)
                                                                                                                                                    aattaAGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCATGTATCTTCCTGATGATTTGCGTTGTCTTGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGatat  <  2:152807/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTACAATTAAAAACACGACATCCGCACCATAAATAGCCTTGCAAAAAATATAACATCGTTGTTTTCAATCTGCCGTTTATGGGATTGACCGTTTTCTTTTGACACGGAGTTCAACAATGTTCGGCATAATTATATCTGTCATCGTATTAATTACGATGGGCTATTTGATCCTGAAAAACTACAAACCTCAGGTGGTGCTGGCTGCCGCAGGTATCTTCCTGATGATGTGCGGTGTCTGGTTAGGGTTCGGTGGTGTACTCGATCCCACCAAAAGCAGCGGCTACTTGATCGTCGATAT  >  NC_000913/3374753‑3375050

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: