Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,057,150 C→T 21.6% intergenic (‑201/+28) ibsC ← / ← serA toxic membrane protein/D‑3‑phosphoglycerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,057,1500CT21.6% 403.7 / 103.6 203intergenic (‑201/+28)ibsC/serAtoxic membrane protein/D‑3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (80/78);  new base T (28/16);  total (109/94)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.28e-01

Only 100 of 213 total aligned reads displayed.

ATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTC  >  NC_000913/3057051‑3057249
                                                                                                   |                                                                                                   
aTCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAcc                                                                                                     <  2:3546710/100‑1 (MQ=255)
aTCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAcc                                                                                                     <  1:3468148/100‑1 (MQ=255)
 tCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCt                                                                                                    >  2:2747075/1‑100 (MQ=255)
      cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTCGAATGCGCAACGCGGCAACGGAGTCGAGCAGGGATAAAAAAACGGGCAAGTCAGTGACCTGc                                                                                                  >  1:2530783/1‑97 (MQ=255)
      cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCg                                                                                               <  1:1863307/100‑1 (MQ=255)
      cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCg                                                                                               <  2:645820/100‑1 (MQ=255)
      cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCg                                                                                               >  1:1178588/1‑100 (MQ=255)
      cTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCg                                                                                               >  2:1163004/1‑100 (MQ=255)
      cTTTGATTTATGTTCTGGAGCAGAATGGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGGTAAAAAATCGGGCAAGTCAGTGGGCTGCCCg                                                                                               <  1:4086412/100‑1 (MQ=255)
        ttGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGtt                                                                                             <  2:1844635/100‑1 (MQ=255)
        ttGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGtt                                                                                             <  1:2603650/100‑1 (MQ=255)
         tGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTg                                                                                            <  1:1340318/100‑1 (MQ=255)
         tGATTTATTTTCGGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACCGTGTGGAGAATGGATAAAAAAACGGGCAGGTCTCTGACTTGCCCGttt                                                                                            <  1:3751374/100‑2 (MQ=255)
          gATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGtttt                                                                                           <  1:3422171/100‑3 (MQ=38)
          gATTTATTTTCTGGAGCAGACTCGAAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttt                                                                                            <  1:1518694/99‑2 (MQ=38)
          gATGTATTTTCTGGAGCAGACTCGAAAAGAAGAAGTAACAACGCGCCAACAGTGTGGAGAAGGGATTAAAAAACTGGCATGTCACTGACTTGCCCGtttt                                                                                           <  1:2910108/100‑3 (MQ=255)
           atttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGACGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGa                                                                                           <  1:797020/99‑1 (MQ=255)
           atttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttttt                                                                                          <  2:337804/100‑4 (MQ=38)
           atatattgTCTGGAGCAGACTCCAAAAGTAGAAAACGCAACGCGGCAACGGTGTAAAAAAAAAAAAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTg                                                                                            <  1:4123620/95‑1 (MQ=255)
             ttatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAGGTCACTGACTTGCCCGttttttt                                                                                        <  2:361198/100‑6 (MQ=37)
             ttattgTCTGGAGTAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGCGAAGGGATAAAAAGACGGGCAAGTCAGTGACCTGCCCGTTgctct                                                                                        >  1:2727847/1‑96 (MQ=255)
              tatttTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGAtttt                                                                                       >  2:2483928/1‑100 (MQ=255)
                ttttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCa                                                                                     <  1:2123580/100‑1 (MQ=255)
                 tttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGGGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                                    >  2:1498014/1‑100 (MQ=255)
                 tttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAACCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                                    <  1:4362365/100‑1 (MQ=255)
                 tttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                                    >  1:3199250/1‑100 (MQ=255)
                  ttCTGGATCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGCGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCaga                                                                                   >  1:1643459/1‑100 (MQ=255)
                  ttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCaga                                                                                   <  2:737398/100‑1 (MQ=255)
                  ttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCaga                                                                                   <  1:1619814/100‑1 (MQ=255)
                  ttCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCaga                                                                                   <  2:2028489/100‑1 (MQ=255)
                   tCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCag                                                                                    <  2:2239701/98‑1 (MQ=255)
                    cTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAgaga                                                                                 <  2:1706873/100‑1 (MQ=255)
                    cTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGAGATAAAAGAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAgaga                                                                                 >  1:1311273/1‑100 (MQ=255)
                      ggAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAg                                                                               <  1:2290462/100‑1 (MQ=255)
                      ggAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAg                                                                               >  2:1003648/1‑100 (MQ=255)
                         gCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAgggg                                                                            >  1:2217944/1‑100 (MQ=255)
                            gACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAt                                                                         >  2:3696337/1‑100 (MQ=255)
                             aCTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAAtt                                                                        <  2:2167006/100‑1 (MQ=255)
                                cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                                     >  2:2691738/1‑100 (MQ=255)
                                cGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGAACTGCCCGTTGATTTTCAGAGAAGGGGAATTAGt                                                                     >  1:298921/1‑100 (MQ=255)
                                 gCAAAGTAGAATACGCAACGCGGCAACGGTGAGGAGAGGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTATCAGAGAAGAGGAATTAGTa                                                                    >  1:1656001/1‑100 (MQ=255)
                                   aaaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAca                                                                  <  1:1720358/100‑1 (MQ=255)
                                    aaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcag                                                                 >  2:3253560/1‑100 (MQ=255)
                                    aaGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcag                                                                 <  2:3753082/100‑1 (MQ=255)
                                    aaGTAGAATACGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTAcag                                                                 <  1:1058146/100‑1 (MQ=255)
                                         gAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGAc                                                            <  2:1946873/100‑1 (MQ=255)
                                         gAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGAc                                                            <  1:3259433/100‑1 (MQ=255)
                                         gAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGAc                                                            >  1:1000867/1‑100 (MQ=255)
                                             gcgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgc                                                        <  2:4397313/100‑1 (MQ=255)
                                              cgcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgcg                                                       >  2:3961812/1‑100 (MQ=255)
                                              cgcAACGCGGCAACGGCGTGGAGAAGGGCTAAAAAAACGGCCAAGTCAGTGACCTGCCCGGTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgcg                                                       >  1:2421011/1‑100 (MQ=255)
                                               gcAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGgcgc                                                      >  1:245538/1‑100 (MQ=255)
                                                cAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGcgcg                                                     <  2:3593225/100‑1 (MQ=255)
                                                cAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGcgcg                                                     >  2:3102467/1‑100 (MQ=255)
                                                  aCGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGaa                                                   >  1:1939398/1‑100 (MQ=255)
                                                     cgGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGCTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATgg                                                >  2:1408112/1‑100 (MQ=255)
                                                       gCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTa                                              >  1:3042458/1‑100 (MQ=255)
                                                        cAACGGCGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAc                                             >  1:2697912/1‑100 (MQ=255)
                                                         aaCGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAcc                                            <  1:4206569/100‑1 (MQ=255)
                                                                tgGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATa                                     <  1:602839/100‑1 (MQ=255)
                                                                 ggAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAg                                    <  2:1326739/100‑1 (MQ=255)
                                                                 ggAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAg                                    <  1:3980659/100‑1 (MQ=255)
                                                                 ggAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAg                                    >  2:795521/1‑100 (MQ=255)
                                                                   agaAGGGATAAAAAAACGGGCAAGTCGGTGACCTGCCCGGTGATTCTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAAt                                      <  1:2806315/96‑1 (MQ=255)
                                                                      aGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGACATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTc                               >  1:2378047/1‑100 (MQ=255)
                                                                       tctgaaaatcaaCGGGCAGGTCACTGACTTGCCCGTTGACTTTCAGAGAAGGGCAACTAGTACATCACACGGGCGCGAATGGTAACCGGAAGTGCTGTCa                              >  1:742223/11‑100 (MQ=12)
                                                                         tgaaaatcaaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCAtt                            >  2:3118894/9‑100 (MQ=37)
                                                                         tgaaaatcaaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCAtt                            >  2:96606/9‑100 (MQ=37)
                                                                           tAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGc                          >  1:81391/1‑100 (MQ=255)
                                                                            aaaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGcc                         <  2:2500838/100‑1 (MQ=255)
                                                                             aaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCt                        >  1:1130950/1‑100 (MQ=255)
                                                                             aaaaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCt                        >  1:3295713/1‑100 (MQ=255)
                                                                              atcaaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTg                       >  2:2377322/4‑100 (MQ=38)
                                                                               aaaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGc                      >  1:46245/1‑100 (MQ=255)
                                                                                caaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCCGAGCCGGGGACTTAGTCCAGCAGACGGGCGCGAAAGGCACCCGGAATAGCTCTCCTTGcc                         >  1:68930/2‑96 (MQ=19)
                                                                                caaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCa                     >  1:3791075/2‑100 (MQ=39)
                                                                                aaaCGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGTCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCa                     <  2:1751793/100‑1 (MQ=255)
                                                                                 aaCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAg                    >  2:890297/1‑100 (MQ=39)
                                                                                  aCGGGCAGGTCACTGTCTTGCCCGTTGATTTTCAGAGACGGGGAATTAGTACAGCAGACGGGCGCGGATGGTACGCGGAATAGCTTTCATTGcc                         >  1:1226022/1‑94 (MQ=255)
                                                                                  aCGGGCAGGTCACTGACTTGCCTGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAGCCGGAATAGCTTTCATTGCCTGGAg                    >  1:3075076/1‑99 (MQ=255)
                                                                                  aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc                   >  2:3136832/1‑100 (MQ=255)
                                                                                  aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc                   >  2:1793464/1‑100 (MQ=39)
                                                                                  aCGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgc                   >  1:4146991/1‑100 (MQ=39)
                                                                                   cGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcg                  >  2:1225284/1‑100 (MQ=255)
                                                                                   cGGGCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAgcg                  >  2:1721899/1‑100 (MQ=255)
                                                                                    gggCAGGTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTAGCCGGAATAGCTTTCATTGCCTGCAgcg                  >  1:3990092/1‑99 (MQ=255)
                                                                                     ggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCt                >  2:1596750/1‑100 (MQ=255)
                                                                                     ggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCt                >  2:1079595/1‑100 (MQ=255)
                                                                                     ggCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGAAATAGCTTTCATTGCCTGCAgcgc                 >  1:3283562/1‑99 (MQ=255)
                                                                                       cgaGTCAGTGCCCTCCGCGTTGATGTTTAGAGAGAGGGAATTAGTACAGCAGACGGGCGGGAATGGTACCCGGATTAGCTGTCATTGCCTGCAGCGCttt              <  1:2048135/98‑1 (MQ=255)
                                                                                        aggTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtttt             <  2:586203/98‑1 (MQ=255)
                                                                                        aaGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCtttt             >  2:2902686/1‑100 (MQ=255)
                                                                                         ggTCACTGACTTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttttt            >  2:3179573/2‑100 (MQ=255)
                                                                                          ggcAGTGACCGGCCCGTTGATTTTCAGGGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCttttt            <  1:1808949/97‑1 (MQ=255)
                                                                                           tCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCg          <  2:1686064/100‑1 (MQ=255)
                                                                                            cAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCgg         <  1:964837/100‑1 (MQ=255)
                                                                                            cAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCgg         <  2:3199267/100‑1 (MQ=255)
                                                                                                gACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAAc     >  2:3861943/1‑100 (MQ=255)
                                                                                                gACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAAc     >  2:2543272/1‑100 (MQ=255)
                                                                                                  ccTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGt   >  2:3265474/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ATCTGACTTTGATTTATTTTCTGGAGCAGACTCGCAAAGTAGAATGCGCAACGCGGCAACGGTGTGGAGAAGGGATAAAAAAACGGGCAAGTCAGTGACCTGCCCGTTGATTTTCAGAGAAGGGGAATTAGTACAGCAGACGGGCGCGAATGGTACCCGGAATAGCTTTCATTGCCTGCAGCGCTTTTTCGGCAACGTC  >  NC_000913/3057051‑3057249

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