Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,390,518 A→T 5.3% intergenic (+69/+65) aaeR → / ← tldD transcriptional regulator for aaeXAB operon/putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,390,5180AT 5.3% 301.7 / 4.9 113intergenic (+69/+65)aaeR/tldDtranscriptional regulator for aaeXAB operon/putative peptidase
Reads supporting (aligned to +/- strand):  ref base A (57/48);  new base T (2/4);  total (59/54)
Fisher's exact test for biased strand distribution p-value = 4.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

Only 100 of 117 total aligned reads displayed.

TGATCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTT  >  NC_000913/3390370‑3390662
                                                                                                                                                    |                                                                                                                                                 
tGATCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCCCGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCCCCCCaa                                                                                                                                                   >  1:821679/1‑149 (MQ=255)
tGATCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCaa                                                                                                                                                   >  2:1369270/1‑149 (MQ=255)
tGATCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCACGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCaa                                                                                                                                                   <  1:1100940/149‑1 (MQ=255)
 gATCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCaaa                                                                                                                                                  >  1:1237011/1‑149 (MQ=255)
     aaCTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAgggg                                                                                                                                              >  1:1481118/1‑149 (MQ=255)
      aCTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGc                                                                                                                                             >  1:44411/1‑149 (MQ=255)
         cGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAg                                                                                                                                          <  1:106533/149‑1 (MQ=255)
         cGCTGACGGATTATTTTGTTGAGGGCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAg                                                                                                                                          <  1:529848/149‑1 (MQ=255)
          ggtGACTGATTATTTTGTTGAGGGCGGTAAATTGTGTCAGGAGATGCACGGGCGCGGTAAAGAGAATTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGTGAgg                                                                                                                                         <  2:1407431/147‑1 (MQ=255)
          gCTGACGGATTATTTTGTTGAGCTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAgg                                                                                                                                         <  1:258076/149‑1 (MQ=255)
           cTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAggg                                                                                                                                        >  1:995100/1‑149 (MQ=255)
            ttactgATTATTTTGTTGAGGTCTGTAAATTGTTTCAGGAGAGGCACGGGCGCGGGAAAGAGATGTAATTCATATTGTACTGTTACGTTGTACATACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAgggg                                                                                                                                       <  1:1485985/144‑1 (MQ=39)
               cGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGAcc                                                                                                                                    >  2:896673/1‑149 (MQ=255)
               cGGATCATTTTGTTGAGGTCGGTAAATTGTATCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCGCCCCTTTGGGGTGAgggttag                                                                                                                                    <  2:622457/149‑5 (MQ=255)
                ggATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCg                                                                                                                                   >  2:650301/1‑149 (MQ=255)
                      ttttGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTg                                                                                                                             >  2:1176177/1‑149 (MQ=255)
                          gttgAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACg                                                                                                                         >  1:341473/1‑149 (MQ=255)
                              aGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACAGTTCGTGCACGTAGa                                                                                                                     >  1:94016/1‑149 (MQ=255)
                                    gTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGatt                                                                                                               <  2:1545618/149‑1 (MQ=255)
                                     tAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGCCCGTTCGTGCACGGAGAAAGatta                                                                                                              >  2:315229/1‑149 (MQ=255)
                                     tAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGagta                                                                                                              >  1:1452346/1‑149 (MQ=255)
                                      aaaTTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGattaa                                                                                                             >  1:175838/1‑149 (MQ=255)
                                         ttGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCCACGGGTTCCCCTCACCCTAACCCTCACGCCAACGGGGCGAGGGGACCGTTCGTGCCCGTAGAAAGattaatta                                                                                                          >  1:887320/1‑149 (MQ=255)
                                          tctatccggAGCTGCACGCGCGCGGGAACGCGAATGACTTCATATTGTACTGTTACGTTGTACAACCCTGTGCCAAGGGGGTCCCCGCACCCTAACCCACACCCCTAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTAt                                                                                                         <  2:1446595/142‑1 (MQ=255)
                                           gTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAACGGGGCGAGGGGACCTTTCGTGCACGTAGAAAGATTAATTATc                                                                                                        <  1:554883/149‑1 (MQ=255)
                                           gTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATc                                                                                                        <  1:650301/149‑1 (MQ=255)
                                             ttCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCt                                                                                                      >  1:1233044/1‑149 (MQ=255)
                                               cAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTc                                                                                                    <  2:669712/149‑1 (MQ=255)
                                                aGGAGATTCACGGGCGCGGGATAGAGAAGTAATTCATATTGTACTGTTACGTTGTAGAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACGCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCt                                                                                                   <  1:289353/149‑1 (MQ=255)
                                                       gCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATa                                                                                            <  1:1033424/149‑1 (MQ=255)
                                                       gCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATa                                                                                            >  1:1306334/1‑149 (MQ=255)
                                                        cACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAg                                                                                           >  1:350083/1‑149 (MQ=255)
                                                        cACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCATCACCCTAACCCTCACCCCAAAGGGGAGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAg                                                                                           >  2:1425780/1‑149 (MQ=255)
                                                         aCGGGCGCGGGAAAGAGAAGTAATTCATCTTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGt                                                                                          <  1:1197851/149‑1 (MQ=255)
                                                         aCGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGt                                                                                          >  1:115644/1‑149 (MQ=255)
                                                         aCGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGt                                                                                          >  1:455842/1‑149 (MQ=255)
                                                            ggCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACCTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGAGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGaa                                                                                       <  1:356081/149‑1 (MQ=255)
                                                             gcgcGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGaaa                                                                                      <  2:145454/149‑1 (MQ=255)
                                                                cgGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAatta                                                                                   >  2:695337/1‑149 (MQ=255)
                                                                  ggAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAGTTATCGTTCTGATAATAGTGAAattatt                                                                                 >  2:1127655/1‑149 (MQ=255)
                                                                   gAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCCCCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAAttatta                                                                                >  2:1547592/1‑149 (MQ=255)
                                                                    aaaGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCACACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATAAATTATCCTTCTGAAAATAGTGAAATTATTAc                                                                               <  2:1189683/149‑1 (MQ=255)
                                                                       gagaAGTAATTCATACTGTACTGGTACGTTGTACACACCTGTGCCAACGGGATCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGGGCACGTAGAAAGATTAATTCTCCTTCTGAAAATAGTGAAATTCTTACGCa                                                                            <  1:1423674/149‑1 (MQ=38)
                                                                        agaAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGGGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAg                                                                           >  1:132134/1‑149 (MQ=255)
                                                                            gTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTAcc                                                                       <  1:321950/149‑1 (MQ=255)
                                                                                 tCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCaa                                                                  >  1:768220/1‑149 (MQ=255)
                                                                                 tCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCaa                                                                  <  1:973573/149‑1 (MQ=255)
                                                                                 tCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCACCTCACACTACCCCTCACCCCAAAGGGTCGAGGGGACCGTTCGTGATCGTAGAAATATTAATTATACTTATGAAAATAGTGAAAGTATTTCGCATTACCgagca                                                                  >  1:1013732/1‑145 (MQ=255)
                                                                                   atatTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACa                                                                <  1:360912/149‑1 (MQ=255)
                                                                                    tatTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAg                                                               >  2:276835/1‑149 (MQ=255)
                                                                                        ttACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAg                                                           <  2:1145722/148‑1 (MQ=255)
                                                                                         tACTGTTACGTTGTACAACCCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTAGCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAgg                                                          >  1:1072068/1‑149 (MQ=255)
                                                                                         tACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAgg                                                          >  2:285381/1‑149 (MQ=255)
                                                                                         tACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAgg                                                          <  1:970723/149‑1 (MQ=255)
                                                                                              ttACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATc                                                     <  2:1481118/149‑1 (MQ=255)
                                                                                              ttACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATc                                                     <  1:896673/149‑1 (MQ=255)
                                                                                              ttACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCGGAAAATAGTGAAATTATTAAGCAGTACCGCCAACAGTCAGGTTATc                                                     >  2:705488/1‑149 (MQ=255)
                                                                                              ttACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGATCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATc                                                     >  2:52526/1‑149 (MQ=255)
                                                                                               tACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCg                                                    <  1:157505/149‑1 (MQ=255)
                                                                                                aCGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGa                                                   <  1:1176177/149‑1 (MQ=255)
                                                                                                 cGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTTCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGAc                                                  <  1:873513/149‑1 (MQ=255)
                                                                                                 cGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGAc                                                  >  2:944769/1‑149 (MQ=255)
                                                                                                  gTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACt                                                 >  1:1529794/1‑149 (MQ=255)
                                                                                                   ttGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCTCGGAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCTACAGTCAGGTTATCGACtt                                                >  1:1276044/1‑149 (MQ=255)
                                                                                                   ttGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGCCtt                                                >  1:1054068/1‑149 (MQ=255)
                                                                                                   ttGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACtt                                                >  2:823441/1‑149 (MQ=255)
                                                                                                     gTACAAACCTGTGCCAACGGGTTCCCCTCACCCTACCCCTCACCCCAAAGGGGCGAGGGGACCGGTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATATGGAAAGTATTACTCAGTACCGCCAACAGTCAGGTTATCGACTTTc                                              >  2:1062631/1‑149 (MQ=39)
                                                                                                      tACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGCCCGTTCGTGCACGTAGAAAGAGTAATTATCCTTCTGAAAATTGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCa                                             >  2:549244/1‑149 (MQ=255)
                                                                                                       aCAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCaa                                            >  2:1480134/1‑149 (MQ=255)
                                                                                                           aCCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGtt                                        <  1:1004931/149‑1 (MQ=255)
                                                                                                           aCCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGtt                                        <  1:818208/149‑1 (MQ=255)
                                                                                                            ccTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTtg                                       >  1:1405431/1‑149 (MQ=255)
                                                                                                            ccTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCATTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTtg                                       >  1:660331/1‑149 (MQ=255)
                                                                                                             cTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCCCGTAGAAAGATTAATTATCCTTCTGAAAAAAGTGAAATTATTACGCAGTACCGCCAACAGACAGGTTATCGACTTTCAACGTtgg                                      >  1:896466/1‑149 (MQ=255)
                                                                                                             cTGTGCCAACGGGTTCCCCTCACCCTAACCCTAACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTtgg                                      <  2:94016/149‑1 (MQ=255)
                                                                                                                 gCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTtggctgg                                 <  2:1237011/149‑1 (MQ=255)
                                                                                                                  ccACCGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCACCAGTCAGGTTATCGACGTTCAACGTttgctggc                                >  2:94706/1‑149 (MQ=255)
                                                                                                                  ccAACGGGTTCCCCTCACCCTACCCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTtggctggc                                <  1:61680/149‑1 (MQ=255)
                                                                                                                  ccAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTtggctggc                                <  2:150178/149‑1 (MQ=255)
                                                                                                                  ccAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGCGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTtggctggc                                <  2:1471362/149‑1 (MQ=255)
                                                                                                                   cAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGcc                               >  1:837474/1‑149 (MQ=255)
                                                                                                                          ttCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCATGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCa                        >  1:371565/1‑149 (MQ=255)
                                                                                                                            ccccTCACCCTCACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAAc                      <  2:1124638/149‑1 (MQ=255)
                                                                                                                              ccTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAATATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCGTGGGCCAGCCAAcgt                    <  2:341473/149‑3 (MQ=38)
                                                                                                                              ccTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAAGTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGCCTTTCAACGTTGGATGGCCCACGCCCACCg                    >  2:1014064/1‑149 (MQ=255)
                                                                                                                               cTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCgg                   <  2:44411/149‑1 (MQ=255)
                                                                                                                               cTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCgg                   <  2:115644/149‑1 (MQ=255)
                                                                                                                               cTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCgg                   >  1:1263756/1‑149 (MQ=255)
                                                                                                                                  aCCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCACCCGCCaa                >  2:957913/1‑149 (MQ=255)
                                                                                                                                    ccTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCAAGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAAc              >  2:966670/1‑149 (MQ=255)
                                                                                                                                     ggtccccTCGCCCCTTTGGGGTGAGGGTACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACt             >  1:365511/5‑149 (MQ=255)
                                                                                                                                       aaCCCTCACCCCACAGGGGCGAGGGGACCGTGCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGGGAAATTATTACGCAGTACCGCCAACAGTCAGGTAATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACttt           <  1:39773/149‑1 (MQ=255)
                                                                                                                                        ccccTCGCCCCTTTGGGGTGAGGGGACCGTTGGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGTAATTATTAGGCAGTACCGCCAACAGTCAGGATAGCGCGTTTCAACGTTGGCTGGCCCACGCCCACCGGCaaccgtag          >  1:1272586/2‑143 (MQ=255)
                                                                                                                                        aCCCTCCCCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGGGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTg          <  1:754983/149‑1 (MQ=255)
                                                                                                                                        aCCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTg          <  1:1562001/149‑1 (MQ=255)
                                                                                                                                        aCCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTg          <  2:350083/149‑1 (MQ=255)
                                                                                                                                         cccTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCATGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGc         <  2:371565/149‑1 (MQ=255)
                                                                                                                                          ccTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGcc        <  2:768220/149‑1 (MQ=255)
                                                                                                                                           cTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGccc       >  2:45186/1‑149 (MQ=255)
                                                                                                                                                cccAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCcttctt  >  1:1001876/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
TGATCAACTCGCTGACGGATTATTTTGTTGAGGTCGGTAAATTGTTTCAGGAGATGCACGGGCGCGGGAAAGAGAAGTAATTCATATTGTACTGTTACGTTGTACAAACCTGTGCCAACGGGTTCCCCTCACCCTAACCCTCACCCCAAAGGGGCGAGGGGACCGTTCGTGCACGTAGAAAGATTAATTATCCTTCTGAAAATAGTGAAATTATTACGCAGTACCGCCAACAGTCAGGTTATCGACTTTCAACGTTGGCTGGCCCACGCCAACCGGCAAACTTTGCCCTTCTT  >  NC_000913/3390370‑3390662

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