Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,752,058 A→C I23L (ATC→CTC)  yiaU → putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,752,0580AC99.6% 816.6 / ‑5.9 234I23L (ATC→CTC) yiaUputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (83/150);  total (84/150)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

Only 100 of 236 total aligned reads displayed.

ATATAAAATTATATATAATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACC  >  NC_000913/3751961‑3752154
                                                                                                 |                                                                                                 
 tataAATTTATATATAATTGGGGAGTTTTAATGATATAATTACATCTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAGTCTCa                                                                                                <  2:2605695/100‑1 (MQ=255)
 tataAAATTATATATAATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCa                                                                                                <  1:2797386/100‑1 (MQ=255)
 tataAAATTATATATAATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCa                                                                                                <  1:1733827/100‑1 (MQ=255)
   taAAATTATATATAATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGc                                                                                              <  1:2992935/100‑1 (MQ=255)
   taAAATTATATATAATTAGGCTGTTGCGTTGAAAAAATTACAAATTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGc                                                                                              <  2:3516422/100‑1 (MQ=255)
    aaaaTTATATATAATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGcc                                                                                             <  1:541297/100‑1 (MQ=255)
    aaaaTTATATATAATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGcc                                                                                             <  1:4242669/100‑1 (MQ=255)
    aaaaTTATATATAATTGGGCTGTTGCGATGAAAAAATTACAACTTAAGTAAAAAAAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGcc                                                                                             <  2:3138523/100‑1 (MQ=255)
     aaaTTATATAGAATTGGGCTGTTGCTATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCa                                                                                            <  2:4813826/100‑1 (MQ=255)
      caTTACATATAAGTAGGCTGTTGCGATGAAAAAAATACAAATTAATTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCACTGAAAATCTCAGCCAt                                                                                           <  2:5124835/99‑1 (MQ=255)
      aaTTATATATAATTGGGATGTTGAGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCAt                                                                                           <  2:3368810/100‑1 (MQ=255)
       aTTATATATAATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATg                                                                                          >  1:2192791/1‑100 (MQ=255)
            atataATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCcgcg                                                                                     <  1:874108/100‑1 (MQ=255)
            atataATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCcgcg                                                                                     >  1:2609832/1‑100 (MQ=255)
            atataATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCcgcg                                                                                     >  1:274609/1‑100 (MQ=255)
             tataATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGa                                                                                    >  2:1744495/1‑100 (MQ=255)
              ataATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGAc                                                                                   <  1:929946/100‑1 (MQ=255)
               aaaTTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACt                                                                                  <  2:1352655/99‑1 (MQ=255)
                aaTTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTg                                                                                 <  1:3275153/100‑1 (MQ=255)
                aaTTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTg                                                                                 >  1:1071840/1‑100 (MQ=255)
                 aTTGGGCTGTTGCAATGACAAAAATACAACTTAAATACAGAAAATTAAAAATTATCTCGGTAATCGCTGCCAGTAAAAATCTAAGACATGACGCGACTGt                                                                                <  2:204446/100‑1 (MQ=255)
                  ttGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTa                                                                               <  1:767912/100‑1 (MQ=255)
                   tGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTAc                                                                              >  2:4098440/1‑100 (MQ=255)
                   tGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTAc                                                                              <  2:2725485/100‑1 (MQ=255)
                   tGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACCGTAc                                                                              <  1:3874212/100‑1 (MQ=255)
                     ggCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACtt                                                                            >  1:1119933/1‑100 (MQ=255)
                      gCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTg                                                                           >  2:2558114/1‑100 (MQ=255)
                       cTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTgg                                                                          <  1:1299130/100‑1 (MQ=255)
                         gTTCCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCa                                                                        >  2:5054576/1‑100 (MQ=255)
                          ttGCGATGACAAAAATACAACTTACGTACCGGGAGATAAATATTATCTCGGTAACCGCCGCCAGTGAAAACCTCAGCCATGCCGCGACTGTACTTGGCAt                                                                       <  2:3485219/100‑1 (MQ=255)
                           tGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATc                                                                      <  1:3423061/100‑1 (MQ=255)
                           tGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATc                                                                      <  2:4310397/100‑1 (MQ=255)
                            gCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCg                                                                     >  1:381364/1‑100 (MQ=255)
                              gATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCACTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGca                                                                   >  2:5402925/1‑100 (MQ=255)
                                tGAGAAAATTACAACTTAAGTACCGGGAGTAAAAAATTATCTCGGCAATCTCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCACCCcaca                                                                 <  2:5181688/100‑1 (MQ=255)
                                   cAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGc                                                              >  1:1106232/1‑100 (MQ=255)
                                   cAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGc                                                              >  2:2274917/1‑100 (MQ=255)
                                    aaaaTTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGcc                                                             <  1:2487705/100‑1 (MQ=255)
                                        ttACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACg                                                         <  2:2151786/100‑1 (MQ=255)
                                         tACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGt                                                        >  1:2847632/1‑100 (MQ=255)
                                           cAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCa                                                      <  1:1886549/100‑1 (MQ=255)
                                           cAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAATGAAAATCTCAGCCAGGCCGCGACTGTACTTAGCATCGCACAGGCCAACGTCa                                                      >  1:5185805/1‑100 (MQ=255)
                                             aCTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGc                                                    >  2:3688378/1‑100 (MQ=255)
                                             aCTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGc                                                    >  1:4387736/1‑100 (MQ=255)
                                             aCTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGc                                                    >  2:3005964/1‑100 (MQ=255)
                                             aCTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGc                                                    <  1:2920243/100‑1 (MQ=255)
                                              cTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCTATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCa                                                   >  1:3086432/1‑100 (MQ=255)
                                              cTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCa                                                   >  2:1214471/1‑100 (MQ=255)
                                                tAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCaaa                                                >  1:4778940/1‑100 (MQ=255)
                                                  aGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAata                                              <  1:4525180/100‑1 (MQ=255)
                                                     aCCGGGATTTAAAAAGTATCACGGTAATCGCTGCTAGTGAAAATCTTAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCt                                           <  2:2317775/100‑1 (MQ=255)
                                                      ccGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCtt                                          >  1:5400458/1‑100 (MQ=255)
                                                      ccGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCtt                                          <  1:5125215/100‑1 (MQ=255)
                                                      ccGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCtt                                          >  1:1218862/1‑100 (MQ=255)
                                                      ccGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCtt                                          <  1:4697372/100‑1 (MQ=255)
                                                       cGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTg                                         <  2:5827/100‑1 (MQ=255)
                                                       cGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTg                                         <  1:461611/100‑1 (MQ=255)
                                                       cGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTg                                         <  2:4004567/100‑1 (MQ=255)
                                                       cGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCAACGCACAGGCCAACGTAAGCAAATAACtta                                         >  2:3269078/1‑99 (MQ=255)
                                                        gggAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGc                                        <  2:5293647/100‑1 (MQ=255)
                                                        gggAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGc                                        <  1:2006030/100‑1 (MQ=255)
                                                        gggAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGc                                        <  1:3918000/100‑1 (MQ=255)
                                                            gTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATATTTGCTGac                                    <  1:1552292/100‑2 (MQ=255)
                                                            gTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGAt                                    >  2:1552596/1‑100 (MQ=255)
                                                            gTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGAt                                    >  1:1655934/1‑100 (MQ=255)
                                                             ttAAAAATTATCTCGGTAATCGCCGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGAtt                                   <  1:4685973/100‑1 (MQ=255)
                                                              tAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAGATATTTGCTGacgt                                  <  1:2036709/100‑4 (MQ=255)
                                                              tAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGAttt                                  >  1:4791956/1‑100 (MQ=255)
                                                              tAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGAttt                                  >  2:1436862/1‑100 (MQ=255)
                                                               aaaaaTTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGAtttt                                 <  2:785149/100‑1 (MQ=255)
                                                                aaaaTTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTg                                <  2:2820537/100‑1 (MQ=255)
                                                                aaaaTTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTg                                >  1:3752881/1‑100 (MQ=255)
                                                                aaaaTTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTCCTGATTTTg                                <  1:5148464/100‑1 (MQ=255)
                                                                 aaaTTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGa                               >  2:4697541/1‑100 (MQ=255)
                                                                   aTTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAAt                             <  1:3196200/100‑1 (MQ=255)
                                                                   aTTATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAAt                             <  2:4505649/100‑1 (MQ=255)
                                                                    ttATCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATc                            <  1:1404336/100‑1 (MQ=255)
                                                                    ttATCTCGGTAATCGCTGCCAGTAAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATc                            >  2:4599886/1‑100 (MQ=255)
                                                                      aTCTCGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCaa                          >  1:4468527/1‑100 (MQ=255)
                                                                       tctcGGTAATCGCTGCCAGTGAAAATCTGAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCaaa                         <  1:5055311/100‑1 (MQ=255)
                                                                       tctcGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCaaa                         <  2:662966/100‑1 (MQ=255)
                                                                       tctcGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCaaa                         <  2:1144463/100‑1 (MQ=255)
                                                                        ctcGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCaaaa                        <  1:4181505/100‑1 (MQ=255)
                                                                         tcGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCaaacg                       >  2:5077887/1‑98 (MQ=255)
                                                                         tcGGTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAg                       >  2:1983367/1‑100 (MQ=255)
                                                                           ggTAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTg                     <  2:2489324/100‑1 (MQ=255)
                                                                             tAATCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTggg                   <  2:508802/100‑1 (MQ=255)
                                                                              aaTCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGt                  <  1:3532057/100‑1 (MQ=255)
                                                                               aTCGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGtt                 <  2:1003713/100‑1 (MQ=255)
                                                                                 cGCTGCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTa               >  1:4921011/1‑100 (MQ=255)
                                                                                  gCTGCAAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTaa              >  2:3891746/1‑100 (MQ=255)
                                                                                     gCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAg           <  1:1014835/100‑1 (MQ=255)
                                                                                     gCCAGTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAg           <  1:1678438/100‑1 (MQ=255)
                                                                                         gTGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCtt       >  2:2985287/1‑100 (MQ=255)
                                                                                          tGAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCttt      <  1:685065/100‑1 (MQ=255)
                                                                                           gAAAATCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTg     <  2:3654895/100‑1 (MQ=255)
                                                                                              aaTCTCAGCCATGCCGCGACTGTACTTGGCATCGCATAGGCCAACGTAAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGAcc  <  2:3871574/100‑1 (MQ=255)
                                                                                              aaTCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCTAATATCTTGCTGATTTTGAAGCAAAAGTGGGTTTAAAAGTCTTTGAcc  <  2:5473579/100‑1 (MQ=255)
                                                                                              aaTCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGAcc  <  2:2794578/100‑1 (MQ=255)
                                                                                              aaTCTCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAAAGTCAGAAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGAcc  <  2:4163778/100‑1 (MQ=255)
                                                                                                 |                                                                                                 
ATATAAAATTATATATAATTGGGCTGTTGCGATGACAAAATTACAACTTAAGTACCGGGAGTTAAAAATTATCTCGGTAATCGCTGCCAGTGAAAATATCAGCCATGCCGCGACTGTACTTGGCATCGCACAGGCCAACGTCAGCAAATATCTTGCTGATTTTGAATCAAAAGTGGGTTTAAAAGTCTTTGACC  >  NC_000913/3751961‑3752154

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: