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)
        gcgcGATGTAGTCGGCTCGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGAGTGCTAATATAAATGTCACGCACGGCGTCACTTTTACTTGCCCCACGGt                                                                                        <  2:2232468/99‑1 (MQ=255)
         cgcgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCg                                                                                      >  2:814725/1‑100 (MQ=255)
         cgcgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCg                                                                                      >  2:1358910/1‑100 (MQ=255)
         cgcgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCg                                                                                      <  1:182179/100‑1 (MQ=255)
          gcgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGt                                                                                     >  1:1168273/1‑100 (MQ=255)
           cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGtt                                                                                    >  1:1049793/1‑100 (MQ=255)
           cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGtt                                                                                    <  2:81498/100‑1 (MQ=255)
                 aGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATg                                                                               <  2:606983/99‑1 (MQ=255)
                  gTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTc                                                                             >  2:991821/1‑100 (MQ=255)
                  gTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTc                                                                             >  1:2140297/1‑100 (MQ=255)
                       cTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTg                                                                        >  1:1989313/1‑100 (MQ=255)
                         ggtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGc                                                                      >  1:673353/1‑100 (MQ=255)
                         ggtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGc                                                                      >  1:2419354/1‑100 (MQ=255)
                          gtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGc                                                                      >  2:966143/1‑99 (MQ=255)
                          gtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCg                                                                     >  1:7600/1‑100 (MQ=255)
                            ggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAAAATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCgat                                                                   >  2:1366310/1‑100 (MQ=255)
                              aaaGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCgatga                                                                 >  1:8835/1‑100 (MQ=255)
                              aaaGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGGCGTTAAATGTCCCGTGGCgatg                                                                  <  2:618835/99‑1 (MQ=255)
                                aGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATATATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGatgatg                                                               <  2:1822982/100‑1 (MQ=255)
                                 gTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGatgatg                                                               <  2:1784833/100‑1 (MQ=255)
                                 gTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATgg                                                              <  1:1123408/100‑1 (MQ=255)
                                  ttttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATggg                                                             >  2:2016797/1‑100 (MQ=255)
                                   tttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGc                                                            <  1:1445863/100‑1 (MQ=255)
                                   tttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGc                                                            <  1:2025012/100‑1 (MQ=255)
                                   tttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGc                                                            <  1:1990620/100‑1 (MQ=255)
                                   tttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGc                                                            <  1:1989879/100‑1 (MQ=255)
                                   tttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGc                                                            <  1:1989787/100‑1 (MQ=255)
                                   tttGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAAGGTCCCGTGGCGATGATGGGc                                                            <  2:1563361/100‑1 (MQ=255)
                                     tgtgCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCag                                                          >  1:976270/1‑100 (MQ=255)
                                     tgtgCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCAAGGTCGTTAAATGTCCCGTGGCGATGATGGGCag                                                          >  2:982210/1‑100 (MQ=255)
                                      gtgCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGCGGCGATGATGGGCaga                                                         <  2:1008117/100‑1 (MQ=255)
                                        gctggAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAgg                                                       >  2:986828/4‑100 (MQ=255)
                                        gCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAgg                                                       >  1:2083744/1‑100 (MQ=255)
                                        gCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAgg                                                       <  2:2083730/100‑1 (MQ=255)
                                        gCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAgg                                                       >  1:1129731/1‑100 (MQ=255)
                                         ccGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGc                                                      >  1:2349712/1‑100 (MQ=255)
                                         ccGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGc                                                      <  2:949606/100‑1 (MQ=255)
                                          cGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCt                                                     >  1:633553/1‑100 (MQ=255)
                                              aaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATc                                                 <  1:1623035/100‑1 (MQ=255)
                                              aaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCAATTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATc                                                 <  2:1794622/100‑1 (MQ=255)
                                                cgcgTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATcgc                                               >  1:2058798/1‑100 (MQ=255)
                                                     ccAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCgcgc                                          >  2:2021928/1‑100 (MQ=255)
                                                       aGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAg                                        <  1:924852/100‑1 (MQ=255)
                                                         cGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAAAGTCCCGGGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGtt                                      >  2:1904878/1‑100 (MQ=255)
                                                          gTGCCAATATAAATGTCACGCACGGCTTCACTTTTACTGGCGCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGttt                                     <  2:636342/100‑1 (MQ=255)
                                                           tGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTg                                    >  2:2027582/1‑100 (MQ=255)
                                                           tGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTg                                    >  1:1604832/1‑100 (MQ=255)
                                                            gCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGc                                   <  2:1989299/100‑1 (MQ=255)
                                                            gCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGc                                   <  1:1192904/100‑1 (MQ=255)
                                                            gCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGACCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGc                                   <  2:1391611/100‑1 (MQ=255)
                                                             ccAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCa                                  >  1:2218148/1‑100 (MQ=255)
                                                             ccAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCa                                  >  1:137792/1‑100 (MQ=255)
                                                              cAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGc                                   <  2:288582/98‑1 (MQ=255)
                                                              cAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAg                                 <  2:1974140/100‑1 (MQ=255)
                                                              cAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAg                                 >  1:2067478/1‑100 (MQ=255)
                                                              cAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAg                                 <  1:2157858/100‑1 (MQ=255)
                                                              cAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAg                                 <  1:792550/100‑1 (MQ=255)
                                                               aaTATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAgg                                >  2:939284/1‑100 (MQ=255)
                                                               aaTATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAgg                                <  1:2489054/100‑1 (MQ=255)
                                                                   taAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATc                            >  1:1389734/1‑100 (MQ=255)
                                                                   taAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATc                            >  1:2447077/1‑100 (MQ=255)
                                                                   taAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATc                            >  1:1387569/1‑100 (MQ=255)
                                                                   taAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATc                            >  1:748951/1‑100 (MQ=255)
                                                                   taAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATc                            >  1:944454/1‑100 (MQ=255)
                                                                     aaTGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCgg                          <  2:498834/100‑1 (MQ=255)
                                                                      aTGTCACGAACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGc                         <  2:109203/100‑1 (MQ=255)
                                                                       tGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCa                        <  1:1616976/100‑1 (MQ=255)
                                                                         tCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATa                      <  1:572708/100‑1 (MQ=255)
                                                                          cacgcacgGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAg                     <  1:328135/100‑1 (MQ=255)
                                                                          cacgcacgGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGAGGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAg                     >  2:1640137/1‑100 (MQ=255)
                                                                           acgcacgGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAgt                    <  2:423679/100‑1 (MQ=255)
                                                                           acgcacgGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAgt                    >  2:926819/1‑100 (MQ=255)
                                                                            cgcacgGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAgtg                   >  2:1443655/1‑100 (MQ=255)
                                                                               acgGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGtgttg                <  2:2216979/100‑1 (MQ=255)
                                                                                 ggCGTCACTTTTACTGGCCCCCACGGTCGTAAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGtgttgt               <  2:1469934/100‑1 (MQ=255)
                                                                                  gCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGGTGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGttgttg             <  2:1162671/100‑1 (MQ=255)
                                                                                  gCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGttgttg             >  1:2513201/1‑100 (MQ=255)
                                                                                   cGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTgg            <  1:2541502/100‑1 (MQ=255)
                                                                                    gTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGt           <  2:1857800/100‑1 (MQ=255)
                                                                                      cACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGACGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaa         >  2:1363408/1‑100 (MQ=255)
                                                                                        cTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaata       >  2:1505830/1‑100 (MQ=255)
                                                                                        cTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaata       <  2:2513186/100‑1 (MQ=255)
                                                                                          tttAGTGGCCCGCCGGTTGTTTAATGTCCCGTGGCGACGATTGGCAGCGTCTGGTCGCCGCGCAGTTTGCTGGCATCTGCAGAGTGTTGTTGGtaa         <  2:212418/96‑1 (MQ=255)
                                                                                           ttACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATc    <  2:897278/100‑1 (MQ=255)
                                                                                           ttACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATc    <  1:2093112/100‑1 (MQ=255)
                                                                                           ttACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATc    <  1:1046485/100‑1 (MQ=255)
                                                                                            tACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGGTCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaat        <  2:2176258/95‑1 (MQ=255)
                                                                                             aCTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCgg  <  1:1363417/100‑1 (MQ=255)
                                                                                             aCTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCgg  <  2:1936280/100‑1 (MQ=255)
                                                                                             aCTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCgg  <  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: