Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,478,879 T→C 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,8790TC97.0% 326.2 / ‑0.7 99V136A (GTT→GCT) rraBprotein inhibitor of RNase E
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (40/56);  total (41/58)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

TTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCA  >  NC_000913/4478783‑4478973
                                                                                                |                                                                                              
ttGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCg                                                                                             <  2:1441369/100‑1 (MQ=255)
 tGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGc                                                                                            <  2:2760224/100‑1 (MQ=255)
  gACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGcc                                                                                           <  2:1223763/100‑1 (MQ=255)
  gACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGcc                                                                                           <  2:2429898/100‑1 (MQ=255)
    cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCAc                                                                                         <  1:41803/100‑1 (MQ=255)
     gTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACt                                                                                        >  2:666374/1‑100 (MQ=255)
     gTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACt                                                                                        >  1:1140659/1‑100 (MQ=255)
     gTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACt                                                                                        <  2:2823556/100‑1 (MQ=255)
         aGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAAtt                                                                                    >  2:949638/1‑100 (MQ=255)
           tACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTaa                                                                                  <  2:959664/100‑1 (MQ=255)
            acgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAg                                                                                 <  1:2962747/100‑1 (MQ=255)
            acgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAg                                                                                 <  2:2943895/100‑1 (MQ=255)
            acgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAg                                                                                 >  1:2591689/1‑100 (MQ=255)
             cgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGc                                                                                >  2:140581/1‑100 (MQ=255)
             cgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGc                                                                                >  2:648174/1‑100 (MQ=255)
             cgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGc                                                                                <  1:2610676/100‑1 (MQ=255)
              gacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCg                                                                               <  2:2293610/100‑1 (MQ=255)
              gacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCg                                                                               <  2:202863/100‑1 (MQ=255)
              gacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCg                                                                               <  2:625206/100‑1 (MQ=255)
              gacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCg                                                                               >  1:1017147/1‑100 (MQ=255)
               acgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGCTGGCGACGATGAAGATCTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCg                                                                               >  1:2159024/1‑99 (MQ=255)
                 ggATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAg                                                                            <  2:2961404/100‑1 (MQ=255)
                  gATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAgc                                                                           <  1:2721784/100‑1 (MQ=255)
                   aTGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAgcg                                                                          <  2:2887910/100‑1 (MQ=255)
                   aTGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAgcg                                                                          >  2:3040281/1‑100 (MQ=255)
                      gggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCAt                                                                       <  1:2680577/100‑1 (MQ=255)
                       ggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATc                                                                      <  1:2477839/100‑1 (MQ=255)
                       ggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATc                                                                      <  1:2563071/100‑1 (MQ=255)
                        gCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCa                                                                     <  2:845398/100‑1 (MQ=255)
                        gCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCa                                                                     <  2:14981/100‑1 (MQ=255)
                        gCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCa                                                                     <  1:2806489/100‑1 (MQ=255)
                         cACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAg                                                                    >  1:2256551/1‑100 (MQ=255)
                         cACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAg                                                                    <  2:584572/100‑1 (MQ=255)
                           cttacttTGAAGATCCTAACGGCGAAGATGGCCACGCTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGc                                                                  >  1:2909447/1‑100 (MQ=255)
                              acttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCa                                                               <  2:48232/100‑1 (MQ=255)
                              acttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCa                                                               >  1:1870261/1‑100 (MQ=255)
                               cttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAg                                                              <  1:180336/100‑1 (MQ=255)
                                  tGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTa                                                           <  1:223552/100‑1 (MQ=255)
                                  tGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTa                                                           <  2:2706431/100‑1 (MQ=255)
                                   gAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAg                                                          >  2:3023891/1‑100 (MQ=255)
                                   gAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAg                                                          <  1:523616/100‑1 (MQ=255)
                                      gATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTc                                                       >  1:2525959/1‑100 (MQ=255)
                                       aTCCTAACGGCGAAGATGGCGACGATGACGATTTTGTCGATGAAGACGATGACGCAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTcc                                                      <  1:1025414/100‑1 (MQ=255)
                                         ccTAAGGTCGAAGATGGCGACGATGAAGATTTAGTCGATGAAGACGATGACGGAGCACGCCACTAATTAAGCGTAGCGAATCAGGCCCAGGTAGTTCCt                                                     <  1:2828827/99‑1 (MQ=255)
                                          cTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGACACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGc                                                   >  2:2900686/1‑100 (MQ=255)
                                           tAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGcn                                                  <  1:3008598/100‑2 (MQ=255)
                                           tAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCa                                                  <  2:1047480/100‑1 (MQ=255)
                                            aaCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCaa                                                 <  1:1516096/100‑1 (MQ=255)
                                             aCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAAc                                                <  1:3071660/100‑1 (MQ=255)
                                               ggCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCt                                              >  1:2272006/1‑100 (MQ=255)
                                                gCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTg                                             >  2:401275/1‑100 (MQ=255)
                                                 cgatcttGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTgg                                            <  2:3039234/94‑1 (MQ=255)
                                                 ccaaGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTgg                                            <  1:2103762/98‑1 (MQ=255)
                                                 cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGTTGTCGGAGCTCGCCACTAATTAAGCGTAACGCATCAGGCCCAGGTAGTTCCTGCAActt                                              >  1:2447558/1‑96 (MQ=255)
                                                 cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTgg                                            >  1:2842463/1‑100 (MQ=255)
                                                 cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTgg                                            <  2:306847/100‑1 (MQ=255)
                                                   aaGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTg                                          <  1:2914317/100‑1 (MQ=255)
                                                    aGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGa                                         <  2:316050/100‑1 (MQ=255)
                                                    aGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGa                                         <  1:2891607/100‑1 (MQ=255)
                                                       tGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGgagatcg                                      >  1:3039185/1‑97 (MQ=255)
                                                       tGGAGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAAtt                                      >  1:2918653/1‑100 (MQ=255)
                                                        ggCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAAttt                                     >  1:1332124/1‑100 (MQ=255)
                                                         gCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTg                                    >  1:214124/1‑100 (MQ=255)
                                                         gCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTg                                    >  1:2170372/1‑100 (MQ=255)
                                                          cgacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGt                                   >  2:2538183/1‑100 (MQ=255)
                                                            acgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTGGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAg                                 <  1:2622276/100‑1 (MQ=255)
                                                            acgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAg                                 >  1:1955980/1‑100 (MQ=255)
                                                            acgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAg                                 <  1:140579/100‑1 (MQ=255)
                                                            acgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAg                                 <  1:2622448/100‑1 (MQ=255)
                                                             cgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGc                                <  1:934457/100‑1 (MQ=255)
                                                               aTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGccc                              >  1:3281865/1‑100 (MQ=255)
                                                               aTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGccc                              <  1:515479/100‑1 (MQ=255)
                                                                tGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCg                             <  1:407648/100‑1 (MQ=255)
                                                                 gAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCgg                            >  1:1922362/1‑100 (MQ=255)
                                                                 gAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCgg                            <  2:649946/100‑1 (MQ=255)
                                                                  aaGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGa                           >  2:114860/1‑100 (MQ=255)
                                                                   aGATTTTGTCGGTGAAGACGATGACGGAGCACGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTCGTTGCTGCAACCTGGTGCATTTGTAGCCCgg                            <  1:1059267/98‑1 (MQ=255)
                                                                   aGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGAt                          <  1:139241/100‑1 (MQ=255)
                                                                     aTTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATaa                        >  2:447524/1‑100 (MQ=255)
                                                                     aTTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTACTGAAACCTGGTGAATTTGTAGCCCGGATaa                        >  2:178784/1‑100 (MQ=255)
                                                                      ttttGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAg                       >  2:2509471/1‑100 (MQ=255)
                                                                      ttttGTCGATGAAGACGATGACGGAGCGCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAg                       <  1:2870281/100‑1 (MQ=255)
                                                                       tttGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAgg                      <  2:1265533/100‑1 (MQ=255)
                                                                       tttGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAgg                      <  1:1584137/100‑1 (MQ=255)
                                                                        ttGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGc                     <  1:2921682/100‑1 (MQ=255)
                                                                          gTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCAt                   >  1:1471211/1‑100 (MQ=255)
                                                                           tCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCAtt                  >  2:2172713/1‑100 (MQ=255)
                                                                             gATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTc                 >  2:1885104/1‑99 (MQ=255)
                                                                             gATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCa                >  2:1949076/1‑100 (MQ=255)
                                                                             gATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCa                >  1:3332142/1‑100 (MQ=255)
                                                                               tGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATg              <  2:314469/100‑1 (MQ=255)
                                                                               tGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATg              <  1:583297/100‑1 (MQ=255)
                                                                                 aaGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTc            <  2:3036408/100‑1 (MQ=255)
                                                                                  aGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTca           >  2:3020403/1‑100 (MQ=255)
                                                                                  aGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTca           >  2:609943/1‑100 (MQ=255)
                                                                                   gactaggACGGGGCTCGCCACTAATTAAGCGTAGTGCATCAGGCCCACGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTcac          <  1:1011906/94‑1 (MQ=255)
                                                                                   gACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTcac          >  2:1022798/1‑100 (MQ=255)
                                                                                     cGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACAt        >  1:1838143/1‑100 (MQ=255)
                                                                                           cGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCa  <  1:1116091/100‑1 (MQ=255)
                                                                                                |                                                                                              
TTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCA  >  NC_000913/4478783‑4478973

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