Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,478,879 T→C 100% 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% 321.8 / ‑1.3 100V136A (GTT→GCT) rraBprotein inhibitor of RNase E
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (42/55);  total (43/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.61e-01

TTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATA  >  NC_000913/4478783‑4478975
                                                                                                |                                                                                                
ttGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCg                                                                                               <  2:2144840/100‑1 (MQ=255)
ttGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCg                                                                                               <  1:978875/100‑1 (MQ=255)
    cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCAc                                                                                           <  2:271898/100‑1 (MQ=255)
    cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCAc                                                                                           >  2:1290549/1‑100 (MQ=255)
    cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCAc                                                                                           <  2:624348/100‑1 (MQ=255)
    cGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCAc                                                                                           <  2:679595/100‑1 (MQ=255)
     gTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACt                                                                                          <  1:2194625/100‑1 (MQ=255)
         aGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAAACGATGACGGAGCTCGCCACTAAtt                                                                                      >  1:594259/1‑100 (MQ=255)
           tACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTa                                                                                     <  1:2427150/99‑1 (MQ=255)
            acgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAg                                                                                   <  1:201921/100‑1 (MQ=255)
            acgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAg                                                                                   >  1:229386/1‑100 (MQ=255)
             cgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGc                                                                                  >  2:2887465/1‑100 (MQ=255)
             cgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGc                                                                                  >  2:1472805/1‑100 (MQ=255)
             cgacgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGc                                                                                  >  2:2757027/1‑100 (MQ=255)
              gacgGATCGGGCACTTACTTTGAAGATCCTATCGGCGAAGAATGCGACGATTAAGATTTTGTTAATGAATATGATGATTGAGCTCGACACTAAGTAAGc                                                                                  <  2:2687433/99‑1 (MQ=255)
               acgGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGt                                                                                <  2:2625581/100‑1 (MQ=255)
                 ggATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAg                                                                              <  1:1431443/100‑1 (MQ=255)
                  gATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAgc                                                                             >  2:522631/1‑100 (MQ=255)
                  gATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAgc                                                                             <  2:316818/100‑1 (MQ=255)
                   aTGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAgcg                                                                            <  2:2144659/100‑1 (MQ=255)
                     ggggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCa                                                                          <  2:1716846/100‑1 (MQ=255)
                     ggggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCa                                                                          >  1:2432878/1‑100 (MQ=255)
                     ggggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCa                                                                          >  2:1546324/1‑100 (MQ=255)
                      gggCACTTACTTTGAAGTTCCTAACGGCGAAGATGGCGATGATGAAGATTTTGTCGATGAAGACGATGACGTAGCTCGCCACTAATTAAGCGTAGCGCa                                                                          <  2:1177112/99‑1 (MQ=255)
                       ggCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATc                                                                        >  2:1967824/1‑100 (MQ=255)
                        gCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATTAAGATTTTGTCGATGAAGACGATGACGGAGCTCACCACTAATTAAGCGTAGCGCATCa                                                                       >  2:1647686/1‑100 (MQ=255)
                         cACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAg                                                                      >  1:2764514/1‑100 (MQ=255)
                           cttacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGc                                                                    >  2:664800/1‑100 (MQ=255)
                             tacttTGAAGATCCTAACGGCGAAGATTGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAgggca                                                                  >  2:1319372/1‑97 (MQ=255)
                             tacttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGccc                                                                  >  2:2090053/1‑100 (MQ=255)
                              acttTGAAGATCCTAACGGCGAAGATGGCGACGCTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCGTCGGGCCCa                                                                 <  2:685223/100‑1 (MQ=255)
                              acttTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCa                                                                 >  2:202046/1‑100 (MQ=255)
                                tttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAgg                                                               >  2:330777/1‑100 (MQ=255)
                                tttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAgg                                                               <  1:664804/100‑1 (MQ=255)
                                 ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGt                                                              <  2:2943244/100‑1 (MQ=255)
                                 ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGt                                                              >  2:2386176/1‑100 (MQ=255)
                                 ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGt                                                              <  2:796244/100‑1 (MQ=255)
                                 ttGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAAGt                                                              <  2:1248029/100‑1 (MQ=255)
                                  tGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCTACTAATTAAGCGTAGCGCATCAGGCCCAGGTa                                                             <  1:2888838/100‑1 (MQ=255)
                                  tGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTa                                                             >  1:425990/1‑100 (MQ=255)
                                   gAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAg                                                            >  1:1020276/1‑100 (MQ=255)
                                   gAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAg                                                            >  1:814814/1‑100 (MQ=255)
                                      gATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAACTCGCCACTAATTAAGCGTAACGCATCAGGCCCAGGTAGTTc                                                         >  2:2663161/1‑100 (MQ=255)
                                       aTCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTcc                                                        <  1:3017971/100‑1 (MQ=255)
                                       aTCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTcc                                                        <  1:747444/100‑1 (MQ=255)
                                        tCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCt                                                       >  1:908558/1‑100 (MQ=255)
                                        tCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCt                                                       >  2:1415688/1‑100 (MQ=255)
                                          cTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGc                                                     <  2:814810/100‑1 (MQ=255)
                                            aaCGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCaa                                                   >  2:1228590/1‑100 (MQ=255)
                                            aaCGGCGAAGATGGAGACGATGAAGATTTTGTCGATGAAGACGATGAAGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCAAGGTAGTTCCttc                                                     >  2:880957/1‑96 (MQ=255)
                                                gcgacgatggcgacgatgAAGATGTTGTCGAAGAAGAAGATGACGGAGCTGGCTACTACTTAAGCGTTACTCATCAGGCCGAGGTAACTCCTGCa                                                    >  2:1843276/1‑95 (MQ=255)
                                                 cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTgg                                              <  2:2058862/100‑1 (MQ=255)
                                                 cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTgg                                              <  1:2586976/100‑1 (MQ=255)
                                                 cGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGATAGTTCCTGCAACCTgg                                              >  2:850584/1‑100 (MQ=255)
                                                  gAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGt                                             <  1:2678570/100‑1 (MQ=255)
                                                   aaGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTg                                            >  2:1685684/1‑100 (MQ=255)
                                                   aaGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTg                                            <  2:272959/100‑1 (MQ=255)
                                                    aGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGa                                           <  1:1962581/100‑1 (MQ=255)
                                                       tGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCTACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAAtt                                        <  2:284848/100‑1 (MQ=255)
                                                       tGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAAtt                                        <  2:2531437/100‑1 (MQ=255)
                                                         gcgcagaTGAAGATTTTGTCGATGAAGACTATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTg                                      <  2:1329416/95‑1 (MQ=255)
                                                         gCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTg                                      >  1:1394584/1‑100 (MQ=255)
                                                          cgacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTCATTAAGCGTAGCGCATCAAGACCAGGAAGTTACTGCAAACTGGTGAATTTGt                                     >  1:343901/1‑100 (MQ=255)
                                                          cgacgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGt                                     <  2:1214090/100‑1 (MQ=255)
                                                            acgaTGAAGATTTTGTCGATGGAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAg                                   <  1:2503396/100‑1 (MQ=255)
                                                            acgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTa                                    <  2:172456/99‑1 (MQ=255)
                                                            acgaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAg                                   <  1:1233962/100‑1 (MQ=255)
                                                              gaTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGcc                                 >  2:303442/1‑100 (MQ=255)
                                                               aTGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGccc                                <  2:2497767/100‑1 (MQ=255)
                                                                tGAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCg                               <  1:522634/100‑1 (MQ=255)
                                                                 gAAGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAACACATCAGGCCCAGGTAGTTCCTGAAACCTGGTGAATTTGTAGCCCgg                              >  1:2500663/1‑100 (MQ=255)
                                                                  aaGATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGa                             <  1:1239475/100‑1 (MQ=255)
                                                                    gATTTTGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATa                           <  1:554352/100‑1 (MQ=255)
                                                                      ttttGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAg                         <  2:2471990/100‑1 (MQ=255)
                                                                      ttttGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAg                         <  2:907305/100‑1 (MQ=255)
                                                                       tttGTGGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGTATCAGGTCCAGGTAGTTCCTGCAACCTGGTGATTTTGTAGCCCGGATAAgg                        <  2:433960/100‑1 (MQ=255)
                                                                       tttGTCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAgg                        <  2:2764482/100‑1 (MQ=255)
                                                                           tCGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGATCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCAtt                    >  2:2129089/1‑100 (MQ=255)
                                                                            cGATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTc                   <  2:1394571/100‑1 (MQ=255)
                                                                             gATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCa                  <  2:2297269/100‑1 (MQ=255)
                                                                             gATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCa                  >  1:1554657/1‑100 (MQ=255)
                                                                             gATGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCa                  >  1:2310805/1‑100 (MQ=255)
                                                                             gATGAAGACGATGACGGAGCTCGCCACTAATCAAGCGGAGCGCATCATGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCa                  <  2:2050775/100‑1 (MQ=255)
                                                                              aTGAAGACGACGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCAt                 >  1:1955914/1‑100 (MQ=255)
                                                                               tGAAGACGATGACGGAGCTCGCCACTAATTAAGCGTAGTGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATg                <  1:952594/100‑1 (MQ=255)
                                                                                 aaGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTc              <  2:2432849/100‑1 (MQ=255)
                                                                                  aGACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTca             >  1:1411367/1‑100 (MQ=255)
                                                                                   gACGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTcac            <  1:1723406/100‑1 (MQ=255)
                                                                                     cGATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTAGTAGACTGGATAAGACATTCATGTCAAAt          >  1:2078524/1‑100 (MQ=255)
                                                                                      gATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTcacaac         >  2:1832217/1‑98 (MQ=255)
                                                                                      gATGACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTcaca           >  2:1832315/1‑98 (MQ=255)
                                                                                         gACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCgg      <  1:1290561/100‑1 (MQ=255)
                                                                                         gACGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCgg      <  1:2710966/100‑1 (MQ=255)
                                                                                          aCGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGc     >  2:33211/1‑100 (MQ=255)
                                                                                          aCGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGc     <  1:2329282/100‑1 (MQ=255)
                                                                                          aCGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGc     <  1:330780/100‑1 (MQ=255)
                                                                                           cGGAGCTCGCCAGTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCa    >  1:1997196/1‑100 (MQ=255)
                                                                                           cGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCa    <  2:2559298/100‑1 (MQ=255)
                                                                                           cGGAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCa    <  2:961379/100‑1 (MQ=255)
                                                                                             gAGCTCGCCACTAATTAAGCGTAGCGCATCAGGCCTAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATa  <  1:816990/100‑1 (MQ=255)
                                                                                                |                                                                                                
TTGACGTTGAGTACGACGGATGGGGCACTTACTTTGAAGATCCTAACGGCGAAGATGGCGACGATGAAGATTTTGTCGATGAAGACGATGACGGAGTTCGCCACTAATTAAGCGTAGCGCATCAGGCCCAGGTAGTTCCTGCAACCTGGTGAATTTGTAGCCCGGATAAGGCATTCATGTCACATCCGGCATA  >  NC_000913/4478783‑4478975

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