Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,020,958 A→T 15.8% intergenic (‑39/‑180) sulA ← / → sxy SOS cell division inhibitor/CRP‑S‑dependent promoter expression factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,020,9580AT15.8% 339.6 / 50.0 133intergenic (‑39/‑180)sulA/sxySOS cell division inhibitor/CRP‑S‑dependent promoter expression factor
Reads supporting (aligned to +/- strand):  ref base A (61/51);  new base T (12/9);  total (73/60)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.69e-01

Only 100 of 142 total aligned reads displayed.

AATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTT  >  NC_000913/1020860‑1021055
                                                                                                  |                                                                                                 
aaTTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacag                                                                                                  <  1:76335/100‑1 (MQ=255)
aaTTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacag                                                                                                  <  1:76378/100‑1 (MQ=255)
   tttACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagtac                                                                                               >  1:4071417/1‑100 (MQ=255)
   tttACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATgtacagtac                                                                                               >  2:3861478/1‑100 (MQ=255)
       ctttcttctgATTAGAACGACAAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATTGATGTACAGTACATcc                                                                                           <  2:3715209/91‑1 (MQ=255)
       cTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATcc                                                                                           <  2:3991860/100‑1 (MQ=255)
         tgctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAg                                                                                         <  1:4111325/100‑1 (MQ=255)
         tgctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAg                                                                                         >  2:3324737/1‑100 (MQ=255)
          gctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTGCATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGt                                                                                        >  2:2326018/1‑100 (MQ=255)
           ctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCatac                                                                                       <  2:4012931/100‑4 (MQ=255)
           ctgcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCatac                                                                                       <  2:2453538/100‑4 (MQ=255)
             gcGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCatacag                                                                                     <  1:13987/100‑6 (MQ=255)
              cGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGaca                                                                                    >  1:1033863/1‑100 (MQ=255)
                gATGAGAACGACGAAGAACGATGTTCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGacaac                                                                                  <  2:2763763/100‑1 (MQ=255)
                 aTGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCatacagtaac                                                                                 <  2:4129591/100‑10 (MQ=39)
                  tGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGGGTTACTGTATGGATGTACAGTACACCCAGTGAcaacaa                                                                                <  2:2588634/100‑1 (MQ=255)
                  tGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACTGTACATCCatacagtaact                                                                                <  1:2569167/100‑11 (MQ=38)
                  tGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGAcaacaa                                                                                >  1:210129/1‑100 (MQ=255)
                  tGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGAcaacaa                                                                                <  1:3096198/100‑1 (MQ=255)
                   gagaACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACaaa                                                                               <  2:3602863/100‑1 (MQ=255)
                    agaACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAg                                                                              >  1:3007558/1‑100 (MQ=255)
                     gaACGACGAAGAACGATGTGCATAGCCTGCAGTGTCCATAATCAATCCAGCCCCTGTGCGTTACTGTAGGGATGTACAGTACATCCCGTGACAACAAGGa                                                                             >  2:3718290/1‑100 (MQ=255)
                     gaACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGa                                                                             <  1:568814/100‑1 (MQ=255)
                       acgacgAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATc                                                                           >  2:1441889/1‑100 (MQ=255)
                        cgacgaAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCa                                                                          >  2:1128913/1‑100 (MQ=255)
                         gacgaAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCaa                                                                         <  2:1635028/100‑1 (MQ=255)
                           cgaAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAAcc                                                                       <  2:1920907/100‑1 (MQ=255)
                            taggtaCGATCTGCAAAGCCTGGAGTGTACATGTTCAGTCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAAc                                                                        <  1:3922032/93‑1 (MQ=255)
                            gaagaaCGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAAccc                                                                      <  1:1691398/100‑1 (MQ=255)
                                 aCGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTAttt                                                                 >  2:3279781/1‑100 (MQ=255)
                                    aTGTGCATAGCCTGAAGTGTATATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCg                                                              >  2:741638/1‑100 (MQ=255)
                                     tgtgCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCgg                                                             <  2:2176851/100‑1 (MQ=255)
                                       tgCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGaa                                                           <  1:3279795/100‑1 (MQ=255)
                                        gCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGaaa                                                          >  1:3364597/1‑100 (MQ=255)
                                        gCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGaaa                                                          >  1:1569262/1‑100 (MQ=255)
                                          aTAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAaga                                                        <  1:657893/100‑1 (MQ=255)
                                           tAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAagag                                                       >  1:625948/1‑100 (MQ=255)
                                            aGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGc                                                      >  1:3009008/1‑100 (MQ=255)
                                             gCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATATTCGGAAAGAGcc                                                     >  1:2436047/1‑100 (MQ=255)
                                                tGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCg                                                  >  1:2930610/1‑100 (MQ=255)
                                                tGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCg                                                  >  1:2272183/1‑100 (MQ=255)
                                                tGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCg                                                  >  2:4150319/1‑100 (MQ=255)
                                                   aGTGTACATAATCAATCCGGCCCCTGTGAGTTACTGTATGGATGTACAGTACATGCAGTGACAACACAGCTCAACACTATTTCCGGTAAGAGCCTCGCaa                                               >  2:3901973/1‑100 (MQ=255)
                                                      gtACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAAtt                                            >  2:3090087/1‑100 (MQ=255)
                                                      gtACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAAtt                                            >  1:1470185/1‑100 (MQ=255)
                                                       tACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAATATAAACCCTATTTTCGGAAAGAGCCTCGCAAAttt                                           >  2:970242/1‑100 (MQ=255)
                                                        aCATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCGAAtttt                                          >  2:2731366/1‑100 (MQ=255)
                                                        aCATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCGAAtttt                                          <  1:2731373/100‑1 (MQ=255)
                                                         cATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTg                                         <  2:2925770/100‑1 (MQ=255)
                                                         cATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTg                                         >  2:3248202/1‑100 (MQ=255)
                                                         cATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTg                                         <  2:1985588/100‑1 (MQ=255)
                                                          aTAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGt                                        >  1:653187/1‑100 (MQ=255)
                                                          aTAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGt                                        >  1:257517/1‑100 (MQ=255)
                                                          aTAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGt                                        <  1:1367198/100‑1 (MQ=255)
                                                           tAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTc                                       >  2:99028/1‑100 (MQ=255)
                                                           tAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAAACCTATTTTCGGAAAGAGCCTCGCAAATTTTGTc                                       >  2:901086/1‑100 (MQ=255)
                                                            aatcaatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCg                                      <  1:1128913/100‑1 (MQ=255)
                                                            aatcaatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCg                                      <  2:3657735/100‑1 (MQ=255)
                                                                aatcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTgg                                  <  2:20884/100‑1 (MQ=255)
                                                                 atcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGt                                 >  2:4035295/1‑100 (MQ=255)
                                                                 atcCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGt                                 >  1:1810339/1‑100 (MQ=255)
                                                                      gCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACgg                            >  2:1388283/1‑100 (MQ=255)
                                                                       ccccTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACggg                           <  1:434271/100‑1 (MQ=255)
                                                                          cTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGaaa                        >  1:520360/1‑100 (MQ=255)
                                                                           tgtgAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGaaaa                       <  1:3598578/100‑1 (MQ=255)
                                                                            gtgAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAAc                      <  2:3012034/100‑1 (MQ=255)
                                                                            gtgAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAAc                      <  2:3054540/100‑1 (MQ=255)
                                                                             tgAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACa                     <  2:244133/100‑1 (MQ=255)
                                                                              gAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACAt                    >  2:489772/1‑100 (MQ=255)
                                                                               aGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATa                   >  1:3041906/1‑100 (MQ=255)
                                                                               aGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATa                   >  2:1232924/1‑100 (MQ=255)
                                                                               aGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATa                   <  1:3651428/100‑1 (MQ=255)
                                                                               aGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATa                   <  1:1730460/100‑1 (MQ=255)
                                                                                 ttACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATaaa                 <  1:1002507/100‑1 (MQ=255)
                                                                                 ttACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATaaa                 <  2:336058/100‑1 (MQ=255)
                                                                                  tACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAAt                <  1:3889691/100‑1 (MQ=255)
                                                                                   aCTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAAtt               >  1:3662098/1‑100 (MQ=255)
                                                                                   aCTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAAt                <  2:949288/99‑1 (MQ=255)
                                                                                    gtcactGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAAtt               >  2:2714807/6‑99 (MQ=255)
                                                                                     tcactGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTaa             >  1:2365610/5‑100 (MQ=255)
                                                                                     tcactGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTaa             >  1:2650682/5‑100 (MQ=255)
                                                                                      gTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTa              >  2:492461/1‑98 (MQ=255)
                                                                                      gTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAAt            >  1:1388679/1‑100 (MQ=255)
                                                                                      gTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAAt            >  2:725014/1‑100 (MQ=255)
                                                                                      cactGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAAt            >  1:2666037/4‑100 (MQ=255)
                                                                                      cactGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAAt            >  1:3350537/4‑100 (MQ=255)
                                                                                       tATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATc           <  2:1569259/100‑1 (MQ=255)
                                                                                       tATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATc           <  2:312374/100‑1 (MQ=255)
                                                                                       actGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATc           >  2:2612616/3‑100 (MQ=255)
                                                                                       actGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATc           >  1:2299032/3‑100 (MQ=255)
                                                                                       actGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATc           >  1:3363484/3‑100 (MQ=255)
                                                                                        ctGGATGTACTGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTATCGTTGGTGACGGGAAAACATAAATTAATCt          >  2:351104/2‑100 (MQ=255)
                                                                                        atcgATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCAGAAAGAGCCTCGCAAATTTCGCCATTAAAAACGGGAACACATAAATTaacat          <  2:779246/97‑4 (MQ=255)
                                                                                           gATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGc       <  1:887784/100‑1 (MQ=255)
                                                                                           gATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGc       >  1:2170195/1‑100 (MQ=255)
                                                                                           gATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGc       >  1:2511365/1‑100 (MQ=255)
                                                                                            aTGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGcc      <  2:1130702/100‑1 (MQ=255)
                                                                                             tGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGccc     <  1:2702266/100‑1 (MQ=255)
                                                                                                acagtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCtt  >  2:1712638/1‑100 (MQ=255)
                                                                                                acagtacaTCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCtt  <  1:944391/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
AATTTTACTTGCTGCGGATGAGAACGACGAAGAACGATGTGCATAGCCTGAAGTGTACATAATCAATCCAGCCCCTGTGAGTTACTGTATGGATGTACAGTACATCCAGTGACAACAAAGATCAACCCTATTTTCGGAAAGAGCCTCGCAAATTTTGTCGTTGGTGACGGGAAAACATAAATTAATCTTGCCCCTT  >  NC_000913/1020860‑1021055

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