Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,478,879 T→C 67.0% V136A (GTT→GCT)  rraB → protein inhibitor of RNase E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,478,8790TC67.0% 137.0 / 96.8 115V136A (GTT→GCT) rraBprotein inhibitor of RNase E
Reads supporting (aligned to +/- strand):  ref base T (16/22);  new base C (29/48);  total (45/70)
Fisher's exact test for biased strand distribution p-value = 6.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 116 total aligned reads displayed.

TTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAAC  >  NC_000913/4478783‑4478978
                                                                                                |                                                                                                   
ttGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCg                                                                                                  >  2:175528/1‑100 (MQ=255)
ttGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCg                                                                                                  <  2:2015574/100‑1 (MQ=255)
ttGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCg                                                                                                  >  1:279939/1‑100 (MQ=255)
  gACGTTGAGTACGACGGATGGGGCACTTACTTTGTAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGTCGGCGTTCGcc                                                                                                <  1:1781109/100‑1 (MQ=255)
  gACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGcc                                                                                                <  2:3165003/100‑1 (MQ=255)
   aCGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCa                                                                                               <  1:2005975/100‑1 (MQ=255)
   aCGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCa                                                                                               >  2:2908336/1‑100 (MQ=255)
    cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCAc                                                                                              <  2:1060118/100‑1 (MQ=255)
    cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCAc                                                                                              <  1:2364070/100‑1 (MQ=255)
    cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCAc                                                                                              <  2:1230991/100‑1 (MQ=255)
    cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCAc                                                                                              <  2:3307962/100‑1 (MQ=255)
    cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCAc                                                                                              >  1:902331/1‑100 (MQ=255)
     gTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGAGGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCgagatc                                                                                             >  1:2161069/1‑95 (MQ=255)
     gTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACt                                                                                             <  1:1570129/100‑1 (MQ=255)
       tGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTaa                                                                                           >  2:1873824/1‑100 (MQ=255)
        gAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAAt                                                                                          <  2:3261988/100‑1 (MQ=255)
         aGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAAtt                                                                                         >  1:981404/1‑100 (MQ=255)
          gTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTa                                                                                        >  1:2708285/1‑100 (MQ=255)
          gTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTa                                                                                        <  2:1009073/100‑1 (MQ=255)
          cgatctcgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTa                                                                                        <  1:2391611/94‑1 (MQ=255)
             cgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGc                                                                                     <  1:2818139/100‑1 (MQ=255)
              gacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCg                                                                                    <  1:1401247/100‑1 (MQ=255)
                cgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAgatcg                                                                                  >  2:2391628/1‑96 (MQ=255)
                cgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTa                                                                                  <  2:1867727/100‑1 (MQ=255)
                  gATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAgc                                                                                >  2:1773084/1‑100 (MQ=255)
                  gATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAgc                                                                                <  2:2139930/100‑1 (MQ=255)
                      gggCACTTACGTTGAAGATCCTAACGGCGAAGAAGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAgcg                                                                               <  1:733453/97‑1 (MQ=255)
                       ggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATc                                                                           <  2:2828236/100‑1 (MQ=255)
                        gCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCGCGCCACTAATTAAGCGTAGCGCATCa                                                                          <  2:317880/100‑1 (MQ=255)
                         cACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCa                                                                          >  1:2924259/1‑99 (MQ=255)
                         cACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAg                                                                         <  2:2291763/100‑1 (MQ=255)
                            ttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGc                                                                       >  1:3263823/1‑99 (MQ=255)
                               cttTGACGATCCTCACGGCGAAGATCACCACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACAAATTAAGCGTAGCGCATCACGCCCa                                                                    >  1:944265/1‑99 (MQ=255)
                               cttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAg                                                                   >  2:2777108/1‑100 (MQ=255)
                               cttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAg                                                                   >  1:1402525/1‑100 (MQ=255)
                                 ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGt                                                                 >  2:1368536/1‑100 (MQ=255)
                                 ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGt                                                                 <  1:532289/100‑1 (MQ=255)
                                 tagAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCGCGCCACTAATTATGCGTAGCGCATCAGGCCCAGGt                                                                 <  2:921103/98‑1 (MQ=255)
                                  tGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGt                                                                 <  1:2834548/99‑1 (MQ=255)
                                  tGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTa                                                                >  1:1056461/1‑100 (MQ=255)
                                  tGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTa                                                                <  2:1170877/100‑1 (MQ=255)
                                  tGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTa                                                                <  2:965379/100‑1 (MQ=255)
                                  tGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTa                                                                <  1:2495617/100‑1 (MQ=255)
                                    aaGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGt                                                              <  1:1657475/100‑1 (MQ=255)
                                       aTCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTcn                                                           <  1:2544974/100‑2 (MQ=255)
                                       aTCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTcc                                                           <  2:84209/100‑1 (MQ=255)
                                       aTCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTcc                                                           <  1:2777085/100‑1 (MQ=255)
                                         ccTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTg                                                         <  2:2419180/100‑1 (MQ=255)
                                           tAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAGGCGTAGCGCATCCGGCCCAGGTAGTTCCTg                                                         <  1:1513087/98‑1 (MQ=255)
                                             aCGGCGAAGATGGTGACGATGAAGTTTTTGTCGATCAAGACGATGACGGAGCTCGCCAATGATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAAc                                                     <  2:434173/100‑1 (MQ=255)
                                             aCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAAc                                                     <  2:28315/100‑1 (MQ=255)
                                             aCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAAc                                                     <  1:2440627/100‑1 (MQ=255)
                                                 cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTgg                                                 <  1:665840/100‑1 (MQ=255)
                                                 cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTgg                                                 <  1:3139075/100‑1 (MQ=255)
                                                 cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTgg                                                 >  1:1776292/1‑100 (MQ=255)
                                                  gAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGt                                                <  1:2325216/100‑1 (MQ=255)
                                                   aaGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTg                                               <  1:2258922/100‑1 (MQ=255)
                                                    aGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGa                                              <  2:1824434/100‑1 (MQ=255)
                                                        ggCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAAttt                                          >  2:2003897/1‑100 (MQ=255)
                                                        ggCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAAttt                                          >  2:2516245/1‑100 (MQ=255)
                                                         gCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTg                                         <  1:222874/100‑1 (MQ=255)
                                                          cgacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGt                                        >  2:1146152/1‑100 (MQ=255)
                                                          cgacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGt                                        <  2:2322925/100‑1 (MQ=255)
                                                           gacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTa                                       <  1:1290173/100‑1 (MQ=255)
                                                           gacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTa                                       <  2:691668/100‑1 (MQ=255)
                                                             cgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGc                                     <  2:2916595/100‑1 (MQ=255)
                                                               aTGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGccc                                   >  1:2272518/1‑100 (MQ=255)
                                                               aTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGccc                                   >  2:1535419/1‑100 (MQ=255)
                                                               aTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGccc                                   >  2:758077/1‑100 (MQ=255)
                                                                tGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCg                                  <  2:1700655/100‑1 (MQ=255)
                                                                tGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGcccn                                  <  2:3362803/100‑2 (MQ=255)
                                                                tGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCg                                  <  1:3145645/100‑1 (MQ=255)
                                                                tGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCg                                  <  2:1801977/100‑1 (MQ=255)
                                                                tGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCg                                  <  1:3392248/100‑1 (MQ=255)
                                                                tGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCg                                  >  2:490016/1‑100 (MQ=255)
                                                                tGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCg                                  >  2:618909/1‑100 (MQ=255)
                                                                  aaGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGa                                <  2:1427006/100‑1 (MQ=255)
                                                                   aGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGAt                               >  2:2908900/1‑100 (MQ=255)
                                                                     aTTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATaa                             >  2:1522009/1‑100 (MQ=255)
                                                                      ttttGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAg                            <  2:279940/100‑1 (MQ=255)
                                                                      ttttGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAg                            >  2:2440188/1‑100 (MQ=255)
                                                                      ttttGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAg                            >  2:1019456/1‑100 (MQ=255)
                                                                          gTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCAt                        <  2:43245/100‑1 (MQ=255)
                                                                              aTGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCAt                    >  2:551930/1‑100 (MQ=255)
                                                                               tGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATg                   <  2:1145130/100‑1 (MQ=255)
                                                                                 aaGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGTTGAATTTGTAGCCCGGATAAGGCATTCATGTc                 <  1:1319213/100‑1 (MQ=255)
                                                                                  aGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTca                <  2:126766/100‑1 (MQ=255)
                                                                                    aCGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTcaca              >  1:1657654/1‑100 (MQ=255)
                                                                                     cGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACAt             >  2:1600739/1‑100 (MQ=255)
                                                                                     cGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACAt             >  1:139030/1‑100 (MQ=255)
                                                                                         gACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCgg         <  1:1239059/100‑1 (MQ=255)
                                                                                          aCGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGc        >  2:3360847/1‑100 (MQ=255)
                                                                                           cGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGc        >  2:3233073/1‑99 (MQ=255)
                                                                                           cGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCa       <  1:2973502/100‑1 (MQ=255)
                                                                                           cGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCa       <  2:1407770/100‑1 (MQ=255)
                                                                                            ggAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCAt      >  1:1995357/1‑100 (MQ=255)
                                                                                             gAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATa     >  2:3175589/1‑100 (MQ=255)
                                                                                              aGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATaa    <  2:2079872/100‑1 (MQ=255)
                                                                                               gTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATaaa   <  1:384894/100‑1 (MQ=255)
                                                                                                ttCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAAc  <  1:83673/100‑1 (MQ=255)
                                                                                                |                                                                                                   
TTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATAAAC  >  NC_000913/4478783‑4478978

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: