Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,397,970 A→C 7.7% Y137D (TAT→GAT) ‡ yhdE ← dTTP/UTP pyrophosphatase; m(5)UTP/m(5)CTP/pseudo‑UTP pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,397,9700AC 7.7% 361.3 / 17.4 126Y137D (TAT→GAT) ‡yhdEdTTP/UTP pyrophosphatase; m(5)UTP/m(5)CTP/pseudo‑UTP pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base A (55/61);  new base C (5/5);  total (60/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-01

Only 100 of 129 total aligned reads displayed.

GTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTGCGCAACATCTGCGCCGCATGCT  >  NC_000913/3397821‑3398119
                                                                                                                                                     |                                                                                                                                                     
gTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATa                                                                                                                                                       >  2:1434898/1‑150 (MQ=255)
  tAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCCCCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGc                                                                                                                                                     >  2:1519262/1‑150 (MQ=255)
   aaaaTTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACCGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGcc                                                                                                                                                    >  2:1460008/1‑150 (MQ=255)
    aaaTTACTTAATAAATCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGccc                                                                                                                                                   <  2:1311970/150‑1 (MQ=255)
     aaTAACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGAATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCcg                                                                                                                                                  <  2:1195975/150‑1 (MQ=255)
       ttACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCcgcg                                                                                                                                                <  2:1427633/150‑1 (MQ=255)
         aCTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCGGGCTATGTCGCCAGCGAt                                                                                                                                              <  2:1218658/150‑1 (MQ=255)
          attAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCGGGCTATGTCGCCAGCGATg                                                                                                                                             <  1:1586601/149‑1 (MQ=255)
            taataaTTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTc                                                                                                                                           >  1:1681159/1‑150 (MQ=255)
            taataaTTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTc                                                                                                                                           >  1:148292/1‑150 (MQ=255)
             aataatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCGGGCTATGTCGCCAGCGATgaac                                                                                                                                          <  1:1376971/150‑4 (MQ=255)
              ataatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACACAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCtt                                                                                                                                         <  1:185634/150‑1 (MQ=255)
               taatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATCGCGGGCTATGTcgcc                                                                                                                                        <  1:755520/150‑4 (MQ=255)
               taatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATCGCGGGCTATGTcgcc                                                                                                                                        <  1:609936/150‑4 (MQ=255)
               taatTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTc                                                                                                                                        >  1:1616256/1‑150 (MQ=255)
                  ttCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTc                                                                                                                                     <  2:1000088/150‑1 (MQ=255)
                    cATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTg                                                                                                                                   >  1:1243752/1‑150 (MQ=255)
                    cATACGTTTCAACCAGCGGTAAGCCGACTACGGCGAGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCTCTGGCGAAATAGCCCGCGATGTCTTCGTCTg                                                                                                                                   >  1:709125/1‑150 (MQ=255)
                       aCGTTTCCACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCGGCCAACACAGCCACCCAGCCCCTGAATACCGGATGCAGCTGCGTTATCTATCGGTTCATCGGTGGCGCCATCGCCCGGGATGTCGTCGTCTGttg                                                                                                                                >  2:1126162/1‑149 (MQ=255)
                          tttCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACg                                                                                                                             <  2:750302/150‑1 (MQ=255)
                             cAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTc                                                                                                                          <  1:1001675/150‑1 (MQ=255)
                              taCCAGCGGTAAGGCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTATCGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCt                                                                                                                         <  1:815070/149‑1 (MQ=255)
                              aaCCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTaa                                                                                                                               <  2:1229371/144‑1 (MQ=255)
                              aaCCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCt                                                                                                                         <  2:216572/150‑1 (MQ=255)
                              aaCCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCt                                                                                                                         <  1:332255/150‑1 (MQ=255)
                               accaccGGTAAGCCGACTACGGCGTGATACCTGCAATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACATGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGAGATGTCTTCGTCTCTTAACGttata                                                                                                                        >  2:815070/1‑147 (MQ=255)
                                 cAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGaa                                                                                                                      <  1:138012/150‑1 (MQ=255)
                                     ggTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTc                                                                                                                  <  1:1157718/150‑1 (MQ=255)
                                          gCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATc                                                                                                             >  1:1466782/1‑150 (MQ=255)
                                          gCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATc                                                                                                             >  1:557450/1‑150 (MQ=255)
                                               cTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGa                                                                                                        <  1:1550272/150‑1 (MQ=255)
                                                tACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGac                                                                                                       <  1:494654/150‑1 (MQ=255)
                                                tACGGAGTGATAGCTGCCATATATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGac                                                                                                       <  2:1584989/150‑1 (MQ=255)
                                                 aCGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGacc                                                                                                      >  1:63551/1‑150 (MQ=255)
                                                 aCGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAAAACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGATAACGTGCGGAAAGTCACATCGGTGccc                                                                                                      >  2:97711/1‑150 (MQ=255)
                                                  cGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGacca                                                                                                     >  1:521759/1‑150 (MQ=255)
                                                    gCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGaccacc                                                                                                   >  2:54508/1‑150 (MQ=255)
                                                       tGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAgg                                                                                                <  1:694482/150‑1 (MQ=255)
                                                         aTAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCa                                                                                              >  1:1000364/1‑150 (MQ=255)
                                                            gCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATc                                                                                           <  2:353685/150‑1 (MQ=255)
                                                            gCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATc                                                                                           <  2:1616256/150‑1 (MQ=255)
                                                             cTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGCCCCCCAGGCAATCg                                                                                          >  1:1629581/1‑150 (MQ=255)
                                                             cTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCg                                                                                          <  2:1310624/150‑1 (MQ=255)
                                                             cTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCg                                                                                          <  1:504228/150‑1 (MQ=255)
                                                               gCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCgag                                                                                        <  1:120207/150‑1 (MQ=255)
                                                                ccATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCgaga                                                                                       >  1:956776/1‑150 (MQ=255)
                                                                ccATTTATCTTCCTGACAAAACAGCCACCCAGCCCCGGAATACCGTATGCACCTGCTTTAGCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCgaga                                                                                       <  1:97711/150‑1 (MQ=255)
                                                                 cATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGaa                                                                                      <  1:1206295/150‑1 (MQ=255)
                                                                       tCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTgct                                                                                <  1:133619/150‑1 (MQ=255)
                                                                        cTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTgctg                                                                               <  1:1453724/150‑1 (MQ=255)
                                                                         ttCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTgctgg                                                                              <  2:148292/150‑1 (MQ=255)
                                                                         ttCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTgctgg                                                                              >  2:465466/1‑150 (MQ=255)
                                                                          tCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCTTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTgctggc                                                                             >  2:677686/1‑150 (MQ=255)
                                                                           ccTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTgctggct                                                                            >  2:646388/1‑150 (MQ=255)
                                                                            cTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAAtgtg                                                                                  <  1:394002/143‑1 (MQ=255)
                                                                            cTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAAtgtg                                                                                  >  2:394002/1‑143 (MQ=255)
                                                                            cTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGATCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTgctggctg                                                                           <  1:1255610/150‑1 (MQ=255)
                                                                                cAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCgg                                                                       >  2:1544157/1‑150 (MQ=255)
                                                                                 aaaaCAGCCACCCAGCCCCTGAATACAGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCTGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGc                                                                      <  1:1224415/150‑1 (MQ=255)
                                                                                  aaaCAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGcc                                                                     >  2:273022/1‑150 (MQ=255)
                                                                                    aCAGCCACCCAGCCCCTGAATACCGTATGCACATGCTTTATCTAACGGTTCATAGCTGGCGACATAGCAAGCGATGTCTTCGACTGTTAACGTTCTGAAAGTAACATCAGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCaa                                                                   >  2:1700127/1‑150 (MQ=255)
                                                                                       gCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGc                                                                <  2:521759/150‑1 (MQ=255)
                                                                                        ccACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGcc                                                               <  1:674597/150‑1 (MQ=255)
                                                                                         cACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCa                                                              <  2:907175/150‑1 (MQ=255)
                                                                                          aCCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCAc                                                             >  1:813452/1‑150 (MQ=255)
                                                                                            ccAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCActg                                                           >  1:1369990/1‑150 (MQ=255)
                                                                                             cAGCCCCTGAATACCGTATGCACCTGCTTTGTCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGAACACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCACAGCCActgc                                                          >  1:206025/1‑150 (MQ=255)
                                                                                              aGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCActgct                                                         <  2:557450/150‑1 (MQ=255)
                                                                                                ccccTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGt                                                       >  1:690720/1‑150 (MQ=255)
                                                                                                 cccTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGAGGTCTTCGTCGGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGtg                                                      >  2:94269/1‑149 (MQ=255)
                                                                                                  ccTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGtca                                                     >  2:126727/1‑150 (MQ=255)
                                                                                                  ccTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGtca                                                     >  2:1163031/1‑150 (MQ=255)
                                                                                                    tGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTGACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGtcatc                                                   <  2:1117487/150‑1 (MQ=255)
                                                                                                       aTACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCAcc                                                >  1:420885/1‑150 (MQ=255)
                                                                                                        tACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCCACGCCACTGCTGTCATCACCt                                               >  2:843028/1‑150 (MQ=255)
                                                                                                          ccGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGa                                             >  2:1160517/1‑150 (MQ=255)
                                                                                                             tATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATgg                                          <  2:1224195/150‑1 (MQ=255)
                                                                                                              aTGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATggg                                         <  1:1544157/150‑1 (MQ=255)
                                                                                                              aTGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATggg                                         <  1:440169/150‑1 (MQ=255)
                                                                                                               tGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGt                                        >  2:1514227/1‑150 (MQ=255)
                                                                                                                 cACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCt                                      >  1:323611/1‑150 (MQ=255)
                                                                                                                   ccTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGa                                    >  1:1309507/1‑150 (MQ=255)
                                                                                                                    cTGCTTTATCTAACGGTGCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGAc                                   <  2:956776/150‑1 (MQ=255)
                                                                                                                     tGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGATCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGAcc                                  >  1:682618/1‑150 (MQ=255)
                                                                                                                      gCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGAccc                                 <  1:1163031/150‑1 (MQ=255)
                                                                                                                         ttATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGCCCCGAt                              >  2:1627548/1‑150 (MQ=255)
                                                                                                                          tATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATa                             >  1:364452/1‑150 (MQ=255)
                                                                                                                           aTCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATaa                            >  2:119803/1‑150 (MQ=255)
                                                                                                                             cTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCAGATGGGTCTGACCCGATAAtt                          >  1:1460139/1‑150 (MQ=255)
                                                                                                                                 gaagacATCGCGGGCTATGTCGCCAGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTgcg                      >  2:419220/6‑150 (MQ=255)
                                                                                                                                  ggTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTgcgc                     <  2:690720/150‑1 (MQ=255)
                                                                                                                                   gTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTGCGCa                    >  2:1034337/1‑150 (MQ=255)
                                                                                                                                   gTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTGCGCa                    >  1:705898/1‑150 (MQ=255)
                                                                                                                                     acATCGCGGGCTATGTCGCCAGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTGCGCAAc                  >  2:1125716/2‑150 (MQ=255)
                                                                                                                                      cATCGCGGGCTATGTCGCCAGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGCCCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTGCGCAACa                 >  1:551116/1‑150 (MQ=255)
                                                                                                                                           cTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTGCGCAACATCTgc            >  1:477859/1‑150 (MQ=255)
                                                                                                                                                 aCATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATAGCCTGATGGGTCTGACCGGTGAATTTGCGCAGCATCTgcgcgggg      >  2:1630066/1‑146 (MQ=255)
                                                                                                                                                    tAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTGCGCAACATCTGCGCCGCATGc   >  1:1546549/1‑150 (MQ=255)
                                                                                                                                                    tAGCCCGCGATGTCTTCGGCTGTTAACGTTCTGACAGTCACATCGGTGACCACCAGGCACTCGAGACTGTGCTGGCTGTCGGCCAACGCCACTGCTGGCATCACCTGATGGGTCTGACCCGATAATTTGCGCAACATCTGCGCCTCATGc   <  1:243415/150‑1 (MQ=255)
                                                                                                                                                     aGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCCCTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTGCGCAACATCTGCGCCGCATGct  >  1:610051/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTAAAATTACTTAATAATTCATACGTTTCAACCAGCGGTAAGCCGACTACGGCGTGATAGCTGCCATTTATCTTCCTGACAAAACAGCCACCCAGCCCCTGAATACCGTATGCACCTGCTTTATCTAACGGTTCATCGCTGGCGACATAGCCCGCGATGTCTTCGTCTGTTAACGTTCTGAAAGTCACATCGGTGACCACCAGGCAATCGAGAATGTGCTGGCTGTCGGCCAACGCCACTGCTGTCATCACCTGATGGGTCTGACCCGATAATTTGCGCAACATCTGCGCCGCATGCT  >  NC_000913/3397821‑3398119

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