Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,015,425 A→T L90Q (CTG→CAG)  fabA ← beta‑hydroxydecanoyl thioester dehydrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,015,4250AT99.4% 665.3 / ‑3.5 178L90Q (CTG→CAG) fabAbeta‑hydroxydecanoyl thioester dehydrase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base T (88/89);  total (89/89)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.13e-01

Only 100 of 188 total aligned reads displayed.

GCAGTACCTGACCAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCAGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATCCGGATTGAT  >  NC_000913/1015336‑1015513
                                                                                         |                                                                                        
gCAGTACCTGACCAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCt                                                                                       >  1:5347022/1‑93 (MQ=255)
 cAGTACCTGACCAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACctgctg                                                                                      >  2:5379037/1‑93 (MQ=255)
  aGTACCTGACCAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCtgctgc                                                                                     >  1:3391066/1‑93 (MQ=255)
  aGTACCTGACCAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCtgctgc                                                                                     >  1:4012738/1‑93 (MQ=255)
  aGTACCTGACCAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCtgctgc                                                                                     <  1:5531179/93‑1 (MQ=255)
   gTACCTGACCAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGcc                                                                                    >  2:4412243/1‑93 (MQ=255)
    tACCTGACCAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCca                                                                                   <  1:4748991/93‑1 (MQ=255)
    tACCTGACCAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCca                                                                                   >  1:3375262/1‑93 (MQ=255)
    tACCTGACCAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCca                                                                                   >  2:2763538/1‑93 (MQ=255)
      ccTGACCAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCcaca                                                                                 <  2:3669134/93‑1 (MQ=255)
      ccTGACCAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCcaca                                                                                 <  2:5818603/93‑1 (MQ=255)
       cTGACCAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACAt                                                                                >  1:3012594/1‑93 (MQ=255)
           ccAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCg                                                                            >  2:3499940/1‑93 (MQ=255)
           ccAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCg                                                                            >  2:5255368/1‑93 (MQ=255)
           ccAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCg                                                                            <  2:831514/93‑1 (MQ=255)
            cAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGGAGAACCCTACCTGCTGCCACATTGCGt                                                                           >  2:3102640/1‑93 (MQ=255)
              gTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTcc                                                                         >  2:4672376/1‑93 (MQ=255)
               tGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTcca                                                                        <  1:864652/93‑1 (MQ=255)
               tGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTcca                                                                        >  1:4907352/1‑93 (MQ=255)
               tGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTcca                                                                        <  2:2210386/93‑1 (MQ=255)
               tGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTcca                                                                        >  2:1836520/1‑93 (MQ=255)
               tGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTcca                                                                        <  2:449457/93‑1 (MQ=255)
               tGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTcca                                                                        >  1:108648/1‑93 (MQ=255)
                gAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTccag                                                                       <  1:4943693/93‑1 (MQ=255)
                 aaTTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTccagg                                                                      >  1:1222293/1‑93 (MQ=255)
                  aTTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTccaggc                                                                     <  2:1700505/93‑1 (MQ=255)
                   tttCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTccaggcc                                                                    >  1:5634676/1‑93 (MQ=255)
                   tttCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTccaggcc                                                                    <  2:2450939/93‑1 (MQ=255)
                     tCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTccaggccca                                                                  >  2:1585706/1‑93 (MQ=255)
                     tCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTccaggccca                                                                  <  1:792232/93‑1 (MQ=255)
                     tCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTccaggccca                                                                  >  1:1095551/1‑93 (MQ=255)
                     tCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTccaggccca                                                                  >  1:1851489/1‑93 (MQ=255)
                      cACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGTGGTAGAACCCTACCTGCTGCCACATTGCGTccaggcccag                                                                 <  2:1370258/93‑1 (MQ=255)
                      cACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTccaggcccag                                                                 >  2:4751402/1‑93 (MQ=255)
                       aCTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTccaggcccagg                                                                <  2:4163670/93‑1 (MQ=255)
                         ttCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCa                                                              <  2:3244394/93‑1 (MQ=255)
                             ccAACGCCCAGCGCGCGGCTTTCACCTTCGCCTCCGAGCCAGCCGAGGTAGAACCCTCCCTGCTGCCACATTGCGTCCAGGCCCAGGCATccc                                                          <  1:2033045/93‑1 (MQ=255)
                              cAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCg                                                         >  2:3197921/1‑93 (MQ=255)
                               aaCGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCgg                                                        >  2:4777841/1‑93 (MQ=255)
                               aaCGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCgg                                                        >  1:2106833/1‑93 (MQ=255)
                               aaCGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCgg                                                        >  1:5614905/1‑93 (MQ=255)
                               aaCGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCgg                                                        >  1:6104209/1‑93 (MQ=255)
                                 cGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCa                                                      <  1:2124279/93‑1 (MQ=255)
                                  gCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCAt                                                     >  2:2460139/1‑93 (MQ=255)
                                  gCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCAt                                                     >  2:3429541/1‑93 (MQ=255)
                                  gCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATACCGGCAt                                                     >  1:4967932/1‑93 (MQ=255)
                                    ccAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATaa                                                   <  1:4274158/93‑1 (MQ=255)
                                     caatgcgcgGCATTTCCCTACTCCAACTAGCCAGCTGAGGTAGAACCCTACCTGCTGCCGCATTGCGTCCAGGCCCAGGCATCCCGGCATAAc                                                  <  2:2256766/89‑1 (MQ=255)
                                     cAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAAc                                                  >  1:1750629/1‑93 (MQ=255)
                                     cAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAAc                                                  <  1:2437229/93‑1 (MQ=255)
                                       gcgcgcGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCg                                                <  2:2972379/93‑1 (MQ=255)
                                        cgcgcgGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCAAGGCATCCCGGCATAACAgg                                               >  1:4617225/1‑93 (MQ=255)
                                         gcgcgGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGa                                              <  2:4967300/93‑1 (MQ=255)
                                          cgcgGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGAt                                             <  1:2842922/93‑1 (MQ=255)
                                             ggCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGc                                          >  1:237711/1‑93 (MQ=255)
                                              gCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGcc                                         >  2:4382850/1‑93 (MQ=255)
                                              gCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGcc                                         >  2:2658083/1‑93 (MQ=255)
                                              gCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGcc                                         <  1:4790384/93‑1 (MQ=255)
                                               ccTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCa                                        >  2:1663627/1‑93 (MQ=255)
                                               ccTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCa                                        >  1:1990015/1‑93 (MQ=255)
                                                 tttACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAAt                                      <  1:4586695/93‑1 (MQ=255)
                                                  ttACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATa                                     <  2:4479657/93‑1 (MQ=255)
                                                    aCCTTCGCCGCACAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATaaa                                   <  1:3833349/93‑1 (MQ=255)
                                                     ccTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAg                                  <  1:3831306/93‑1 (MQ=255)
                                                     ccTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAg                                  <  1:1996338/93‑1 (MQ=255)
                                                     ccTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCAGCATAACCGGATCGCCAATaacg                                  >  1:4492364/1‑91 (MQ=255)
                                                      cTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAGAGt                                 >  2:1090821/1‑93 (MQ=255)
                                                      cTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGt                                 <  1:998877/93‑1 (MQ=255)
                                                      cTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGt                                 >  1:957790/1‑93 (MQ=255)
                                                      cTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGt                                 <  1:1414547/93‑1 (MQ=255)
                                                          gccaccCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCa                             <  1:3203535/93‑1 (MQ=255)
                                                             gccCAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATcc                          <  2:3591623/93‑1 (MQ=255)
                                                             ccccAGCCAGCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATcc                          <  1:5140614/92‑1 (MQ=255)
                                                               ccagccagCCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCga                        >  2:5740674/1‑93 (MQ=255)
                                                                cagccagcCGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCgaa                       >  1:1661327/1‑93 (MQ=255)
                                                                 agccagccGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCgaag                      <  1:5753561/93‑1 (MQ=255)
                                                                 agccagccGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCgaag                      <  1:5244347/93‑1 (MQ=255)
                                                                    cagccGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAAc                   <  2:1391099/93‑1 (MQ=255)
                                                                      gccGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACca                 >  2:2891140/1‑93 (MQ=255)
                                                                       ccGAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACcac                <  1:2509300/93‑1 (MQ=255)
                                                                         gAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACcacat              >  2:64322/1‑92 (MQ=255)
                                                                         gAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAg              >  2:5338099/1‑93 (MQ=255)
                                                                         gAGGTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAg              <  2:3556226/93‑1 (MQ=255)
                                                                           ggTAGAACCCTACCAGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGAt            >  1:2236305/1‑93 (MQ=255)
                                                                           agTAGAACCCTACCTGCTGCCACATTGCGCCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGAt            <  2:2136953/92‑1 (MQ=255)
                                                                            gTAGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATc           >  1:2855063/1‑93 (MQ=255)
                                                                              aGACCCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATCCg         <  1:3441604/93‑1 (MQ=255)
                                                                              aGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATCCg         <  2:4687788/93‑1 (MQ=255)
                                                                              aGAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATCCg         >  2:5542784/1‑93 (MQ=255)
                                                                               gAACCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATCCgg        <  2:1719884/93‑1 (MQ=255)
                                                                                aaCCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATCCGGa       >  2:5624341/1‑93 (MQ=255)
                                                                                 aCCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATCCGGAt      <  2:4881094/93‑1 (MQ=255)
                                                                                 aCCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATCCGGAt      >  2:376996/1‑93 (MQ=255)
                                                                                 aCCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATCCGGAt      >  1:4703813/1‑93 (MQ=255)
                                                                                 aCCCTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATCCGGAt      >  2:723000/1‑93 (MQ=255)
                                                                                  cccTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATCCGGAtt     <  2:3409847/93‑1 (MQ=255)
                                                                                  cccTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATCCGGAtt     >  2:3473994/1‑93 (MQ=255)
                                                                                  cccTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCCGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATCCGCAtt     >  1:5905983/1‑93 (MQ=255)
                                                                                    cTACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATCCGGATTGa   >  2:3939339/1‑93 (MQ=255)
                                                                                     tACCTGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATCCGGATTGat  >  2:744820/1‑93 (MQ=255)
                                                                                         |                                                                                        
GCAGTACCTGACCAGTGAATTTCACTTCGCCAACGCCCAGCGCGCGGCCTTTACCTTCGCCGCCCAGCCAGCCGAGGTAGAACCCTACCAGCTGCCACATTGCGTCCAGGCCCAGGCATCCCGGCATAACCGGATCGCCAATAAAGTGGCATCCGAAGAACCACAGATCCGGATTGAT  >  NC_000913/1015336‑1015513

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