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.

GAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCAT  >  NC_000913/2365916‑2366113
                                                                                                    |                                                                                                   
gaaggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCAACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATca                                                                                                     <  2:136500/100‑1 (MQ=255)
 aaggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATcac                                                                                                    <  2:2371893/100‑1 (MQ=255)
  aggaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATcaca                                                                                                   <  1:2173519/100‑1 (MQ=255)
    gaaAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAAc                                                                                                 <  2:2661057/100‑1 (MQ=255)
     aaaaGCAATGTCAGAATTTTTTTTTTTTTTTTTACAAGCAATGGGCGTGGAGGAACTTGCTGCAGTTAAAGAGGTTCTCGAATTCGGTTG‑GATCACAAc                                                                                                 >  1:3683919/1‑99 (MQ=255)
     aaaaGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAAcc                                                                                                <  2:1929585/100‑1 (MQ=255)
      aaaGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCg                                                                                               >  2:1368154/1‑100 (MQ=255)
        aGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑t                                                                                            <  1:456673/100‑1 (MQ=255)
         gCAATGTCAGAATTTTTTTTTTTTTTTTTTCCAGCAATGGGCGGGGAAAAACTTTCTTTTGTTAAAGAGGTTCTCGAATCTGGTTC‑TATCACAACCg                                                                                               >  1:1378590/1‑97 (MQ=255)
         gCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATC‑CAACCggat                                                                                            <  2:1301983/100‑3 (MQ=255)
         gCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGCTGTGATC‑CAACCggat                                                                                            <  2:1894933/100‑3 (MQ=255)
          cAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATC‑CAACCggatt                                                                                           <  2:1847567/100‑4 (MQ=255)
          cAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATC‑CAACCggatt                                                                                           <  1:1384671/100‑4 (MQ=255)
          cAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATC‑CAACCggatt                                                                                           <  1:248040/100‑4 (MQ=255)
          cAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑Tcc                                                                                          <  2:1320215/100‑1 (MQ=255)
           aaTGTCAGAATTTTTGCCTTTTTCGCGACCATCCATGGGCCTGGCGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGTTTGTGATC‑CAACCggattc                                                                                          <  1:501440/100‑5 (MQ=255)
           aaTGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATC‑CAACCggattc                                                                                          <  1:2233067/100‑5 (MQ=255)
           aaTGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATC‑CAACCggattc                                                                                          <  2:532858/100‑5 (MQ=255)
           aaTGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATC‑CAACCggattc                                                                                          <  2:988642/100‑5 (MQ=255)
             tGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCgaa                                                                                       <  2:1951329/100‑1 (MQ=255)
             tGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCgaa                                                                                       <  2:1758088/100‑1 (MQ=255)
               tCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCgaaga                                                                                     >  2:3651929/1‑100 (MQ=255)
                cAGAATTTTTTTTTTTTTTTTGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑Tcttaa                                                                                       >  1:1238424/1‑93 (MQ=255)
                cAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATC‑CAACCggattcgagaa                                                                                     <  2:369376/100‑10 (MQ=38)
                cAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCgaagaa                                                                                    >  1:1016197/1‑100 (MQ=255)
                cAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCgaagaa                                                                                    >  2:718769/1‑100 (MQ=255)
                cAGAATGTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGGTGCAGTTAAAGAGGTTCTCGAATCCGGTTGTGATC‑CAACCggattcgagaa                                                                                     <  1:345113/100‑10 (MQ=37)
                 aGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAAt                                                                                   >  2:2814673/1‑100 (MQ=255)
                  gAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATc                                                                                  >  2:824671/1‑100 (MQ=255)
                   aaTTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCa                                                                                 <  2:1824352/100‑1 (MQ=255)
                   aaTTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCa                                                                                 <  2:2702272/100‑1 (MQ=255)
                   aaTTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCa                                                                                 <  2:343011/100‑1 (MQ=255)
                      ttttGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGgc                                                                              <  2:766207/100‑1 (MQ=255)
                          gCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCtt                                                                          <  2:216470/100‑1 (MQ=255)
                           ccTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTg                                                                         <  2:1727324/100‑1 (MQ=255)
                           ccTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTg                                                                         >  2:3250274/1‑100 (MQ=255)
                            cTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGa                                                                        >  2:2305343/1‑100 (MQ=255)
                             tttttCGCGACCAGTAATGGGCGTGGAGGAACTCGTTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGa                                                                        <  1:3182518/99‑1 (MQ=255)
                             tttttCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGa                                                                        >  2:834126/1‑99 (MQ=255)
                             tttttCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAg                                                                       >  1:427186/1‑100 (MQ=255)
                             tttttCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAg                                                                       <  2:942826/100‑1 (MQ=255)
                              ttttCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGc                                                                      <  2:427189/100‑1 (MQ=255)
                              ttttCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGc                                                                      <  2:79233/100‑1 (MQ=255)
                                 tCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAg                                                                   <  1:1970717/100‑1 (MQ=255)
                                  cgcgACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGTGCTTGAGCAAGc                                                                  >  2:3083519/1‑100 (MQ=255)
                                  cgcgACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGc                                                                  <  1:3193591/100‑1 (MQ=255)
                                  cgcgACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGc                                                                  >  1:2302380/1‑100 (MQ=255)
                                       ccATCAATGGGCGTGGAGGAACTCGCTCCAGTTAAAGAGGTTCTCGAACCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCttttt                                                             <  1:3661272/100‑1 (MQ=255)
                                       ccAGCAATGGNAGNGGNGGNACTCGCTGCANTTNAAGAGGTNCTCGNATCCGGTTG‑GATCNCCACCGG‑TNCGAAGAATCAGGCGCTTGAGCAAGCatnt                                                              >  1:3010789/1‑99 (MQ=255)
                                       ccAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCttttt                                                             <  1:1007711/100‑1 (MQ=255)
                                        cAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTg                                                            <  2:3690488/100‑1 (MQ=255)
                                         aGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGc                                                           >  1:2801528/1‑100 (MQ=255)
                                         aGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGc                                                           >  2:2068587/1‑100 (MQ=255)
                                         aGCAATGGGCGTGGAGGAACTCGCCACAGCTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAATAATCAGGCGCTTGAGCAAGCTTTTTg                                                            <  1:3646659/99‑1 (MQ=255)
                                          gCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGcc                                                          <  1:1780166/100‑1 (MQ=255)
                                          gCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGcc                                                          <  2:3566710/100‑1 (MQ=255)
                                              tGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCCGGCGCTTGAGCAAGCTTTTTGCCAGtt                                                      >  1:815433/1‑100 (MQ=255)
                                              tGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGtt                                                      >  1:2179932/1‑100 (MQ=255)
                                              tGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGtt                                                      >  2:3434302/1‑100 (MQ=255)
                                                ggCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGGGCAAGCTTTTTGCCAGTTGa                                                    >  1:2031712/1‑100 (MQ=255)
                                                   gTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACgg                                                 >  2:1662751/1‑100 (MQ=255)
                                                    ttgaggaACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACggg                                                <  1:2359756/98‑1 (MQ=255)
                                                    tGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACggg                                                <  1:574955/100‑1 (MQ=255)
                                                     ggaggaACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGa                                               <  2:513937/100‑1 (MQ=255)
                                                        ggaACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAAt                                            <  2:1668359/100‑1 (MQ=255)
                                                        ggaACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAAt                                            <  2:3612496/100‑1 (MQ=255)
                                                          aaCTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCa                                          >  2:3781639/1‑100 (MQ=255)
                                                               gCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATg                                     <  2:1449006/100‑1 (MQ=255)
                                                                cTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATg                                     >  1:1728068/1‑99 (MQ=255)
                                                                 ngCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGcc                                   >  2:3506505/2‑100 (MQ=255)
                                                                   cAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCAt                                 >  2:449520/1‑100 (MQ=255)
                                                                    aGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATc                                >  1:535639/1‑100 (MQ=255)
                                                                    aGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATc                                >  2:3193896/1‑100 (MQ=255)
                                                                     gTTAAAGAGGTTCTCGAATCCGGTTG‑GATAACAACCGG‑TCCGAAGAATCAGGCGCGTGAGCAAGATTTTTGCCATTTGACGGGAAATCAGCATGCCAAcg                               <  1:3301540/100‑1 (MQ=255)
                                                                       tAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATcgcg                             >  1:1181677/1‑100 (MQ=255)
                                                                         aaGAGGTTCTCGAATCCGGTTG‑GATCAAAAACGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGt                           >  1:1651660/1‑100 (MQ=255)
                                                                              attcttcggaccGGTTGTGATC‑CAACCGGATTCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGt                       >  1:275383/11‑100 (MQ=21)
                                                                               ttcttcggaccGGTTGTGATC‑CAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTc                     >  2:304505/10‑100 (MQ=38)
                                                                               ttCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTc                     <  1:2497061/100‑1 (MQ=255)
                                                                                tcttcggaccGGTTGTGATC‑CAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATGGCGGTCAGTTc                     >  2:2344082/9‑99 (MQ=37)
                                                                                tctcGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCa                    >  1:952317/1‑100 (MQ=255)
                                                                                 cttcggaccGGTTGTGATC‑CAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAg                   >  2:2570564/8‑100 (MQ=38)
                                                                                   tcggaccGGTTGTGATC‑CAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGcc                 >  2:2254919/6‑100 (MQ=255)
                                                                                   cGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGcc                 >  2:2828731/1‑100 (MQ=255)
                                                                                   cGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGcc                 >  1:1145143/1‑100 (MQ=255)
                                                                                   cGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGcc                 <  1:2700866/100‑1 (MQ=255)
                                                                                    cggaccGGTTGTGATC‑CAACCGG‑TCCGAAGAATCCGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCa                >  2:650592/5‑100 (MQ=255)
                                                                                    cggaccGGTTGTGATC‑CAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCa                >  1:2956036/5‑100 (MQ=255)
                                                                                     ggaccGGTTGTGATC‑CAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCAc               >  2:2154599/4‑100 (MQ=255)
                                                                                        ccGGTTGTGATC‑CAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCAccgc            >  1:1660131/1‑100 (MQ=255)
                                                                                        ccGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCAccgc            >  2:2051510/1‑100 (MQ=255)
                                                                                         cGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGTGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCAccgcc           >  1:2394693/1‑100 (MQ=255)
                                                                                          ggTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCAccgccg          <  1:3561107/100‑1 (MQ=255)
                                                                                           gTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCgg         <  2:2070779/100‑1 (MQ=255)
                                                                                             tG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGaa       <  1:3778468/100‑1 (MQ=255)
                                                                                             tG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGaa       <  2:2017937/100‑1 (MQ=255)
                                                                                              g‑gATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTTCCTGTTGATTGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAAt      >  1:1569780/1‑100 (MQ=255)
                                                                                              g‑gATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAAt      <  1:3372427/100‑1 (MQ=255)
                                                                                                 aTCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAAt      <  1:687149/98‑1 (MQ=255)
                                                                                                 aTCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGc    <  2:3683949/100‑1 (MQ=255)
                                                                                                    |                                                                                                   
GAAGGAAAAGCAATGTCAGAATTTTTGCCTTTTTCGCGACCAGCAATGGGCGTGGAGGAACTCGCTGCAGTTAAAGAGGTTCTCGAATCCGGTTG‑GATCACAACCGG‑TCCGAAGAATCAGGCGCTTGAGCAAGCTTTTTGCCAGTTGACGGGAAATCAGCATGCCATCGCGGTCAGTTCAGCCACCGCCGGAATGCAT  >  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: ‑