Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,366,015 Δ1 bp 16.5% coding (88/1140 nt) arnB → uridine 5'‑(beta‑1‑threo‑pentapyranosyl‑4‑ulose diphosphate) aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,366,0150A.16.5% 406.4 / 78.3 139coding (88/1140 nt)arnBuridine 5'‑(beta‑1‑threo‑pentapyranosyl‑4‑ulose diphosphate) aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base A (59/57);  new base . (11/12);  total (70/69)
Fisher's exact test for biased strand distribution p-value = 8.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.59e-01

Only 100 of 145 total aligned reads displayed.

GAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCAT  >  NC_000913/2365916‑2366113
                                                                                                    |                                                                                                   
gaaggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCAACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATca                                                                                                     <  2:136500/100‑1 (MQ=255)
 aaggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATcac                                                                                                    <  2:2371893/100‑1 (MQ=255)
  aggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATcaca                                                                                                   <  1:2173519/100‑1 (MQ=255)
    gaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAAc                                                                                                 <  2:2661057/100‑1 (MQ=255)
     aaaaGCAATGTCAGAATTTTTTTTTTTTTTTTTACAAGCAATGGGCGTGGAGGAACTTGCTGCAGTTAAAGAGGTTCTCGAATTCGGTTGGATCACAAc                                                                                                 >  1:3683919/1‑99 (MQ=255)
     aaaaGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAAcc                                                                                                <  2:1929585/100‑1 (MQ=255)
      aaaGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCg                                                                                               >  2:1368154/1‑100 (MQ=255)
        aGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGt                                                                                            <  1:456673/100‑1 (MQ=255)
         gCAATGTCAGAATTTTTTTTTTTTTTTTTTCCAGCAATGGGCGGGGAAAAACTTTCTTTTGTTAAAGAGGTTCTCGAATCTGGTTCTATCACAACCg                                                                                               >  1:1378590/1‑97 (MQ=255)
         gCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATCCAACCggat                                                                                            <  2:1301983/100‑3 (MQ=255)
         gCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGCTGTGATCCAACCggat                                                                                            <  2:1894933/100‑3 (MQ=255)
          cAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATCCAACCggatt                                                                                           <  2:1847567/100‑4 (MQ=255)
          cAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATCCAACCggatt                                                                                           <  1:1384671/100‑4 (MQ=255)
          cAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATCCAACCggatt                                                                                           <  1:248040/100‑4 (MQ=255)
          cAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTcc                                                                                          <  2:1320215/100‑1 (MQ=255)
           aaTGTCAGAATTTTTGCCTTTTTCGCGACCATCCATGGGCCTGGCGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGTTTGTGATCCAACCggattc                                                                                          <  1:501440/100‑5 (MQ=255)
           aaTGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATCCAACCggattc                                                                                          <  1:2233067/100‑5 (MQ=255)
           aaTGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATCCAACCggattc                                                                                          <  2:532858/100‑5 (MQ=255)
           aaTGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATCCAACCggattc                                                                                          <  2:988642/100‑5 (MQ=255)
             tGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCgaa                                                                                       <  2:1951329/100‑1 (MQ=255)
             tGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCgaa                                                                                       <  2:1758088/100‑1 (MQ=255)
               tCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCgaaga                                                                                     >  2:3651929/1‑100 (MQ=255)
                cAGAATTTTTTTTTTTTTTTTGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTcttaa                                                                                       >  1:1238424/1‑93 (MQ=255)
                cAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATCCAACCggattcgagaa                                                                                     <  2:369376/100‑10 (MQ=38)
                cAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCgaagaa                                                                                    >  1:1016197/1‑100 (MQ=255)
                cAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCgaagaa                                                                                    >  2:718769/1‑100 (MQ=255)
                cAGAATGTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGGTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATCCAACCggattcgagaa                                                                                     <  1:345113/100‑10 (MQ=37)
                 aGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAAt                                                                                   >  2:2814673/1‑100 (MQ=255)
                  gAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATc                                                                                  >  2:824671/1‑100 (MQ=255)
                   aaTTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCa                                                                                 <  2:1824352/100‑1 (MQ=255)
                   aaTTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCa                                                                                 <  2:2702272/100‑1 (MQ=255)
                   aaTTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCa                                                                                 <  2:343011/100‑1 (MQ=255)
                      ttttGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGgc                                                                              <  2:766207/100‑1 (MQ=255)
                          gCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCtt                                                                          <  2:216470/100‑1 (MQ=255)
                           ccTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTg                                                                         <  2:1727324/100‑1 (MQ=255)
                           ccTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTg                                                                         >  2:3250274/1‑100 (MQ=255)
                            cTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGa                                                                        >  2:2305343/1‑100 (MQ=255)
                             tttttCGCGACCAGTAATGGGCGTGGAGGAACTCGTTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGa                                                                        <  1:3182518/99‑1 (MQ=255)
                             tttttCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGa                                                                        >  2:834126/1‑99 (MQ=255)
                             tttttCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAg                                                                       >  1:427186/1‑100 (MQ=255)
                             tttttCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAg                                                                       <  2:942826/100‑1 (MQ=255)
                              ttttCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGc                                                                      <  2:427189/100‑1 (MQ=255)
                              ttttCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGc                                                                      <  2:79233/100‑1 (MQ=255)
                                 tCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAg                                                                   <  1:1970717/100‑1 (MQ=255)
                                  cgcgACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGTGCTTGAGCAAGc                                                                  >  2:3083519/1‑100 (MQ=255)
                                  cgcgACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGc                                                                  <  1:3193591/100‑1 (MQ=255)
                                  cgcgACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGc                                                                  >  1:2302380/1‑100 (MQ=255)
                                       ccATCAATGGGCGTGGAGGAACTCGCTCCAGTTAAAGAGGTTCTCGAACCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCttttt                                                             <  1:3661272/100‑1 (MQ=255)
                                       ccAGCAATGGNAGNGGNGGNACTCGCTGCANTTNAAGAGGTNCTCGNATCCGGTTGGATCNCCACCGGTNCGAAGAATCAGGCGCTTGAGCAAGCatnt                                                              >  1:3010789/1‑99 (MQ=255)
                                       ccAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCttttt                                                             <  1:1007711/100‑1 (MQ=255)
                                        cAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                            <  2:3690488/100‑1 (MQ=255)
                                         aGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGc                                                           >  1:2801528/1‑100 (MQ=255)
                                         aGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGc                                                           >  2:2068587/1‑100 (MQ=255)
                                         aGCAATGGGCGTGGAGGAACTCGCCACAGCTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAATAATCAGGCGCTTGAGCAAGCTTTTTg                                                            <  1:3646659/99‑1 (MQ=255)
                                          gCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGcc                                                          <  1:1780166/100‑1 (MQ=255)
                                          gCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGcc                                                          <  2:3566710/100‑1 (MQ=255)
                                              tGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCCGGCGCTTGAGCAAGCTTTTTGCCAGtt                                                      >  1:815433/1‑100 (MQ=255)
                                              tGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGtt                                                      >  1:2179932/1‑100 (MQ=255)
                                              tGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGtt                                                      >  2:3434302/1‑100 (MQ=255)
                                                ggCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGGGCAAGCTTTTTGCCAGTTGa                                                    >  1:2031712/1‑100 (MQ=255)
                                                   gTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACgg                                                 >  2:1662751/1‑100 (MQ=255)
                                                    ttgaggaACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACggg                                                <  1:2359756/98‑1 (MQ=255)
                                                    tGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACggg                                                <  1:574955/100‑1 (MQ=255)
                                                     ggaggaACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGa                                               <  2:513937/100‑1 (MQ=255)
                                                        ggaACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAAt                                            <  2:1668359/100‑1 (MQ=255)
                                                        ggaACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAAt                                            <  2:3612496/100‑1 (MQ=255)
                                                          aaCTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCa                                          >  2:3781639/1‑100 (MQ=255)
                                                               gCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATg                                     <  2:1449006/100‑1 (MQ=255)
                                                                cTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATg                                     >  1:1728068/1‑99 (MQ=255)
                                                                 ngCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGcc                                   >  2:3506505/2‑100 (MQ=255)
                                                                   cAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCAt                                 >  2:449520/1‑100 (MQ=255)
                                                                    aGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATc                                >  1:535639/1‑100 (MQ=255)
                                                                    aGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATc                                >  2:3193896/1‑100 (MQ=255)
                                                                     gTTAAAGAGGTTCTCGAATCCGGTTGGATAACAACCGGTCCGAAGAATCAGGCGCGTGAGCAAGATTTTTGCCATTTGACGGGAAATCAGCATGCCAAcg                               <  1:3301540/100‑1 (MQ=255)
                                                                       tAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATcgcg                             >  1:1181677/1‑100 (MQ=255)
                                                                         aaGAGGTTCTCGAATCCGGTTGGATCAAAAACGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGt                           >  1:1651660/1‑100 (MQ=255)
                                                                              attcttcggaccGGTTGTGATCCAACCGGATTCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGt                       >  1:275383/11‑100 (MQ=21)
                                                                               ttcttcggaccGGTTGTGATCCAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTc                     >  2:304505/10‑100 (MQ=38)
                                                                               ttCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTc                     <  1:2497061/100‑1 (MQ=255)
                                                                                tcttcggaccGGTTGTGATCCAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATGGCGGTCAGTTc                     >  2:2344082/9‑99 (MQ=37)
                                                                                tctcGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCa                    >  1:952317/1‑100 (MQ=255)
                                                                                 cttcggaccGGTTGTGATCCAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAg                   >  2:2570564/8‑100 (MQ=38)
                                                                                   tcggaccGGTTGTGATCCAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGcc                 >  2:2254919/6‑100 (MQ=255)
                                                                                   cGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGcc                 >  2:2828731/1‑100 (MQ=255)
                                                                                   cGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGcc                 >  1:1145143/1‑100 (MQ=255)
                                                                                   cGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGcc                 <  1:2700866/100‑1 (MQ=255)
                                                                                    cggaccGGTTGTGATCCAACCGGTCCGAAGAATCCGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCa                >  2:650592/5‑100 (MQ=255)
                                                                                    cggaccGGTTGTGATCCAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCa                >  1:2956036/5‑100 (MQ=255)
                                                                                     ggaccGGTTGTGATCCAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCAc               >  2:2154599/4‑100 (MQ=255)
                                                                                        ccGGTTGTGATCCAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCAccgc            >  1:1660131/1‑100 (MQ=255)
                                                                                        ccGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCAccgc            >  2:2051510/1‑100 (MQ=255)
                                                                                         cGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGTGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCAccgcc           >  1:2394693/1‑100 (MQ=255)
                                                                                          ggTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCAccgccg          <  1:3561107/100‑1 (MQ=255)
                                                                                           gTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCgg         <  2:2070779/100‑1 (MQ=255)
                                                                                             tGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGaa       <  1:3778468/100‑1 (MQ=255)
                                                                                             tGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGaa       <  2:2017937/100‑1 (MQ=255)
                                                                                              ggATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTTCCTGTTGATTGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAAt      >  1:1569780/1‑100 (MQ=255)
                                                                                              ggATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAAt      <  1:3372427/100‑1 (MQ=255)
                                                                                                 aTCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAAt      <  1:687149/98‑1 (MQ=255)
                                                                                                 aTCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGc    <  2:3683949/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
GAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGGATCACAACCGGTCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCAT  >  NC_000913/2365916‑2366113

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