Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,679,503 G→A F1081F (TTC→TTT yphG ← DUF4380 domain‑containing TPR repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,679,5030GA96.7% 289.6 / ‑1.4 90F1081F (TTC→TTTyphGDUF4380 domain‑containing TPR repeat protein
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (46/41);  total (48/42)
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

TAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAG  >  NC_000913/2679407‑2679598
                                                                                                 |                                                                                               
tAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAaat                                                                                              >  1:940579/1‑100 (MQ=255)
tAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACACGGTTAAAAATAAAAGGCATCACGGTGGTAaat                                                                                              >  2:110724/1‑100 (MQ=255)
 aaTTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAaata                                                                                             <  2:693275/100‑1 (MQ=255)
    ctgtaatatttgaGCGGCACGCGCTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATg                                                                                          >  2:1741464/11‑100 (MQ=25)
     attTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAATTAATAAGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGc                                                                                         >  1:1042218/1‑100 (MQ=255)
      ttttGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGcc                                                                                        <  2:1594255/100‑1 (MQ=255)
      ttttGTGATGAGCGGCAGCGCGTGCAGCTCAAATATTACAGTTAAGAGGGTTTAGTGAAAGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGc                                                                                         <  1:1347885/99‑1 (MQ=255)
       tttGTTATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCg                                                                                       <  2:197171/100‑1 (MQ=255)
        ttGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGATGAATAATGCCGc                                                                                      >  2:10862/1‑100 (MQ=255)
           tgatgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCttt                                                                                   >  2:881585/1‑100 (MQ=255)
           tgatgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCttt                                                                                   >  2:172108/1‑100 (MQ=255)
           tgatgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCttt                                                                                   >  1:1574292/1‑100 (MQ=255)
           tgatgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCtt                                                                                    >  2:488894/1‑99 (MQ=255)
            tttgaGCGGCACGCGCTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTc                                                                                  >  2:470070/3‑100 (MQ=255)
            gatgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTc                                                                                  <  2:1495379/100‑1 (MQ=255)
             ttgaGCGGCACGCGCTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCg                                                                                 >  2:1027393/2‑100 (MQ=255)
             attaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCg                                                                                 >  1:336389/4‑100 (MQ=255)
              tgactggCAGCGCGTTCTGCTCAAATATTACAGTTAATAGGGTTTAGTGTACTCAGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCg                                                                                 <  1:445598/94‑1 (MQ=255)
              tgaGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCgg                                                                                <  1:913497/100‑1 (MQ=255)
                 gCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCa                                                                             <  2:390984/100‑1 (MQ=255)
                 gCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCa                                                                             >  1:1409500/1‑100 (MQ=255)
                 gCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCa                                                                             >  2:875065/1‑100 (MQ=255)
                 gCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCa                                                                             >  2:621240/1‑100 (MQ=255)
                  cGGCAGCGCGTGCTGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCa                                                                             <  2:1510128/99‑1 (MQ=255)
                  cGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGTTAAATAATGCCGCTTTCGGc                                                                               <  1:376351/97‑1 (MQ=255)
                   ggCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCaa                                                                            >  1:244437/1‑99 (MQ=255)
                   ggCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAg                                                                           <  2:787029/100‑1 (MQ=255)
                    gtaGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACTCGGTTAAAAATAAAAGGCTTCACGGTGGTAAATAATGCCGCTTTCGGcc                                                                              <  1:522424/94‑1 (MQ=255)
                    gCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAg                                                                           <  1:1779032/99‑1 (MQ=255)
                      aGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCa                                                                             <  1:212743/95‑1 (MQ=255)
                       gcgcGTGCCGCTCAAATATTACAGATAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAATCCgg                                                                       >  1:422253/1‑100 (MQ=255)
                          cgTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGAtt                                                                    >  2:1747021/1‑100 (MQ=255)
                          cgTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGAtt                                                                    >  2:1613062/1‑100 (MQ=255)
                           gTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGAtt                                                                    <  1:881579/99‑1 (MQ=255)
                             tGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTaa                                                                  >  2:1507505/1‑100 (MQ=255)
                               tcgttCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATc                                                                <  1:1149409/96‑1 (MQ=255)
                               ccGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATc                                                                >  1:184844/1‑100 (MQ=255)
                                    cAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGt                                                           <  2:1317074/100‑1 (MQ=255)
                                    cAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGt                                                           >  1:1816846/1‑100 (MQ=255)
                                    cAAATATGACAGTTAATAGGGTTTTGTGAACGCGGTTAAAAATAAAAGGCATCACGGTTGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCg                                                               <  1:1603341/96‑1 (MQ=255)
                                     aaaTATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTc                                                          <  2:1504881/100‑1 (MQ=255)
                                     aaaTATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTCGGCCAAGCCGGATTAATCGCCGTc                                                          <  1:518305/99‑1 (MQ=255)
                                      aaTATTAAAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTc                                                          <  1:260449/99‑1 (MQ=255)
                                          ttACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCACAGTCAAt                                                       >  2:1141343/1‑98 (MQ=255)
                                           tACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCa                                                    >  2:1665230/1‑100 (MQ=255)
                                            aCAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAg                                                   <  2:505980/100‑1 (MQ=255)
                                             cgattAATAGGGTTTAGTTAACTCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAgc                                                  <  1:943010/97‑1 (MQ=255)
                                              aGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAgcg                                                 >  1:782780/1‑100 (MQ=255)
                                              aGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAgcg                                                 >  1:355559/1‑100 (MQ=255)
                                               gTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAgcgc                                                >  2:858379/1‑100 (MQ=255)
                                               gTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAgcgc                                                >  2:72581/1‑100 (MQ=255)
                                               gTTAATAGGGTTTAGTGAACGCGGATAAAAATAAAAGGCAGCACGGTGGTAAATAATGCCGTTTTCGGCCAAGCCGGAATAATCTCAGTCGATTCAgcg                                                 >  1:1307028/1‑99 (MQ=255)
                                                ttAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGcgcg                                               <  2:1299813/100‑1 (MQ=255)
                                                 tAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATACCGCTTTCGGCCAAGCCGGATTAATCGCGGTCAATTCAGCGCGa                                              >  1:322024/1‑100 (MQ=255)
                                                    tAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGc                                            <  1:725048/99‑1 (MQ=255)
                                                    tAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCa                                           >  2:809464/1‑100 (MQ=255)
                                                        gTTTAGTGAACGCGGTTAAAAATAAAAGGCTTTATGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACt                                       <  2:38817/100‑1 (MQ=255)
                                                        gTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACt                                       <  1:1417600/100‑1 (MQ=255)
                                                        gTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACt                                       >  1:1252628/1‑100 (MQ=255)
                                                        gTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACt                                       >  1:1671234/1‑100 (MQ=255)
                                                         tttAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGCTTAATCGCCGTCAATTCAGCGCGAGCAGAcac                                      >  2:703124/1‑98 (MQ=255)
                                                         tttAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTc                                      <  2:1090706/100‑1 (MQ=255)
                                                         tttAGTGAACGCCGTTAAAAATAAACGGCATCGCGCTGGTAAATAATGCCGCTTTCGGCCAACCCGTATTAATCGCCGTTAATTCAGCGCGAGCAGAATc                                      <  1:752645/100‑1 (MQ=255)
                                                         tttAGCGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCAGATTAATCGCCGTCAATTCAGCGCGAGCAGACTc                                      <  1:1507494/100‑1 (MQ=255)
                                                           tAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCat                                    >  1:126025/1‑100 (MQ=255)
                                                            aGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCata                                   >  1:39466/1‑100 (MQ=255)
                                                                aaCGCGGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgc                               <  2:756646/100‑1 (MQ=255)
                                                                  cgcgGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgcgc                             >  1:361355/1‑100 (MQ=255)
                                                                  cgcgGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGATTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTgcgc                             <  1:1271587/100‑1 (MQ=255)
                                                                    cgGTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCaa                           >  2:724809/1‑100 (MQ=255)
                                                                     ggTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAAc                          <  1:675475/100‑1 (MQ=255)
                                                                      gTTAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAAcc                         <  2:1353027/100‑1 (MQ=255)
                                                                       ttAAAAATAAAAGGCTTCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCTCGAGCAGACTCATATTGCGCa                            <  1:1540669/96‑1 (MQ=255)
                                                                       ttAAAAATAAAAGGCATCACGGTGGTATATAATGCCGCTTTCGGCCAAGCCGTATTAATCTCCTTCAATTTAGTTTTTGTAGACTCATATTGCGaaa                           <  1:448415/97‑1 (MQ=255)
                                                                       ttAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCt                        <  2:1750565/100‑1 (MQ=255)
                                                                       tgaaaaaTAAAAGACATCACGGTGGTAAATAGTGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCt                        >  2:1227026/3‑100 (MQ=255)
                                                                        tAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTc                       >  1:477806/1‑100 (MQ=255)
                                                                        tAAAAATAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTc                       >  2:1332838/1‑100 (MQ=255)
                                                                          aaaaTAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTccc                     >  2:139039/1‑100 (MQ=255)
                                                                           aaaTAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTcccc                    <  2:1115260/100‑1 (MQ=255)
                                                                           aaaTAAAAGGCATCACGGTAGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTAATATTGCGCAACATcccc                    <  1:1420369/100‑1 (MQ=255)
                                                                            aaTAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCg                   >  1:1247907/1‑100 (MQ=255)
                                                                              tAAAAGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAg                 <  1:1775261/100‑1 (MQ=255)
                                                                               aaaaGGCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGt                <  1:1741451/100‑1 (MQ=255)
                                                                                   ggCATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCcg            >  2:725962/1‑100 (MQ=255)
                                                                                     cATCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCcgcg          >  2:1368315/1‑100 (MQ=255)
                                                                                       tCACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGa        >  2:1400896/1‑100 (MQ=255)
                                                                                        cACGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGAc       <  1:772900/100‑1 (MQ=255)
                                                                                          cGGTGGTAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGc     >  2:720946/1‑100 (MQ=255)
                                                                                           ggtggtAAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGcc    >  2:1396176/1‑100 (MQ=255)
                                                                                             tggtaAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCa   <  1:101145/99‑1 (MQ=255)
                                                                                             tggtaAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAg  <  1:1027383/100‑1 (MQ=255)
                                                                                                 |                                                                                               
TAATTATTTTGTGATGAGCGGCAGCGCGTGCCGCTCAAATATTACAGTTAATAGGGTTTAGTGAACGCGGTTAAAAATAAAAGGCATCACGGTGGTGAATAATGCCGCTTTCGGCCAAGCCGGATTAATCGCCGTCAATTCAGCGCGAGCAGACTCATATTGCGCAACCTCCCCGAGTCCCGCGGACGCCAG  >  NC_000913/2679407‑2679598

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: