Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,950,468 G→T 100% intergenic (+48/‑39) ilvL → / → ilvX ilvG operon leader peptide/uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,950,4680GT98.9% 556.8 / ‑3.0 171intergenic (+48/‑39)ilvL/ilvXilvG operon leader peptide/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (80/89);  total (81/90)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.52e-01

Only 100 of 178 total aligned reads displayed.

GATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑GTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAA  >  NC_000913/3950372‑3950564
                                                                                                 |                                                                                                
gATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑tttt                                                                                               <  2:1544639/100‑1 (MQ=255)
   tattatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑Ttttttt                                                                                            >  2:1655190/1‑100 (MQ=255)
    attatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑tttttttt                                                                                           <  1:3739275/100‑1 (MQ=255)
     ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑Ttttttttt                                                                                          <  1:2480022/100‑1 (MQ=255)
       atCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACTAACTAAGACCGCCGCACCGAAAGGTCCGTGG‑TTTTTTTTTGa                                                                                        <  1:2788806/100‑1 (MQ=255)
       atCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGa                                                                                        >  1:2152769/1‑100 (MQ=255)
         cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGTcc                                                                                      >  1:3808361/1‑100 (MQ=255)
         cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGAcc                                                                                      <  1:1220343/100‑1 (MQ=255)
         cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGAcc                                                                                      >  1:3940355/1‑100 (MQ=255)
         cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGAcc                                                                                      <  2:4114064/100‑1 (MQ=255)
          ccaccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCgggg‑ttttttttttgacc                                                                                     <  1:3295350‑M2/95‑6 (MQ=255)
          ccaccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCt                                                                                     >  2:1507273‑M2/6‑100 (MQ=255)
          ccaccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCt                                                                                     >  1:1069917‑M2/6‑100 (MQ=255)
          ccaccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCt                                                                                     >  1:1930272‑M2/6‑100 (MQ=255)
            accGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTa                                                                                   <  1:3489148‑M2/97‑1 (MQ=255)
            accGTGCGGGGCTGCACTTCGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTa                                                                                   <  2:4037719‑M2/97‑1 (MQ=255)
             ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGT‑TTTTTTTTTGACCTTaa                                                                                  >  1:328240‑M2/3‑100 (MQ=255)
             ccGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTaa                                                                                  >  1:3173947‑M2/3‑100 (MQ=255)
              cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTaaa                                                                                 <  1:2953207‑M2/99‑1 (MQ=255)
               gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTTACCtt                                                                                    >  2:3413460/1‑96 (MQ=255)
               gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGCACTTa                                                                                   >  2:3820296/1‑97 (MQ=255)
                 gCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTaaaaa                                                                               <  2:3759482/99‑1 (MQ=255)
                     ggCTGCACTTGGACGAGGAAAGGCTTAGAGATCCAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACa                                                                             <  2:2493998/97‑1 (MQ=255)
                      gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAAc                                                                         <  2:1318671/100‑1 (MQ=255)
                      gCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAAc                                                                         >  1:159892/1‑100 (MQ=255)
                       cTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTTTCCTTAAAAACaata                                                                          >  2:3874325/1‑95 (MQ=255)
                        tGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCg                                                                       <  1:2082923/100‑1 (MQ=255)
                           aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCgagg                                                                    <  2:3134538/100‑1 (MQ=255)
                           aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCgagg                                                                    >  1:679378/1‑100 (MQ=255)
                            cTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTGACC‑TAAAAAAATAAACgaag                                                                    >  1:232802/1‑97 (MQ=255)
                              tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                                 <  2:1959182/100‑1 (MQ=255)
                              tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                                 >  2:1100229/1‑100 (MQ=255)
                               ggACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                                 <  2:1232208/99‑1 (MQ=255)
                               ggACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                                >  1:840476/1‑100 (MQ=255)
                               ggACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                                >  1:115656/1‑100 (MQ=255)
                                gACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAg                                                               <  1:698777/100‑1 (MQ=255)
                                 aCGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGa                                                              <  1:2001123/100‑1 (MQ=255)
                                  cGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTTTCCTTTTAATCATATCCGAGGAGCAGAc                                                             >  2:1133449/1‑100 (MQ=255)
                                  cGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACAGAACCGAGGAGCAGAc                                                             >  1:4067245/1‑100 (MQ=255)
                                     ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                          <  2:3193380/100‑1 (MQ=255)
                                     ggAAAGGCTTAGAGATCAAGCCTGAACAAACAAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                          <  1:6947/100‑1 (MQ=255)
                                      gAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                         >  1:1493287/1‑100 (MQ=255)
                                      gAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                         <  1:3864532/100‑1 (MQ=255)
                                       aaaGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGa                                                        >  1:3124212/1‑100 (MQ=255)
                                          ggCTTAGAGATTAAGCCTTAACGAACTAAGACCCCGGGACCAAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATa                                                     <  2:1152163/100‑1 (MQ=255)
                                            cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                   >  1:1041007/1‑100 (MQ=255)
                                            cTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAAc                                                   >  2:357177/1‑100 (MQ=255)
                                             ttAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACa                                                  <  1:3656137/100‑1 (MQ=255)
                                             ttAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACa                                                  >  1:3378131/1‑100 (MQ=255)
                                              tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                                 <  1:3188556/100‑1 (MQ=255)
                                               agagATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAATCGAGGAGCAGACAATGAATAACAGc                                                >  2:683892/1‑100 (MQ=255)
                                                 agaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGcac                                              <  2:4067089/100‑1 (MQ=255)
                                                    ttaaGCCGTAACGAACTACGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAAt                                                      <  2:2191148/88‑1 (MQ=255)
                                                    tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGTTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaa                                           >  1:1854631/1‑100 (MQ=255)
                                                    tCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACaaa                                           >  1:2533486/1‑100 (MQ=255)
                                                        gcgttAACGAACTAAGACCCCCGCACCGAAAGGCCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTc                                       <  2:1859424/97‑1 (MQ=255)
                                                        gCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTc                                       >  2:1009247/1‑100 (MQ=255)
                                                         ccTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTc                                       >  1:1056249/1‑99 (MQ=255)
                                                         ccTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCt                                      <  1:3039613/100‑1 (MQ=255)
                                                         ccTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCt                                      >  1:2632817/1‑100 (MQ=255)
                                                           ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                    <  2:1166857/100‑1 (MQ=255)
                                                             aaCGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGttt                                  <  1:1664059/100‑1 (MQ=255)
                                                              aCGAACTAAGACCCCCGCACCGAGAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                                 <  1:3531747/100‑1 (MQ=255)
                                                              aCGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                                 <  2:2936790/100‑1 (MQ=255)
                                                              aCGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                                 >  1:458744/1‑100 (MQ=255)
                                                               cgtaCTAAGACCCCCGCACCGAAATTTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtct                                <  2:3656730/97‑1 (MQ=255)
                                                                gAACTAAGAACCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                                 <  2:519046/98‑1 (MQ=255)
                                                                 aaCTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                              >  2:3434619/1‑100 (MQ=255)
                                                                 aaCTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                              <  2:3163555/100‑1 (MQ=255)
                                                                 aaCTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                              >  1:3695941/1‑100 (MQ=255)
                                                                 aaCGAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                              <  2:2476824/100‑1 (MQ=255)
                                                                  aCTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAACATTCTGTTTCTCa                              >  1:496834/1‑99 (MQ=255)
                                                                   cTAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAg                            <  2:72376/100‑1 (MQ=255)
                                                                    tAAGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGa                           <  1:2766033/100‑1 (MQ=255)
                                                                    tAAGACCCCCGCACCGAAAGGGCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                             <  1:1103692/98‑1 (MQ=255)
                                                                     aaGACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAt                          >  2:1925542/1‑100 (MQ=255)
                                                                     aaGACCCCCGCACCGAAAGGGCC‑GAG‑GCTTTTTTTGACCTTAAAAACATAACCAAGGAGCCGACACAGAATCGTAGCTCAAAAtccagacg                                  >  2:3903442/1‑84 (MQ=255)
                                                                       gACCCCCGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTc                        >  1:1718742/1‑100 (MQ=255)
                                                                         cccccGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAg                      >  1:2581673/1‑100 (MQ=255)
                                                                         cccccGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAg                      >  1:785058/1‑100 (MQ=255)
                                                                          ccccGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                     <  2:2034695/100‑1 (MQ=255)
                                                                           cccGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGa                    >  1:1263459/1‑100 (MQ=255)
                                                                           cccGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGa                    <  1:1833243/100‑1 (MQ=255)
                                                                            ccGCACCGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGAc                   >  1:1752184/1‑100 (MQ=255)
                                                                                aCCGAAAGGTCCGGGGTTTTTTTTTTTTTTTTTTAAACATAAACGAGGAGCAGACAATCAATTATTGCTTAAAATTCTGTTTTTGAAGTTTCAGGa                    >  2:2966715/1‑96 (MQ=255)
                                                                                aCCGAAAGGTCC‑GGG‑GTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg               <  1:1545441/99‑1 (MQ=255)
                                                                                 ccGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg               >  1:3729674/1‑99 (MQ=255)
                                                                                 ccGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGa              <  1:2948800/100‑1 (MQ=255)
                                                                                 ccGAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGa              <  2:2434760/100‑1 (MQ=255)
                                                                                   gAAAGGTCCGGGT‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac            >  1:1045392/1‑100 (MQ=255)
                                                                                   gAAAGGTCCGGGG‑TTTTTTTTTTTCTTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGcac            >  1:503075/1‑100 (MQ=255)
                                                                                   gAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac            >  1:3005219/1‑100 (MQ=255)
                                                                                   gAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac            >  1:549848/1‑100 (MQ=255)
                                                                                   gAAAGGTCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaac            <  2:2000035/100‑1 (MQ=255)
                                                                                         tCCGGGG‑TTTTTTTTTTTTTTTTTAAACATAACCGAGGAGCACACAATTCATCACAGCCCAGAATTCTGTTTCTTAAGATTCAGGAcccgccacta          >  2:1831272/1‑87 (MQ=255)
                                                                                         tCCGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAactaacta      >  1:1776240/1‑100 (MQ=255)
                                                                                          ccGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAATATTCAGGACGGGGAactaacta      <  2:2395151/99‑1 (MQ=255)
                                                                                          ccGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTAt     <  2:3464863/100‑1 (MQ=255)
                                                                                          ccGGGG‑TTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTAt     >  2:3209345/1‑100 (MQ=255)
                                                                                             ggg‑tTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGaa  <  1:2568143/100‑1 (MQ=255)
                                                                                                 |                                                                                                
GATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGG‑GTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAA  >  NC_000913/3950372‑3950564

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