Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,885,388 A→G intergenic (‑146/+213) ygcB ← / ← cysH hypothetical protein/3'‑phosphoadenosine 5'‑phosphosulfate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,885,3880AG99.4% 497.5 / ‑4.7 155intergenic (‑146/+213)ygcB/cysHhypothetical protein/3'‑phosphoadenosine 5'‑phosphosulfate reductase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (59/95);  total (59/96)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.00e-01

Only 100 of 159 total aligned reads displayed.

ATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACC  >  NC_000913/2885298‑2885477
                                                                                          |                                                                                         
aTGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCg                                                                                         >  2:1659244/1‑93 (MQ=255)
 atctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAgcgc                                                                                        <  2:3086699/91‑1 (MQ=255)
  gctgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCt                                                                                       >  1:4560045/1‑93 (MQ=255)
    tgctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTgg                                                                                     <  1:4527088/93‑1 (MQ=255)
      ctCATAAATATCGACTTACAGTGATGTACACATGTCTTGCCTCCCTGCAAAGTAGAGCCTAACATTCGAATTGTTTGGGTTCGAGCGCTGGcc                                                                                   <  1:3615447/93‑1 (MQ=255)
      ctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGcc                                                                                   <  2:1438885/93‑1 (MQ=255)
      ctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGcc                                                                                   >  1:3990081/1‑93 (MQ=255)
      ctCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGcc                                                                                   <  1:3202330/93‑1 (MQ=255)
       tCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCt                                                                                  >  2:3400040/1‑93 (MQ=255)
       tCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCt                                                                                  <  2:3451987/93‑1 (MQ=255)
       tCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCt                                                                                  >  1:529800/1‑93 (MQ=255)
       tCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCt                                                                                  >  1:676697/1‑93 (MQ=255)
        cATAAATATCGAATTACAGTGATATACACACGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTc                                                                                 >  1:3723433/1‑93 (MQ=255)
         aTAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCa                                                                                <  2:2763862/93‑1 (MQ=255)
         aTAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCa                                                                                <  2:2595169/93‑1 (MQ=255)
          tAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAg                                                                               <  2:9467/93‑1 (MQ=255)
          tAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAg                                                                               <  1:4609184/93‑1 (MQ=255)
            aaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGt                                                                             <  2:1452213/93‑1 (MQ=255)
            aaTATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGt                                                                             >  1:4654926/1‑93 (MQ=255)
              tatCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTg                                                                           >  2:2277724/1‑93 (MQ=255)
              tatCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTg                                                                           >  2:1684243/1‑93 (MQ=255)
               atCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGa                                                                          <  2:582830/93‑1 (MQ=255)
                 cGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATa                                                                        <  2:4463630/93‑1 (MQ=255)
                 cGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATa                                                                        <  1:2499722/93‑1 (MQ=255)
                  gAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAg                                                                       >  1:2976629/1‑93 (MQ=255)
                  gAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAg                                                                       <  2:4663761/93‑1 (MQ=255)
                    aTTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGaa                                                                     >  2:1485927/1‑93 (MQ=255)
                    aTTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGaa                                                                     <  1:161717/93‑1 (MQ=255)
                    aTTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGaa                                                                     >  2:645079/1‑93 (MQ=255)
                     ttACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGaaa                                                                    >  1:2664625/1‑93 (MQ=255)
                     atACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGaaa                                                                    <  1:1374651/92‑1 (MQ=255)
                        cAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAatat                                                                 >  1:1709902/1‑93 (MQ=255)
                         aatGATATACCAATGGCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATc                                                                <  1:3762852/91‑1 (MQ=255)
                         aGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATc                                                                <  2:1210498/93‑1 (MQ=255)
                          gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCg                                                               <  2:3105273/93‑1 (MQ=255)
                          gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCg                                                               >  1:4717953/1‑93 (MQ=255)
                          gTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCg                                                               <  1:2536029/93‑1 (MQ=255)
                           tGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGc                                                              <  2:1021677/93‑1 (MQ=255)
                             atatACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCt                                                            <  2:3315884/93‑1 (MQ=255)
                             atatACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCt                                                            <  2:3472893/93‑1 (MQ=255)
                             atatACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCt                                                            <  1:2862908/93‑1 (MQ=255)
                               ataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTgg                                                          >  2:3193132/1‑93 (MQ=255)
                               ataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTgg                                                          >  2:1051088/1‑93 (MQ=255)
                               ataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTgg                                                          <  1:98307/93‑1 (MQ=255)
                               ataCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTgg                                                          >  2:4437789/1‑93 (MQ=255)
                                taCACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTggg                                                         <  2:2869808/93‑1 (MQ=255)
                                   acaTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCt                                                      <  2:1449820/93‑1 (MQ=255)
                                    caTGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCtt                                                     >  1:558738/1‑93 (MQ=255)
                                        tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGt                                                 <  1:1245871/93‑1 (MQ=255)
                                        tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGt                                                 >  2:1847871/1‑93 (MQ=255)
                                        tCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGt                                                 <  2:874060/93‑1 (MQ=255)
                                         cTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTc                                                <  1:599098/93‑1 (MQ=255)
                                          ttGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTcc                                               >  2:4136323/1‑93 (MQ=255)
                                          ttGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTcc                                               <  1:2939258/93‑1 (MQ=255)
                                             ccTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATc                                            <  1:4640359/93‑1 (MQ=255)
                                             ccTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATc                                            <  1:1976206/93‑1 (MQ=255)
                                             ccTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATc                                            >  2:1996369/1‑93 (MQ=255)
                                             ccTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATc                                            <  2:2084022/93‑1 (MQ=255)
                                                  cTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGaa                                       >  2:336509/1‑93 (MQ=255)
                                                  cTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGaa                                       >  2:420559/1‑93 (MQ=255)
                                                   tGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAc                                      >  2:995583/1‑93 (MQ=255)
                                                   tGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAAc                                      >  1:3219273/1‑93 (MQ=255)
                                                      aaaGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTc                                   <  2:1328829/93‑1 (MQ=255)
                                                      aaaGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTc                                   <  1:1829417/93‑1 (MQ=255)
                                                      aaaGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTc                                   >  1:2573020/1‑93 (MQ=255)
                                                      aaaGAAGAGACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGACTCAGGGTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTc                                   >  1:2573019/1‑93 (MQ=255)
                                                           agagACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGaa                              <  2:3349660/93‑1 (MQ=255)
                                                           agagACTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGaa                              <  2:2682174/93‑1 (MQ=255)
                                                              gaCTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATga                           >  2:2752037/1‑92 (MQ=255)
                                                               aCTAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCt                          >  1:3446266/1‑93 (MQ=255)
                                                                 tAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTa                        >  1:604531/1‑93 (MQ=255)
                                                                 tAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTa                        <  2:3841408/93‑1 (MQ=255)
                                                                 tAACATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTa                        >  2:4323050/1‑93 (MQ=255)
                                                                  aaCATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTaa                       <  1:796289/93‑1 (MQ=255)
                                                                  aaCATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTaa                       <  2:3239567/93‑1 (MQ=255)
                                                                  aaCATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTaa                       <  2:1306703/93‑1 (MQ=255)
                                                                  aaCATTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTaa                       <  1:3788797/93‑1 (MQ=255)
                                                                     aTTCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGc                    <  1:3397108/93‑1 (MQ=255)
                                                                     aTTCGAATTGTTTGGGTTCGAGCGCTGGCATCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGc                    <  1:4203313/93‑1 (MQ=255)
                                                                      ttCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc                   <  2:2667408/93‑1 (MQ=255)
                                                                      ttCGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGcc                   <  1:1277246/93‑1 (MQ=255)
                                                                        cGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAg                 >  2:3509498/1‑93 (MQ=255)
                                                                        cGAATTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAg                 >  1:108771/1‑93 (MQ=255)
                                                                          aaTTTTTTTGGTTCGAGCGCTGGCCTCAAGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCCCGCGAAAGCTTAACCCCAGAc               <  2:4260473/93‑1 (MQ=255)
                                                                          aaTTGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGAc               <  1:3858801/93‑1 (MQ=255)
                                                                           aTTGTTTGGGTGCGAGCGCTGTCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACa              <  2:1837246/93‑1 (MQ=255)
                                                                            ttGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAg             >  1:1543582/1‑93 (MQ=255)
                                                                             tGTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGc            <  2:2851048/93‑1 (MQ=255)
                                                                              gTTTGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGcc           <  1:31101/93‑1 (MQ=255)
                                                                               tttGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCt          >  2:2385187/1‑93 (MQ=255)
                                                                                 tGGGTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCa        >  1:4070981/1‑93 (MQ=255)
                                                                                  gggTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCaa       <  2:2930394/93‑1 (MQ=255)
                                                                                  gggTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCaa       <  2:1225616/93‑1 (MQ=255)
                                                                                   ggTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAg      <  1:2242090/93‑1 (MQ=255)
                                                                                    gTTCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGc     >  1:4181976/1‑93 (MQ=255)
                                                                                    gTTCGAGCGCTGGCCTCAGGTTAATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGc     <  2:2260414/93‑1 (MQ=255)
                                                                                     ttCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCa    >  2:3971562/1‑93 (MQ=255)
                                                                                     ttCGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCa    <  2:4622767/93‑1 (MQ=255)
                                                                                       cGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCAcc  <  1:4420148/93‑1 (MQ=255)
                                                                                       cGAGCGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCAcc  <  2:779847/93‑1 (MQ=255)
                                                                                          |                                                                                         
ATGCTGCTCATAAATATCGAATTACAGTGATATACACATGTCTTGCCTCCCTGCAAAGAAGAGACTAACATTCGAATTGTTTGGGTTCGAACGCTGGCCTCAGGTTGATAGAAATATCGCCTGGGGCTTTTGTCCATCTGGAACCTCGCGAATGCTTAACGCCAGACAGCCTCAAGCACC  >  NC_000913/2885298‑2885477

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