Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,099,661 (A)5→6 intergenic (‑19/+164) yggN ← / ← b2959 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,099,6561.A99.0% 688.8 / ‑1.6 194intergenic (‑14/+169)yggN/b2959hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (88/104);  total (89/105)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 214 total aligned reads displayed.

TGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGC  >  NC_000913/3099564‑3099743
                                                                                             |                                                                                        
 gACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGa                                                                                         <  2:2584823/93‑1 (MQ=255)
 gACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGa                                                                                         <  1:1699542/93‑1 (MQ=255)
    gCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaa                                                                                     <  2:537835/93‑1 (MQ=255)
      tGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaa                                                                                   <  1:4091542/93‑2 (MQ=255)
      tGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaa                                                                                   >  1:655729/1‑92 (MQ=255)
      tGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaa                                                                                   >  2:3784625/1‑92 (MQ=255)
        cACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaatc                                                                                 <  1:3757350/93‑4 (MQ=255)
         aCTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGaaaaaatca                                                                                <  2:3698549/93‑5 (MQ=255)
           tGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAgg                                                                               >  1:3946242/1‑93 (MQ=255)
           tGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAgg                                                                               >  1:3019104/1‑93 (MQ=255)
             gTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGcc                                                                             <  1:2921886/93‑1 (MQ=255)
             gTAGTCGGCGTGAGCGCTCATTGCCGTCCCTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGcc                                                                             <  2:1293281/93‑1 (MQ=255)
              tAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCa                                                                            <  2:2089088/93‑1 (MQ=255)
               aGTCGTCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAg                                                                           <  2:3246228/93‑1 (MQ=255)
                 tCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCa                                                                         <  1:1954855/93‑1 (MQ=255)
                   ggCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAtt                                                                       <  1:1616728/93‑1 (MQ=255)
                    gCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCAttt                                                                      >  2:1815703/1‑93 (MQ=255)
                     cGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTg                                                                     <  2:637581/93‑1 (MQ=255)
                      gTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGc                                                                    >  2:1341102/1‑93 (MQ=255)
                      gTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGc                                                                    <  2:448630/93‑1 (MQ=255)
                       tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCt                                                                   >  2:1557303/1‑93 (MQ=255)
                       tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCt                                                                   <  1:2653813/93‑1 (MQ=255)
                       tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCt                                                                   <  2:1065870/93‑1 (MQ=255)
                       tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCt                                                                   >  2:103009/1‑93 (MQ=255)
                       tGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCt                                                                   >  1:1194576/1‑93 (MQ=255)
                        gAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTg                                                                  >  2:1258366/1‑93 (MQ=255)
                        gAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTg                                                                  <  2:4268728/93‑1 (MQ=255)
                         aGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTgg                                                                 <  2:3150850/93‑1 (MQ=255)
                         aGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTgc                                                                 >  1:1548090/1‑92 (MQ=255)
                            ggTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGccc                                                              >  2:4056663/1‑93 (MQ=255)
                            ggTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGccc                                                              >  2:4056665/1‑93 (MQ=255)
                             gTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGcccc                                                             <  1:3792357/93‑1 (MQ=255)
                             gTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGcccc                                                             <  1:3409539/93‑1 (MQ=255)
                              tCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCa                                                            <  1:3632389/93‑1 (MQ=255)
                              tCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCa                                                            <  1:176948/93‑1 (MQ=255)
                               cATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAg                                                           <  1:2238315/93‑1 (MQ=255)
                                aTTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGa                                                          >  2:2552693/1‑93 (MQ=255)
                                aTTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGa                                                          <  2:2250109/93‑1 (MQ=255)
                                aTTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGa                                                          <  2:3971643/93‑1 (MQ=255)
                                 ttGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAt                                                         <  2:2894486/93‑1 (MQ=255)
                                  tGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGAtt                                                        <  2:4083466/93‑1 (MQ=255)
                                    ccGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                                                      <  2:3436923/93‑1 (MQ=255)
                                    ccGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGa                                                      >  2:3776687/1‑93 (MQ=255)
                                      gTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATa                                                    <  2:4309505/93‑1 (MQ=255)
                                       tCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATaa                                                   >  1:5686/1‑93 (MQ=255)
                                       tCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATaa                                                   >  1:1902370/1‑93 (MQ=255)
                                         aCTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACa                                                 >  1:4479851/1‑93 (MQ=255)
                                          cTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaa                                                >  1:4674989/1‑93 (MQ=255)
                                           tGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACaaa                                               <  1:3557043/93‑1 (MQ=255)
                                            gAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg                                              <  2:1082159/93‑1 (MQ=255)
                                            gAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg                                              >  2:4114354/1‑93 (MQ=255)
                                            gAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAg                                              <  2:3108133/93‑1 (MQ=255)
                                             aaaGTGCTGCCGCCAGCAGCATGTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTAATAACAAAGt                                             <  1:3805353/93‑1 (MQ=255)
                                              aaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTa                                            <  1:3358434/93‑1 (MQ=255)
                                              aaGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTa                                            >  1:3326602/1‑93 (MQ=255)
                                               aGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAc                                           <  1:3150568/93‑1 (MQ=255)
                                               aGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAc                                           <  2:2403373/93‑1 (MQ=255)
                                               aGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAc                                           >  2:214406/1‑93 (MQ=255)
                                               aGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAc                                           <  2:2039181/93‑1 (MQ=255)
                                                gTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcg                                          <  1:3211821/93‑1 (MQ=255)
                                                 tgctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTAcgc                                         <  1:2662969/93‑1 (MQ=255)
                                                   ctgcCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGt                                       <  2:265267/93‑1 (MQ=255)
                                                     gccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGttt                                     <  2:1741231/93‑1 (MQ=255)
                                                      ccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGtttt                                    <  1:3998283/93‑1 (MQ=255)
                                                      ccgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGtttt                                    <  2:4673239/93‑1 (MQ=255)
                                                       cgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTc                                   >  2:623685/1‑93 (MQ=255)
                                                       cgccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTc                                   <  1:4452715/93‑1 (MQ=255)
                                                        gccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTcc                                  >  2:2992609/1‑93 (MQ=255)
                                                        gccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTcc                                  >  2:2500036/1‑93 (MQ=255)
                                                         ccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCa                                 <  1:2631513/93‑1 (MQ=255)
                                                         ccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCa                                 <  2:3582001/93‑1 (MQ=255)
                                                         ccAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCa                                 <  1:2586261/93‑1 (MQ=255)
                                                          cagcagCATTTTGCGCTTCATAGTCTTCCCTCAGGAAACAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCAt                                <  2:4200835/93‑1 (MQ=255)
                                                            gcagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGc                              >  1:1400394/1‑93 (MQ=255)
                                                            gcagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGc                              >  1:3283146/1‑93 (MQ=255)
                                                             cagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGcc                             <  1:1153573/93‑1 (MQ=255)
                                                             cagcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGcc                             <  2:1274488/93‑1 (MQ=255)
                                                              agcaTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCg                            >  1:2711693/1‑93 (MQ=255)
                                                               gcaTTTTGCGCATAATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCgg                           <  1:1005115/93‑1 (MQ=255)
                                                                caTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGa                          >  2:1608768/1‑93 (MQ=255)
                                                                 aTTTTGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGAt                         >  2:3869176/1‑93 (MQ=255)
                                                                  ttttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATg                        <  1:4208881/93‑1 (MQ=255)
                                                                    ttGCGCATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcg                      <  1:981684/93‑1 (MQ=255)
                                                                       cgcATCATAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATgcggcg                   >  1:242225/1‑93 (MQ=255)
                                                                          atcatAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTaa                >  2:803330/1‑93 (MQ=255)
                                                                            catAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAac              >  1:917110/1‑93 (MQ=255)
                                                                             atAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAaca             <  1:2881923/93‑1 (MQ=255)
                                                                              tAGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAacac            >  2:3881773/1‑93 (MQ=255)
                                                                               aGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACAcc           >  1:502620/1‑93 (MQ=255)
                                                                               aGTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACAcc           >  2:485457/1‑93 (MQ=255)
                                                                                gTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCt          >  2:1948870/1‑93 (MQ=255)
                                                                                gTCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCt          >  1:2214329/1‑93 (MQ=255)
                                                                                gTCTTCACTCAAGAAAAAATCAGGCCCGCATTTTCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCt          <  1:3804997/93‑1 (MQ=37)
                                                                                 tCTTCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCtt         <  2:3209195/93‑1 (MQ=255)
                                                                                   ttCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTAt       >  1:732647/1‑93 (MQ=255)
                                                                                   ttCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTAt       >  1:2532498/1‑93 (MQ=255)
                                                                                   ttCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTAt       >  2:1980931/1‑93 (MQ=255)
                                                                                   ttCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCAGGATGCGGCGTAAACACCTTAt       >  2:1385475/1‑93 (MQ=255)
                                                                                    tCCCTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATc      >  1:2759315/1‑93 (MQ=255)
                                                                                     cccTCAAGAAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATcc     <  2:2382219/93‑1 (MQ=255)
                                                                                             |                                                                                        
TGACGCTGCACTGGTAGTCGGCGTGAGCGGTCATTGCCGTCACTGAAAGTGCTGCCGCCAGCAGCATTTTGCGCATCATAGTCTTCCCTCAAGAAAAATCAGGCCAGCATTTGCTGGCCCCAGATTGATAACAAAGTACGCGTTTTCCATGCCGGATGCGGCGTAAACACCTTATCCCGC  >  NC_000913/3099564‑3099743

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