Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 416,389 (C)5→4 100% coding (564/1203 nt) sbcD ← exonuclease, dsDNA, ATP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913416,3850C.96.7% 394.8 / 5.2 90coding (568/1203 nt)sbcDexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base . (41/46);  total (41/49)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGG  >  NC_000913/416287‑416480
                                                                                                  |                                                                                               
atGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGcc                                                                                                >  2:648853/1‑100 (MQ=255)
atGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGcc                                                                                                <  2:276249/100‑1 (MQ=255)
  gCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGcc                                                                                                <  1:1790727/98‑1 (MQ=255)
   cccGAGCGCGCTGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCCGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGcccc                                                                                              >  2:751321/1‑99 (MQ=255)
   cccGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGcccca                                                                                             >  2:96482/1‑99 (MQ=255)
   cccGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGcccca                                                                                             >  2:915354/1‑99 (MQ=255)
        gcgcGATGTAGTCGGCTCGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGAGTGCTAATATAAATGTCACGCACGGCGTCACTTTTACTTG‑CCCCACGGt                                                                                        <  2:2232468/99‑1 (MQ=255)
         cgcgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCg                                                                                      >  2:814725/1‑100 (MQ=255)
         cgcgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCg                                                                                      >  2:1358910/1‑100 (MQ=255)
         cgcgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCg                                                                                      <  1:182179/100‑1 (MQ=255)
          gcgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGt                                                                                     >  1:1168273/1‑100 (MQ=255)
           cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGtt                                                                                    >  1:1049793/1‑100 (MQ=255)
           cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGtt                                                                                    <  2:81498/100‑1 (MQ=255)
                 aGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATg                                                                               <  2:606983/99‑1 (MQ=255)
                  gTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTc                                                                             >  2:991821/1‑100 (MQ=255)
                  gTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTc                                                                             >  1:2140297/1‑100 (MQ=255)
                       cTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTg                                                                        >  1:1989313/1‑100 (MQ=255)
                         ggtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGc                                                                      >  1:673353/1‑100 (MQ=255)
                         ggtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGc                                                                      >  1:2419354/1‑100 (MQ=255)
                          gtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGc                                                                      >  2:966143/1‑99 (MQ=255)
                          gtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCg                                                                     >  1:7600/1‑100 (MQ=255)
                            ggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAAAATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCgat                                                                   >  2:1366310/1‑100 (MQ=255)
                              aaaGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCgatga                                                                 >  1:8835/1‑100 (MQ=255)
                              aaaGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGGCGTTAAATGTCCCGTGGCgatg                                                                  <  2:618835/99‑1 (MQ=255)
                                aGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATATATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGatgatg                                                               <  2:1822982/100‑1 (MQ=255)
                                 gTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGatgatg                                                               <  2:1784833/100‑1 (MQ=255)
                                 gTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATgg                                                              <  1:1123408/100‑1 (MQ=255)
                                  ttttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATggg                                                             >  2:2016797/1‑100 (MQ=255)
                                   tttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGc                                                            <  1:1445863/100‑1 (MQ=255)
                                   tttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGc                                                            <  1:2025012/100‑1 (MQ=255)
                                   tttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGc                                                            <  1:1990620/100‑1 (MQ=255)
                                   tttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGc                                                            <  1:1989879/100‑1 (MQ=255)
                                   tttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGc                                                            <  1:1989787/100‑1 (MQ=255)
                                   tttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAAGGTCCCGTGGCGATGATGGGc                                                            <  2:1563361/100‑1 (MQ=255)
                                     tgtgCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCag                                                          >  1:976270/1‑100 (MQ=255)
                                     tgtgCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCAAGGTCGTTAAATGTCCCGTGGCGATGATGGGCag                                                          >  2:982210/1‑100 (MQ=255)
                                      gtgCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGCGGCGATGATGGGCaga                                                         <  2:1008117/100‑1 (MQ=255)
                                        gctggAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAgg                                                       >  2:986828/4‑100 (MQ=255)
                                        gCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAgg                                                       >  1:2083744/1‑100 (MQ=255)
                                        gCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAgg                                                       <  2:2083730/100‑1 (MQ=255)
                                        gCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAgg                                                       >  1:1129731/1‑100 (MQ=255)
                                         ccGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGc                                                      >  1:2349712/1‑100 (MQ=255)
                                         ccGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGc                                                      <  2:949606/100‑1 (MQ=255)
                                          cGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCt                                                     >  1:633553/1‑100 (MQ=255)
                                              aaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATc                                                 <  1:1623035/100‑1 (MQ=255)
                                              aaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCAATTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATc                                                 <  2:1794622/100‑1 (MQ=255)
                                                cgcgTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATcgc                                               >  1:2058798/1‑100 (MQ=255)
                                                     ccAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCgcgc                                          >  2:2021928/1‑100 (MQ=255)
                                                       aGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAg                                        <  1:924852/100‑1 (MQ=255)
                                                         cGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAAAGTCCCGGGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGtt                                      >  2:1904878/1‑100 (MQ=255)
                                                          gTGCCAATATAAATGTCACGCACGGCTTCACTTTTACTGG‑CGCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGttt                                     <  2:636342/100‑1 (MQ=255)
                                                           tGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTg                                    >  2:2027582/1‑100 (MQ=255)
                                                           tGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTg                                    >  1:1604832/1‑100 (MQ=255)
                                                            gCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGc                                   <  2:1989299/100‑1 (MQ=255)
                                                            gCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGc                                   <  1:1192904/100‑1 (MQ=255)
                                                            gCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑ACCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGc                                   <  2:1391611/100‑1 (MQ=255)
                                                             ccAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCa                                  >  1:2218148/1‑100 (MQ=255)
                                                             ccAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCa                                  >  1:137792/1‑100 (MQ=255)
                                                              cAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGc                                   <  2:288582/98‑1 (MQ=255)
                                                              cAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAg                                 <  2:1974140/100‑1 (MQ=255)
                                                              cAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAg                                 >  1:2067478/1‑100 (MQ=255)
                                                              cAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAg                                 <  1:2157858/100‑1 (MQ=255)
                                                              cAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAg                                 <  1:792550/100‑1 (MQ=255)
                                                               aaTATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAgg                                >  2:939284/1‑100 (MQ=255)
                                                               aaTATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAgg                                <  1:2489054/100‑1 (MQ=255)
                                                                   taAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATc                            >  1:1389734/1‑100 (MQ=255)
                                                                   taAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATc                            >  1:2447077/1‑100 (MQ=255)
                                                                   taAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATc                            >  1:1387569/1‑100 (MQ=255)
                                                                   taAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATc                            >  1:748951/1‑100 (MQ=255)
                                                                   taAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATc                            >  1:944454/1‑100 (MQ=255)
                                                                     aaTGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCgg                          <  2:498834/100‑1 (MQ=255)
                                                                      aTGTCACGAACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGc                         <  2:109203/100‑1 (MQ=255)
                                                                       tGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCa                        <  1:1616976/100‑1 (MQ=255)
                                                                         tCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATa                      <  1:572708/100‑1 (MQ=255)
                                                                          cacgcacgGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAg                     <  1:328135/100‑1 (MQ=255)
                                                                          cacgcacgGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGAGGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAg                     >  2:1640137/1‑100 (MQ=255)
                                                                           acgcacgGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAgt                    <  2:423679/100‑1 (MQ=255)
                                                                           acgcacgGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAgt                    >  2:926819/1‑100 (MQ=255)
                                                                            cgcacgGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAgtg                   >  2:1443655/1‑100 (MQ=255)
                                                                               acgGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGtgttg                <  2:2216979/100‑1 (MQ=255)
                                                                                 ggCGTCACTTTTACTGGCCCCCACGGTCGTAAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGtgttgt               <  2:1469934/100‑1 (MQ=255)
                                                                                  gCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGGTGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGttgttg             <  2:1162671/100‑1 (MQ=255)
                                                                                  gCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGttgttg             >  1:2513201/1‑100 (MQ=255)
                                                                                   cGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTgg            <  1:2541502/100‑1 (MQ=255)
                                                                                    gTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGt           <  2:1857800/100‑1 (MQ=255)
                                                                                      cACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGACGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaa         >  2:1363408/1‑100 (MQ=255)
                                                                                        cTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaata       >  2:1505830/1‑100 (MQ=255)
                                                                                        cTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaata       <  2:2513186/100‑1 (MQ=255)
                                                                                          tttAGTGGCCCGC‑CGGTTGTTTAATGTCCCGTGGCGACGATTGGCAGCGTCTGGTCGCCGCGCAGTTTGCTGGCATCTGCAGAGTGTTGTTGGtaa         <  2:212418/96‑1 (MQ=255)
                                                                                           ttACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATc    <  2:897278/100‑1 (MQ=255)
                                                                                           ttACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATc    <  1:2093112/100‑1 (MQ=255)
                                                                                           ttACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATc    <  1:1046485/100‑1 (MQ=255)
                                                                                            tACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGGTCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaat        <  2:2176258/95‑1 (MQ=255)
                                                                                             aCTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCgg  <  1:1363417/100‑1 (MQ=255)
                                                                                             aCTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCgg  <  2:1936280/100‑1 (MQ=255)
                                                                                             aCTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCgg  <  1:2024115/100‑1 (MQ=255)
                                                                                                  |                                                                                               
ATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGG  >  NC_000913/416287‑416480

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: ‑