Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 996,684 A→G S18P (TCC→CCC)  ycbP ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913996,6840AG99.4% 553.7 / NA 175S18P (TCC→CCC) ycbPhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (60/113);  total (60/115)
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 176 total aligned reads displayed.

GTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGAGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCTTTTTTTGTT  >  NC_000913/996594‑996770
                                                                                          |                                                                                      
gTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGc                                                                                      >  1:267431/1‑93 (MQ=255)
 tAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCt                                                                                     <  2:6087453/93‑1 (MQ=255)
  aGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTg                                                                                    <  1:3013094/93‑1 (MQ=255)
    aTCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGa                                                                                  >  1:5716415/1‑93 (MQ=255)
      cTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGc                                                                                <  1:3487642/93‑1 (MQ=255)
      cTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGc                                                                                >  2:2399914/1‑93 (MQ=255)
        tCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCga                                                                              <  2:628329/93‑1 (MQ=255)
          cGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCgaga                                                                            >  2:3629421/1‑93 (MQ=255)
           ggggCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGaa                                                                           <  1:797208/93‑1 (MQ=255)
            gggCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAg                                                                          <  1:1567635/93‑1 (MQ=255)
            gggCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGCGCTGGAGCGAGAAg                                                                          <  1:4067828/93‑1 (MQ=255)
             ggCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAgg                                                                         <  1:1555304/93‑1 (MQ=255)
             ggCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGCCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAgg                                                                         >  2:3567824/1‑93 (MQ=255)
              gCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGa                                                                        <  1:2294067/93‑1 (MQ=255)
                gAAGTTTTTCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGGAGGaaa                                                                      >  2:3530723/1‑93 (MQ=255)
                gAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGaaa                                                                      >  2:1347380/1‑93 (MQ=255)
                gAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGaaa                                                                      <  2:4237410/93‑1 (MQ=255)
                   gTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCg                                                                   >  1:526328/1‑93 (MQ=255)
                   gTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCg                                                                   >  1:4927662/1‑93 (MQ=255)
                    ttttGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGa                                                                  <  2:1703355/93‑1 (MQ=255)
                      ttGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAgg                                                                <  2:1880647/93‑1 (MQ=255)
                      ttGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAgg                                                                <  2:2501273/93‑1 (MQ=255)
                      ttGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAgg                                                                >  1:1957058/1‑93 (MQ=255)
                       tGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGa                                                               <  1:4605450/93‑1 (MQ=255)
                        gCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGAc                                                              <  2:896252/93‑1 (MQ=255)
                         cAGATTCGAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACt                                                             <  2:754879/93‑1 (MQ=255)
                         cAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACt                                                             <  2:5492104/93‑1 (MQ=255)
                         cAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACt                                                             <  1:3642315/93‑1 (MQ=255)
                         cAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACt                                                             <  2:1009487/93‑1 (MQ=255)
                          aGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTa                                                            >  1:92616/1‑93 (MQ=255)
                          aGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTa                                                            <  2:3354611/93‑1 (MQ=255)
                           gATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTAc                                                           <  1:1955867/93‑1 (MQ=255)
                            aTTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTAcc                                                          <  2:4287369/93‑1 (MQ=255)
                            aTTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTAcc                                                          <  1:5752930/93‑1 (MQ=255)
                             ttCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTAccc                                                         <  1:6061217/93‑1 (MQ=255)
                               ccAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGc                                                       <  1:1635370/93‑1 (MQ=255)
                                cAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGcc                                                      <  2:1260007/93‑1 (MQ=255)
                                  gTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAg                                                    >  2:6326738/1‑93 (MQ=255)
                                  gTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAg                                                    <  2:2064583/93‑1 (MQ=255)
                                    gATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAggg                                                  <  2:5884334/93‑1 (MQ=255)
                                    gATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAggg                                                  >  2:4968586/1‑93 (MQ=255)
                                      tAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAggggg                                                >  2:2957972/1‑91 (MQ=255)
                                      tAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGtg                                                >  1:6142399/1‑93 (MQ=255)
                                           ctctACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGac                                           >  1:2478418/1‑93 (MQ=255)
                                           ctctACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGac                                           >  2:1116783/1‑93 (MQ=255)
                                            tctACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGaca                                          <  2:4741322/93‑1 (MQ=255)
                                            tctACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGaca                                          <  1:6379163/93‑1 (MQ=255)
                                             ctACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGacac                                         >  2:2812895/1‑93 (MQ=255)
                                              tACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACg                                        >  1:4613087/1‑93 (MQ=255)
                                               aCATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGc                                       >  1:3035101/1‑93 (MQ=255)
                                                cATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCa                                      <  1:1767595/93‑1 (MQ=255)
                                                cATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCa                                      <  1:2325014/93‑1 (MQ=255)
                                                 aTCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCAt                                     <  2:675513/93‑1 (MQ=255)
                                                  tCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATa                                    >  2:3957095/1‑93 (MQ=255)
                                                   ccAGTCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATAc                                   <  2:15758/93‑1 (MQ=255)
                                                   ccAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATAc                                   <  1:5706902/93‑1 (MQ=255)
                                                    cAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATAct                                  <  1:5162417/93‑1 (MQ=255)
                                                    cAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGCTTGGAGCGAGAAGGAAAGCGAGGACTCCCCGCCAGGGTGATGACACGCATAct                                  <  1:598918/93‑1 (MQ=255)
                                                      ggCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATActct                                >  1:4351757/1‑93 (MQ=255)
                                                       gCCATTTATTTTTTCCCGCGCATATTCCAGCAAGAGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATActctc                               <  1:1426247/93‑1 (MQ=255)
                                                       gCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATActctc                               <  2:6659235/93‑1 (MQ=255)
                                                        ccATTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATActctct                              <  2:2619522/93‑1 (MQ=255)
                                                          aTTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTcc                            <  1:3828621/93‑1 (MQ=255)
                                                          aTTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTcc                            <  1:1334568/93‑1 (MQ=255)
                                                          aTTTAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTcc                            <  2:3894344/93‑1 (MQ=255)
                                                           tttAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTCCCCGCCAGGGTGATGACACGCATACTCTCTCCt                           <  1:1066041/93‑1 (MQ=255)
                                                           tttAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCt                           <  2:3419513/93‑1 (MQ=255)
                                                            ttAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCtt                          <  1:2283712/93‑1 (MQ=255)
                                                             tAGTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTta                         >  2:3677714/1‑93 (MQ=255)
                                                               gTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTtata                       <  1:5883093/93‑1 (MQ=255)
                                                               gTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTtata                       >  2:4179176/1‑93 (MQ=255)
                                                               gTTTTTCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTtata                       >  2:5299481/1‑93 (MQ=255)
                                                                tttttCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCCAGGACTACCCGCCAGGGTGATGACTCGCATACTCTCTGCTTATaa                      >  1:3283964/1‑93 (MQ=255)
                                                                  tttCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAAcc                    <  1:287141/93‑1 (MQ=255)
                                                                  tttCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAACGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAAcc                    <  1:118459/93‑1 (MQ=255)
                                                                    tCCCGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCaa                  <  1:4003202/93‑1 (MQ=255)
                                                                     cccGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAAt                 <  1:5735389/93‑1 (MQ=255)
                                                                      ccGCGCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAAtt                <  1:2171070/93‑1 (MQ=255)
                                                                       cgcgCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTg               >  1:5000237/1‑93 (MQ=255)
                                                                       cgcgCATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTg               <  1:3633664/93‑1 (MQ=255)
                                                                        gcgcATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGt              <  2:4900705/93‑1 (MQ=255)
                                                                        gcgcATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGt              <  1:6033550/93‑1 (MQ=255)
                                                                         cgcATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGtt             <  2:4558135/93‑1 (MQ=255)
                                                                         cgcATATTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGtt             <  1:4724910/93‑1 (MQ=255)
                                                                         cgcATATTCCAGCAAGGAGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGtt             <  2:3964386/93‑1 (MQ=255)
                                                                            atatTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCtt          <  2:4333895/93‑1 (MQ=255)
                                                                            atatTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCtt          <  2:5883669/93‑1 (MQ=255)
                                                                             tatTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCttt         >  1:2066957/1‑93 (MQ=255)
                                                                             tatTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCttt         <  2:6548443/93‑1 (MQ=255)
                                                                             tatTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCttt         <  1:2318444/93‑1 (MQ=255)
                                                                             tatTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCttt         <  1:2469730/93‑1 (MQ=255)
                                                                             tatTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCttt         >  1:3126558/1‑93 (MQ=255)
                                                                              atTCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCtttt        >  1:3439033/1‑93 (MQ=255)
                                                                                tCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCtttttt      <  2:5895396/93‑1 (MQ=255)
                                                                                tCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCtttttt      <  2:6718401/93‑1 (MQ=255)
                                                                                tCCAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCtttttt      <  2:698017/93‑1 (MQ=255)
                                                                                 ccAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCttttttt     <  2:2657379/93‑1 (MQ=255)
                                                                                 ccAGAAAGGGGCTGGAGCGAGAAGAAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCttttttt     <  1:3530723/93‑1 (MQ=255)
                                                                                  cAGCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCTTTTTTTg    <  2:2168932/93‑1 (MQ=255)
                                                                                    gCAAGGGGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCTTTTTTTGtt  >  2:3124929/1‑93 (MQ=255)
                                                                                          |                                                                                      
GTAGATCTTCCGGGGCGAAGTTTTGCAGATTCCAGTGATAAACCTCTACATCCAGGCCATTTAGTTTTTCCCGCGCATATTCCAGCAAGGAGCTGGAGCGAGAAGGAAAGCGAGGACTACCCGCCAGGGTGATGACACGCATACTCTCTCCTTATAACCAATTGTTCTTTTTTTGTT  >  NC_000913/996594‑996770

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