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,4680GT97.9% 280.0 / ‑1.2 95intergenic (+48/‑39)ilvL/ilvXilvG operon leader peptide/uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base T (48/44);  total (50/45)
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.59e-01

Only 100 of 101 total aligned reads displayed.

GATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAA  >  NC_000913/3950372‑3950564
                                                                                                |                                                                                                
gATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTATAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGtttt                                                                                               >  2:1232222/1‑100 (MQ=255)
 attattATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttt                                                                                              <  1:1856484/100‑1 (MQ=255)
 attattATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttt                                                                                              <  1:1435830/100‑1 (MQ=255)
  taattatCCCACCATGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGtttttt                                                                                             >  1:1740704/3‑100 (MQ=255)
   tattatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttttt                                                                                            >  1:1448199/1‑100 (MQ=255)
    attatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGtttttttt                                                                                           <  2:2083422/100‑1 (MQ=255)
     ttatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGGTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTtttttttt                                                                                          >  2:277535/1‑100 (MQ=255)
      tatCCCACCGTGCGGGGCTGCACTTGGACTAGGAAAGGCTTAGAGATCAAGCCTTATCGAACTAAGACCCTCGCACTTAAATGTCCGGGGTTTTTTTTTg                                                                                         <  2:1010045/100‑1 (MQ=255)
      tatCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTg                                                                                         <  1:1179846/100‑1 (MQ=255)
       atCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGa                                                                                        >  2:544920/1‑100 (MQ=255)
         cccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGAcc                                                                                      >  1:1757957/1‑100 (MQ=255)
          ccACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCt                                                                                     >  2:826635/1‑100 (MQ=255)
           cACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCtt                                                                                    >  1:2067011/1‑100 (MQ=255)
              cGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTaaa                                                                                 <  2:1947562/100‑1 (MQ=255)
               gTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTaaaa                                                                                >  2:809167/1‑100 (MQ=255)
                tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTaaaaa                                                                               <  1:933859/100‑1 (MQ=255)
                tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTaaaaa                                                                               >  1:161355/1‑100 (MQ=255)
                tGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTaaaaa                                                                               >  1:47275/1‑100 (MQ=255)
                 gCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAAc                                                                              <  1:1149858/100‑1 (MQ=255)
                  cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACTTTAAAAACa                                                                             >  2:125613/1‑100 (MQ=255)
                  cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACa                                                                             <  2:642938/100‑1 (MQ=255)
                  cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACa                                                                             <  1:601319/100‑1 (MQ=255)
                  cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACa                                                                             <  2:213700/100‑1 (MQ=255)
                  cGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACa                                                                             >  1:1693702/1‑100 (MQ=255)
                   ggggCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACAt                                                                            <  1:2064616/100‑1 (MQ=255)
                    ggggTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGAACCGAAAGGTCCGGGGTTTTTTTTTTACCTTAAAAACAt                                                                            <  2:1298327/99‑1 (MQ=255)
                          cACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgag                                                                     <  2:1736121/100‑1 (MQ=255)
                           aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgagg                                                                    <  2:570350/100‑1 (MQ=255)
                           aCTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTATTTTTTTGACCTTAAAAACATAACCgagg                                                                    >  1:1781915/1‑100 (MQ=255)
                            cTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgagga                                                                   >  1:462344/1‑100 (MQ=255)
                            cTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgagga                                                                   >  1:1587634/1‑100 (MQ=255)
                            cTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgagga                                                                   >  1:1136655/1‑100 (MQ=255)
                            cTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgagg                                                                    <  1:1066870/99‑1 (MQ=255)
                             ttGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgagg                                                                    <  1:2021450/98‑1 (MQ=255)
                              tGGGCGGGGAAAGGCTTAGAGAGCAAGCGTTAACGAGGTAAGACCCGCGCACCGAAGGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                                 <  2:1172735/100‑1 (MQ=255)
                              tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                                 <  1:451743/100‑1 (MQ=255)
                              tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGc                                                                 <  2:828002/100‑1 (MQ=255)
                              tGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGAACAGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCgagg                                                                    <  2:2364567/97‑1 (MQ=255)
                               ggACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                                <  1:1153639/100‑1 (MQ=255)
                                gACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAg                                                               >  1:1991395/1‑100 (MQ=255)
                                   gAGGAAAGGCATAAAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCa                                                                <  2:2255884/96‑1 (MQ=255)
                                    aGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACaa                                                           >  1:787531/1‑100 (MQ=255)
                                     ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                          >  1:418268/1‑100 (MQ=255)
                                     ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                          <  2:1022697/100‑1 (MQ=255)
                                     ggAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAAt                                                          >  1:1206297/1‑100 (MQ=255)
                                      gAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGACAGGTCCGGGGTTTTTTTTTGACCTTAAAAACCTAACCAAGGAGCAGACCATg                                                         >  2:2089732/1‑100 (MQ=255)
                                      gAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                         <  2:1447899/100‑1 (MQ=255)
                                       aaaGGCTTAGAGCTCACGCCTTAGCGAACAAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATg                                                         <  1:2261096/99‑1 (MQ=255)
                                        acggCTTAGAGATAAAGACTTAACGAAATAAGCCACCCGTAGAGAAAGGCCCGGGGTTTTTTGTTGACCTTAAAAACATAACCTAGCAGCAGACTAt                                                          <  2:211991/95‑1 (MQ=255)
                                          ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATa                                                     <  2:80657/100‑1 (MQ=255)
                                          ggCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATa                                                     >  1:1559561/1‑100 (MQ=255)
                                             ttAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACa                                                  <  2:121627/100‑1 (MQ=255)
                                              tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                                 >  1:419551/1‑100 (MQ=255)
                                              tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                                 >  1:1686113/1‑100 (MQ=255)
                                              tAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                                 >  1:1646817/1‑100 (MQ=255)
                                               agagATCCAGCCTTACCGAGCGAGGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGc                                                <  1:826642/100‑1 (MQ=255)
                                               agagATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGc                                                >  1:33843/1‑100 (MQ=255)
                                               agagATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGc                                                >  1:2510689/1‑100 (MQ=255)
                                                gagaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAg                                                 >  1:2158748/1‑98 (MQ=255)
                                                gagaTCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGca                                               >  1:2451159/1‑100 (MQ=255)
                                                         ccTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAAACGAGGAACAGACAAAGAATAACAGCACAAAATTCt                                      >  1:752382/1‑100 (MQ=255)
                                                          cTTAACGAACTAAGACCCCCGAACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTg                                     <  2:1241812/100‑1 (MQ=255)
                                                           ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTTTTTTATAAACATAACCGAGGAGCAGACAATGAATACCAGCACAAAATTCTGt                                    >  2:1497290/1‑100 (MQ=255)
                                                           ttAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGt                                    >  1:1795969/1‑100 (MQ=255)
                                                            tAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAGATTCTGtt                                   <  2:81356/100‑1 (MQ=255)
                                                            tAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGtt                                   >  1:1598728/1‑100 (MQ=255)
                                                             aaCGAACTCAGCCCCCCGCAGTGAAAGGACCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTc                                       <  1:2412729/95‑1 (MQ=255)
                                                              aCGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtc                                 <  1:747855/100‑1 (MQ=255)
                                                               cGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTTTCCTTAATAACATAACCGGGGAGCAGACAATGCATAACAGCACAAAATTCTGTTtct                                >  2:1113914/1‑100 (MQ=255)
                                                               cGAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtct                                >  2:1425313/1‑100 (MQ=255)
                                                                gAACTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTtctc                               >  2:2459447/1‑100 (MQ=255)
                                                                 aaCTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                              >  1:555286/1‑100 (MQ=255)
                                                                 aaCGAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCa                              <  2:1206289/100‑1 (MQ=255)
                                                                   cTAAGACCCCCGCACCGAAATGTCCGGGGGTTTTTTTTGACCTTAAGAACCTAACCGAGTAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCaa                             <  1:574473/99‑1 (MQ=255)
                                                                   cTAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAg                            >  1:812796/1‑100 (MQ=255)
                                                                    tAAGACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGa                           >  1:2439342/1‑100 (MQ=255)
                                                                      aGACCCCCGCACCGAAAGATCCGAGGATTTTGTTTGACCTTAAAAACATAACGCACGGGAAGATAATGAATAACAGAAAAAACTACTGTTtgtcgatcat                         >  2:964055/1‑94 (MQ=255)
                                                                       ggcccccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTc                        <  1:1990448/98‑1 (MQ=255)
                                                                       gACCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTc                        <  2:300202/100‑1 (MQ=255)
                                                                        aCCCCCGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCa                       >  2:647230/1‑100 (MQ=255)
                                                                         ccccgGCACCAAAAGGCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGTGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCa                       <  2:1271482/98‑1 (MQ=255)
                                                                         cccccGCACCGAAAGGTCCGGGGTTTTTTTTTTTCTTTTTTTTCATAACCGAAGACCAGACAATGAATAACACCACAAAATTCTGGTTCTCAAGAtt                         >  2:609910/1‑97 (MQ=255)
                                                                         cccccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGAtt                         <  2:631974/97‑1 (MQ=255)
                                                                         cccccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAg                      >  2:1981432/1‑100 (MQ=255)
                                                                         cccccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAg                      >  1:486985/1‑100 (MQ=255)
                                                                          ccccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAgg                     <  2:1296790/100‑1 (MQ=255)
                                                                           cccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGa                    >  1:2169463/1‑100 (MQ=255)
                                                                           cccGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGa                    >  1:1041493/1‑100 (MQ=255)
                                                                             cGCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACg                  >  1:2501723/1‑100 (MQ=255)
                                                                              gCACCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgg                 >  1:1420458/1‑100 (MQ=255)
                                                                                aCCGAAAGGTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACgggg               <  2:812789/100‑1 (MQ=255)
                                                                                 cggAAAGGTCGGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGa              <  2:866295/98‑1 (MQ=255)
                                                                                 ccGAAAGGTCCGGGGTTTTTTTTTTTTTTTTTAAACATAAAGGAGGGGCAGACAATGAATCTCTGCACAGAGTTCTGTTTCTTTAGTTTCAggtc                   >  2:266561/1‑93 (MQ=255)
                                                                                      aGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaa         >  1:2509251/1‑100 (MQ=255)
                                                                                       ggTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGaactaac        <  1:47826/100‑1 (MQ=255)
                                                                                        gTCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCGGGACGGGGaactaact       >  1:2518799/1‑100 (MQ=255)
                                                                                         tCCGGGGTTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAactaacta      <  2:2279238/100‑1 (MQ=255)
                                                                                            ggggtTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGa   >  1:1659081/1‑100 (MQ=255)
                                                                                            ggggtTTTTTTATTACCTTAAAAACATAACCGAGGAGCAGACAATGAAGAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACGAACTATGa   <  2:1424047/100‑1 (MQ=255)
                                                                                             ggggTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGaa  <  2:2016029/100‑1 (MQ=255)
                                                                                                |                                                                                                
GATTATTATCCCACCGTGCGGGGCTGCACTTGGACGAGGAAAGGCTTAGAGATCAAGCCTTAACGAACTAAGACCCCCGCACCGAAAGGTCCGGGGGTTTTTTTTGACCTTAAAAACATAACCGAGGAGCAGACAATGAATAACAGCACAAAATTCTGTTTCTCAAGATTCAGGACGGGGAACTAACTATGAA  >  NC_000913/3950372‑3950564

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: