Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,286,142 T→A 21.5% intergenic (+6/‑69) yejM → / → proL essential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor/tRNA‑Pro

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,286,1420TA21.5% 194.4 / 50.3 93intergenic (+6/‑69)yejM/proLessential inner membrane DUF3413 domain‑containing protein; lipid A production and membrane permeability factor/tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base T (39/34);  new base A (9/11);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

Only 100 of 105 total aligned reads displayed.

CAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCG  >  NC_000913/2286043‑2286240
                                                                                                   |                                                                                                  
cAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatt                                                                                                    >  2:1767197/1‑100 (MQ=255)
cAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGattaatt                                                                                                    >  2:145377/1‑100 (MQ=255)
   cgtACGTGCGGAAAGATTGAAAGGTGAAAAGCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                 <  1:1566155/98‑1 (MQ=255)
   cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACATTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                 <  1:497931/100‑1 (MQ=255)
   cTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATtata                                                                                                 >  2:2074342/1‑100 (MQ=255)
     tACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATaaa                                                                                               >  2:571089/1‑100 (MQ=255)
       cGTGGTGAAAGAGTGTAAGATGAAAAAGCACAGTTAAGATTGTTATTGCAAGTGATGACAGACGAGATGCGTTTTATCGCTAACTGATTAATTATAAATc                                                                                             >  1:1598433/1‑100 (MQ=255)
         tggtgAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAg                                                                                           >  1:285302/1‑100 (MQ=255)
          gatgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGt                                                                                          >  2:1746558/3‑100 (MQ=255)
            tgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTa                                                                                        >  1:1264811/1‑100 (MQ=255)
            tgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTa                                                                                        >  2:2455791/1‑100 (MQ=255)
            tgaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATAAGTTa                                                                                        >  1:494318/1‑100 (MQ=255)
             gaaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAg                                                                                       <  2:635791/100‑1 (MQ=38)
              aaagagtgaaagaTGAAAAACCCCAGTTCAGTTTGTTATTGCAATTGCTGAAAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAg                                                                                       <  1:1560368/99‑1 (MQ=37)
              aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTCTTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGc                                                                                      >  1:1382080/1‑100 (MQ=255)
              aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTACTCAGTTAGc                                                                                      <  2:322705/100‑1 (MQ=38)
              aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGc                                                                                      <  1:534037/100‑1 (MQ=38)
              aaagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATCAATCAGTTAGc                                                                                      <  1:1354256/100‑1 (MQ=255)
               aagagtgaaagaTGAAAAGCGACAGTTAAGTTTGTTATTGCATGTGCTGACAGACGAGAAGTGTTTTATCGCTAACTGATTTATAATTAATTAGTTAGCg                                                                                     <  1:1813965/100‑1 (MQ=37)
               aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCg                                                                                     <  2:1590911/100‑1 (MQ=255)
               aagagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGACAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAgtg                                                                                     <  2:1745452/100‑3 (MQ=38)
                 gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCGat                                                                                   <  2:667395/100‑2 (MQ=38)
                 gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCGat                                                                                   <  2:1722492/100‑2 (MQ=38)
                 gagtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTTATAATTAATCAGTTAGCGat                                                                                   <  2:666701/100‑2 (MQ=38)
                   gtgaaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAat                                                                                 >  1:626691/1‑100 (MQ=255)
                      aaagaTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATc                                                                              <  2:1302959/100‑1 (MQ=255)
                          aTGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTAc                                                                          >  2:180829/1‑100 (MQ=255)
                           tGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACt                                                                         >  2:2488184/1‑100 (MQ=255)
                           tGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACt                                                                         >  2:173224/1‑100 (MQ=255)
                             aaaaaCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTg                                                                       >  1:163851/1‑100 (MQ=255)
                                 aCCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAAt                                                                   <  1:772358/100‑1 (MQ=255)
                                      aGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgtg                                                              <  2:2116484/100‑1 (MQ=255)
                                       gTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGgtgt                                                             >  1:731965/1‑100 (MQ=255)
                                        ttAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGtgtg                                                            <  2:1019669/100‑1 (MQ=255)
                                        ttAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGtgtg                                                            <  1:919768/100‑1 (MQ=255)
                                        ttAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGtgtg                                                            <  2:607437/100‑1 (MQ=255)
                                        ttAAGATTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGtgtg                                                            <  1:721795/100‑1 (MQ=255)
                                          aaGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGa                                                          >  1:1688494/1‑100 (MQ=255)
                                            gTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGaaa                                                        >  2:318357/1‑100 (MQ=255)
                                            gTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGaaa                                                        >  2:683219/1‑100 (MQ=255)
                                            atctGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGaaa                                                        <  2:223384/97‑1 (MQ=255)
                                               tGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGTGAAATATCTTACTTGCAATCGGTGTGGAAAACg                                                     <  1:2244224/100‑1 (MQ=255)
                                                gTTATTTCAAATGCTGAAAGACGAGAAACGGTTTATCGCTAATTAATTAATTATAAATCAGTTAGCGAAATATCTTATTTGCAATCGGTGTGGAAAACgg                                                    <  1:1114597/100‑1 (MQ=255)
                                                gTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGaaaagat                                                    >  1:223382/1‑97 (MQ=255)
                                                cgatctGCAAGTGCTGACAGACGGGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACgg                                                    <  2:2270/95‑1 (MQ=255)
                                                 ttATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                   <  1:1767063/100‑1 (MQ=255)
                                                 ttATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgt                                                   <  1:755560/100‑1 (MQ=255)
                                                   aTTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtag                                                 <  1:1342135/100‑1 (MQ=255)
                                                   aTTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtag                                                 <  1:1972535/100‑1 (MQ=255)
                                                   aTTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtag                                                 <  1:2486758/100‑1 (MQ=255)
                                                    cgatctgTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtagt                                                <  1:44195/94‑1 (MQ=255)
                                                     tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtagta                                               >  2:1800347/1‑100 (MQ=255)
                                                     tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtagta                                               >  2:2231889/1‑100 (MQ=255)
                                                     tGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtagta                                               >  2:1262338/1‑100 (MQ=255)
                                                      gCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACggagatcg                                              >  1:2270/1‑94 (MQ=255)
                                                      gCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAt                                              >  1:1058776/1‑100 (MQ=255)
                                                      gCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTAt                                              >  1:1291186/1‑100 (MQ=255)
                                                        aaGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTa                                            <  2:2086975/100‑1 (MQ=255)
                                                         aGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTag                                           <  1:1605166/100‑1 (MQ=255)
                                                          gTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGgtagtagatcg                                          >  2:44196/1‑95 (MQ=255)
                                                           tGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagca                                         <  2:580864/100‑1 (MQ=255)
                                                             cTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTagcagc                                       >  1:200721/1‑100 (MQ=255)
                                                              tGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGcc                                      >  1:22591/1‑100 (MQ=255)
                                                                 cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                   <  2:1232315/100‑1 (MQ=255)
                                                                 cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                   >  2:678110/1‑100 (MQ=255)
                                                                 cAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACg                                   >  1:1571953/1‑100 (MQ=255)
                                                                   gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGa                                  >  1:2273784/1‑99 (MQ=255)
                                                                   gACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAg                                 <  2:17881/100‑1 (MQ=255)
                                                                      gagaAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCg                              <  2:2238470/100‑1 (MQ=255)
                                                                       agaAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTc                               >  1:2458755/1‑98 (MQ=255)
                                                                        gaACCGTTTTAACGCTAACTGAGTAATTATAAGTCAGTCAGCGAAATATCTTACTTGCAATCGGTTTGGAAAACGATAGTATTAGCAGCCAccattc                               >  1:1030548/1‑92 (MQ=255)
                                                                         aaGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCa                           <  1:508481/100‑1 (MQ=255)
                                                                         aaGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCa                           <  2:2271376/100‑1 (MQ=255)
                                                                          aGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCAc                          >  2:445025/1‑100 (MQ=255)
                                                                           gCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACg                         <  1:2285856/100‑1 (MQ=255)
                                                                             gTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTa                       <  1:298268/100‑1 (MQ=255)
                                                                               tttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgc                     <  2:200723/100‑1 (MQ=255)
                                                                               tttATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgc                     >  2:32277/1‑100 (MQ=255)
                                                                                 tATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAgcgc                   >  2:673339/1‑100 (MQ=255)
                                                                                   tCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATATTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAg                 >  1:1594947/1‑100 (MQ=255)
                                                                                     gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc               >  2:1082021/1‑100 (MQ=38)
                                                                                     gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc               >  1:288284/1‑100 (MQ=38)
                                                                                     gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGc                >  1:1456953/1‑99 (MQ=38)
                                                                                     gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGGCGGCACGTAGCGCAGcc               >  2:2097269/1‑100 (MQ=38)
                                                                                     gCTAACTGATTTATAATTAATCAGTTAGCGAAATATCCTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGACGGCTCGTAGCGCAGcc               >  1:924871/1‑100 (MQ=37)
                                                                                     gCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGcc               <  1:1341021/100‑1 (MQ=255)
                                                                                      cTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTCGCGCAGCCt              >  1:1018672/1‑100 (MQ=38)
                                                                                      cTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt              >  2:2531697/1‑100 (MQ=38)
                                                                                      cTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt              >  1:752847/1‑100 (MQ=255)
                                                                                      cTAACTGATTTATAATTAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt              >  1:1309901/1‑100 (MQ=38)
                                                                                      cTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCt              <  2:2176927/100‑1 (MQ=255)
                                                                                         aCTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGCGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGt           <  1:2560073/100‑1 (MQ=255)
                                                                                          cTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTa          <  1:2575041/100‑1 (MQ=255)
                                                                                            gATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgc        <  1:2457767/100‑1 (MQ=255)
                                                                                             attaattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAgcg       >  1:2430829/1‑100 (MQ=255)
                                                                                               taattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCa     >  2:378293/1‑100 (MQ=255)
                                                                                                aattaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAc    >  2:1901354/1‑100 (MQ=255)
                                                                                                 attaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCAcc   >  2:638522/1‑100 (MQ=255)
                                                                                                  ttaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCg  <  1:1746549/100‑1 (MQ=255)
                                                                                                  ttaTAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCg  <  2:694227/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
CAACTTACGTGGTGAAAGAGTGAAAGATGAAAAACCACAGTTAAGTTTGTTATTGCAAGTGCTGACAGACGAGAAGCGTTTTATCGCTAACTGATTAATTATAAATCAGTTAGCGAAATATCTTACTTGCAATCGGTGTGGAAAACGGTAGTATTAGCAGCCACGAGTCGGCACGTAGCGCAGCCTGGTAGCGCACCG  >  NC_000913/2286043‑2286240

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: