Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,966,721 G→A R102H (CGC→CAC)  rho → transcription termination factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,966,7210GA99.1% 386.8 / ‑4.7 114R102H (CGC→CAC) rhotranscription termination factor
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (63/50);  total (63/51)
Fisher's exact test for biased strand distribution p-value = 4.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

Only 100 of 121 total aligned reads displayed.

AGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGA  >  NC_000913/3966624‑3966817
                                                                                                 |                                                                                                
aGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcacc                                                                                                <  2:2033618/100‑1 (MQ=255)
 gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcaccc                                                                                               >  2:141606/1‑100 (MQ=255)
 gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcaccc                                                                                               >  2:2179870/1‑100 (MQ=255)
  cTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcacccg                                                                                              >  2:2845438/1‑100 (MQ=255)
   tCCTACCTCGCCGTTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcacccg                                                                                              <  1:351083/99‑1 (MQ=255)
   tCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcacccgc                                                                                             >  2:1549604/1‑100 (MQ=255)
     cTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTcacccgcc                                                                                            <  1:2286550/99‑1 (MQ=255)
      tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGa                                                                                          >  2:2321975/1‑100 (MQ=255)
      tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGa                                                                                          >  1:2076385/1‑100 (MQ=255)
       aCCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGaa                                                                                         >  2:251948/1‑100 (MQ=255)
        ccTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGaaa                                                                                        >  1:885870/1‑100 (MQ=255)
          tCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAaaga                                                                                      <  1:1385558/100‑1 (MQ=255)
           cgctgtaccTGATGACATCTACGTTTCCCCCTACCAAATCCCGCGTATCAACCTCCACACAGGTCATACCATCTCTGTTATGATTCACCCACCGAaag                                                                                       <  1:2723344/91‑1 (MQ=255)
             ccGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAgg                                                                                   >  2:737561/1‑100 (MQ=255)
              cGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGt                                                                                  <  2:1900881/100‑1 (MQ=255)
                gTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCCGGTAAGATTCACCCGCCGAAAGAAGGTGa                                                                                >  1:2763769/1‑100 (MQ=255)
                 tCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGa                                                                                >  2:2939018/1‑99 (MQ=255)
                  ccTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAAc                                                                              >  1:15043/1‑100 (MQ=255)
                   cTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACg                                                                             <  2:1457968/100‑1 (MQ=255)
                   cTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACg                                                                             >  1:644463/1‑100 (MQ=255)
                    tgatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGc                                                                            >  2:1635646/1‑100 (MQ=255)
                    tgatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGc                                                                            <  1:2321966/100‑1 (MQ=255)
                      atgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTa                                                                          >  2:1677752/1‑100 (MQ=255)
                         aCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTAttt                                                                       <  2:92012/100‑1 (MQ=255)
                         aCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTAtt                                                                        <  2:1029050/99‑1 (MQ=255)
                          cATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTAtttt                                                                      >  1:824018/1‑100 (MQ=255)
                           aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTg                                                                     >  1:1110694/1‑100 (MQ=255)
                           aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTg                                                                     >  1:783279/1‑100 (MQ=255)
                           aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTg                                                                     >  1:496970/1‑100 (MQ=255)
                           aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTg                                                                     <  1:278022/100‑1 (MQ=255)
                           aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTg                                                                     <  1:2733447/100‑1 (MQ=255)
                            tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                    >  2:2121256/1‑100 (MQ=255)
                            tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTAGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTAtttt                                                                      >  1:2937808/1‑98 (MQ=255)
                             cTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTgcg                                                                   >  1:819509/1‑100 (MQ=255)
                              tACGTTTCCCCTACCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                  <  2:483694/100‑1 (MQ=255)
                               aCGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCgct                                                                 >  1:1762108/1‑100 (MQ=255)
                                   ttCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGctgctg                                                             <  2:2454832/100‑1 (MQ=255)
                                     ccccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaa                                                           >  1:2425410/1‑100 (MQ=255)
                                      cccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                          >  2:2245415/1‑100 (MQ=255)
                                      cccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                          >  1:1294309/1‑100 (MQ=255)
                                      cccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                          <  2:3032762/100‑1 (MQ=255)
                                       ccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAg                                                         <  1:575686/100‑1 (MQ=255)
                                       ccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAg                                                         <  2:2762355/100‑1 (MQ=255)
                                        cTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGt                                                        <  1:1743510/100‑1 (MQ=255)
                                         tAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGtt                                                       >  2:1852445/1‑100 (MQ=255)
                                            ccAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAAc                                                    >  2:2698612/1‑100 (MQ=255)
                                             cAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                   <  2:2076393/100‑1 (MQ=255)
                                                 tCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGt                                               <  1:1162838/100‑1 (MQ=255)
                                                    gccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTaa                                            >  1:2165983/1‑100 (MQ=255)
                                                    gccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTaa                                            >  1:2569438/1‑100 (MQ=255)
                                                      cgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACt                                          >  2:1894928/1‑100 (MQ=255)
                                                       gTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt                                         >  2:2501824/1‑100 (MQ=255)
                                                        tttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTc                                        <  2:2680728/100‑1 (MQ=255)
                                                           cAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                     >  1:1914572/1‑100 (MQ=255)
                                                           cAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                     <  2:2221048/100‑1 (MQ=255)
                                                             aCCTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaa                                   <  2:1085210/100‑1 (MQ=255)
                                                              ccTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaaa                                  >  2:87502/1‑100 (MQ=255)
                                                              ccTCCGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaaa                                  >  2:2741601/1‑100 (MQ=255)
                                                                 ccGCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                               >  2:629458/1‑100 (MQ=255)
                                                                   gCACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGa                             >  1:123177/1‑100 (MQ=255)
                                                                    cACTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaa                            <  1:2669275/100‑1 (MQ=255)
                                                                     aCTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaa                           <  1:482066/100‑1 (MQ=255)
                                                                     aCTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaa                           <  2:1102334/100‑1 (MQ=255)
                                                                     aCTGGTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaa                           <  2:1033726/100‑1 (MQ=255)
                                                                        ggTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaaa                          >  1:1464722/1‑98 (MQ=255)
                                                                        ggTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACg                        >  1:2898073/1‑100 (MQ=255)
                                                                         gTGATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGc                       >  2:1726353/1‑100 (MQ=255)
                                                                           ggtACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGccc                     >  2:1555172/3‑100 (MQ=255)
                                                                           gATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGccc                     >  2:260943/1‑100 (MQ=255)
                                                                           gATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGccc                     >  1:2885020/1‑100 (MQ=255)
                                                                           gATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGccc                     >  1:695713/1‑100 (MQ=255)
                                                                           gATACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGccc                     <  2:1762115/100‑1 (MQ=255)
                                                                            aTACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                    <  1:1502207/100‑1 (MQ=255)
                                                                             tACCATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGc                   <  2:1084300/100‑1 (MQ=255)
                                                                               ccATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGca                  >  2:2547141/1‑99 (MQ=255)
                                                                                cATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGCAAGTTAACGAACTTACCTTAGACAACCCTGAAAACGCCCGcaac                >  1:2986561/1‑100 (MQ=255)
                                                                                cATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACTCCCGcaac                <  1:1015298/100‑1 (MQ=255)
                                                                                cATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaac                <  1:1555167/100‑1 (MQ=255)
                                                                                cATCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaac                <  2:187908/100‑1 (MQ=255)
                                                                                 aTCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaaca               >  2:1534282/1‑100 (MQ=255)
                                                                                 aTCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaaca               >  2:721906/1‑100 (MQ=255)
                                                                                 aTCTCTGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaaca               <  1:911108/100‑1 (MQ=255)
                                                                                     ctGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt           >  1:2519899/1‑100 (MQ=255)
                                                                                     ctGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt           <  1:2149277/100‑1 (MQ=255)
                                                                                     ctGGTAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACAAAGTTAACTTCGACAAACCTGAAAACGCCCGcaacaa              <  1:1635639/97‑1 (MQ=255)
                                                                                         tAAGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc       >  2:398764/1‑100 (MQ=255)
                                                                                          aaGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct      <  2:1179096/100‑1 (MQ=255)
                                                                                           aGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt     <  2:2569452/100‑1 (MQ=255)
                                                                                           aGATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCtt     <  2:1630905/100‑1 (MQ=255)
                                                                                            gATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt    >  2:1027411/1‑100 (MQ=255)
                                                                                            gATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt    >  2:529813/1‑100 (MQ=255)
                                                                                            gATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt    <  1:2121247/100‑1 (MQ=255)
                                                                                            gATTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt    >  2:2737160/1‑100 (MQ=255)
                                                                                             aTTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTATCTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg   <  1:1677745/100‑1 (MQ=255)
                                                                                             aTTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg   >  2:2239639/1‑100 (MQ=255)
                                                                                             aTTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg   <  1:771666/100‑1 (MQ=255)
                                                                                             aTTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg   >  2:1242716/1‑100 (MQ=255)
                                                                                             aTTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg   >  2:2950452/1‑100 (MQ=255)
                                                                                             aTTCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTg   >  2:2811218/1‑100 (MQ=255)
                                                                                              ttCACCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTga  >  2:2853166/1‑100 (MQ=255)
                                                                                                 |                                                                                                
AGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTTGA  >  NC_000913/3966624‑3966817

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