Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 416,389 (C)5→4 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.92.9% 353.0 / 23.3 86coding (568/1203 nt)sbcDexonuclease, dsDNA, ATP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (3/3);  new base . (36/43);  total (39/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAAT  >  NC_000913/416286‑416484
                                                                                                   |                                                                                                   
tatGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGc                                                                                                     <  2:262183/100‑1 (MQ=255)
 atGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGTCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGcc                                                                                                    <  2:242121/100‑1 (MQ=255)
 atGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGcc                                                                                                    <  2:246270/100‑1 (MQ=255)
  tGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGccc                                                                                                   >  2:1483693/1‑100 (MQ=255)
   gCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGcccc                                                                                                  >  1:346243/1‑100 (MQ=255)
     ccGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGccccac                                                                                                >  1:780096/1‑98 (MQ=255)
      cGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGccccacg                                                                                               >  2:1409534/1‑97 (MQ=255)
          cgcgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCg                                                                                          <  2:1045937/100‑1 (MQ=255)
            cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAATGCGTCTAGCGTGCCAATATAAATGTTACGCATGGTGTTACTTTTACTGGCCCCTCGGTCGtt                                                                                        >  2:1616871/1‑100 (MQ=255)
            cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGtt                                                                                        <  1:219123/100‑1 (MQ=255)
            cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGtt                                                                                        <  1:896851/100‑1 (MQ=255)
            cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGtt                                                                                        >  2:1916307/1‑100 (MQ=255)
            cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGtt                                                                                        >  2:125500/1‑100 (MQ=255)
            cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGtt                                                                                        <  1:1565339/100‑1 (MQ=255)
             gATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTa                                                                                       >  2:127500/1‑100 (MQ=255)
              aTGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTaa                                                                                      <  2:674824/100‑1 (MQ=255)
              aTGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTaa                                                                                      >  1:694989/1‑100 (MQ=255)
                gTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGAGTCACTTTTACTGGCCCCACGGTCGTTAAAt                                                                                    >  2:1697802/1‑100 (MQ=255)
                 tAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATg                                                                                   <  2:640304/100‑1 (MQ=255)
                 tAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATg                                                                                   <  2:1744091/100‑1 (MQ=255)
                 tAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATg                                                                                   >  2:1720173/1‑100 (MQ=255)
                 tAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCCTTAAATg                                                                                   <  2:1022124/100‑1 (MQ=255)
                   gTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTc                                                                                 >  1:204411/1‑100 (MQ=255)
                   gTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTc                                                                                 >  2:1689282/1‑100 (MQ=255)
                   gTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTc                                                                                 >  2:1372285/1‑100 (MQ=255)
                         tggtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTgg                                                                           <  2:777083/100‑1 (MQ=255)
                          ggtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGc                                                                          <  1:1663025/100‑1 (MQ=255)
                          ggtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGc                                                                          <  2:1513849/100‑1 (MQ=255)
                          ggtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGc                                                                          <  1:1720162/100‑1 (MQ=255)
                           gtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCg                                                                         <  2:1703405/100‑1 (MQ=255)
                             ggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCgat                                                                       <  2:1096152/100‑1 (MQ=255)
                             ggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCgat                                                                       >  2:1392473/1‑100 (MQ=255)
                              gAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCg                                                                         >  1:367121/1‑98 (MQ=255)
                               aaaGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCgatga                                                                     >  1:355128/1‑100 (MQ=255)
                               aaaGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCgatga                                                                     >  1:1706755/1‑100 (MQ=255)
                               aaaGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCgatga                                                                     >  2:153503/1‑100 (MQ=255)
                               aaaGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGGCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgatg                                                                      <  1:295712/100‑1 (MQ=255)
                                aaGTTTTGTGTCGAAAATGCGTCCATCATGCAAAGAGACATGTCACGCACGGCGACACTTTTACTGGCCCCACGGTCGTTGAATGTCCCGTGGCgatg                                                                      <  1:871541/98‑1 (MQ=255)
                                 aGTTTTATGCCGGAAACGCGTCCAACGTGCCAATATAAATGTCACGCACGACGTCACTGTTACTTGCCCCACGGGCGGTAAATGTCCCGTGGCTatgatg                                                                   >  1:1621579/1‑100 (MQ=255)
                                  gTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATgg                                                                  >  2:1833312/1‑100 (MQ=255)
                                      tgtgCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCag                                                              <  1:1916294/100‑1 (MQ=255)
                                        tgCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCagag                                                            >  2:1044444/1‑100 (MQ=255)
                                         gCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACTGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAgg                                                           <  2:1623941/100‑1 (MQ=255)
                                         gCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAgg                                                           >  1:1211276/1‑100 (MQ=255)
                                          ccGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGc                                                          >  2:854791/1‑100 (MQ=255)
                                          ccGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGc                                                          >  2:313475/1‑100 (MQ=255)
                                           cGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCt                                                         >  2:119005/1‑100 (MQ=255)
                                              aaaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGAt                                                      <  2:1012246/100‑1 (MQ=255)
                                              aaaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGAt                                                      <  2:1781590/100‑1 (MQ=255)
                                               aaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATc                                                     <  1:127499/100‑1 (MQ=255)
                                               aaCGCGTCCAACGTGCCAATATAAATCTCACACACACCATCACTTTTACTAACCCCACAATCGTTAAATATCCCATAGCAATAATAAGCAGAGGCTGATc                                                     <  1:1616861/100‑1 (MQ=255)
                                                aCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATcg                                                    <  2:912139/100‑1 (MQ=255)
                                                aCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATcg                                                    <  1:1184626/100‑1 (MQ=255)
                                                    gTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATcgccgc                                                >  1:1224900/1‑100 (MQ=255)
                                                    gTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATcgccgc                                                >  2:1810123/1‑100 (MQ=255)
                                                    gTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATcgccgc                                                <  2:426688/100‑1 (MQ=255)
                                                     tCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCcgcg                                               <  2:1760885/100‑1 (MQ=255)
                                                          cGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGtt                                          <  1:1831132/100‑1 (MQ=255)
                                                            tGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTg                                        <  1:976813/100‑1 (MQ=255)
                                                             gCCAATATGAATGTCACGCACGGCGTCAGTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGc                                       <  1:481311/100‑1 (MQ=255)
                                                              ccAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCa                                      >  2:10964/1‑100 (MQ=255)
                                                              ccAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCa                                      <  1:905213/100‑1 (MQ=255)
                                                               cAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAg                                     <  2:1227936/100‑1 (MQ=255)
                                                                aaTATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAg                                     >  2:155493/1‑100 (MQ=255)
                                                                aaTATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAgg                                    >  2:136075/1‑100 (MQ=255)
                                                                aaTATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAgg                                    >  2:257892/1‑100 (MQ=255)
                                                                   ataAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCAt                                 >  1:1394023/1‑100 (MQ=255)
                                                                    taAATGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATc                                >  1:1788246/1‑100 (MQ=255)
                                                                      aaTGTCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCgg                              <  2:355130/100‑1 (MQ=255)
                                                                          tCACGCACGGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATa                          >  2:1713728/1‑100 (MQ=255)
                                                                           cacgcacgGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAg                         >  1:1478233/1‑100 (MQ=255)
                                                                           cacgcacgGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAg                         <  2:1112051/100‑1 (MQ=255)
                                                                            acgcacgGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAgt                        >  1:1428242/1‑100 (MQ=255)
                                                                            acgcacgGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAgt                        >  2:1268159/1‑100 (MQ=255)
                                                                             cgcacgGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAgt                        >  2:569877/1‑100 (MQ=255)
                                                                               cacgGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAgtgt                      <  2:1436748/100‑1 (MQ=255)
                                                                               cacgGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGtgtt                     <  1:263925/100‑1 (MQ=255)
                                                                                acgGCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGtgttg                    <  2:1706766/100‑1 (MQ=255)
                                                                                   gCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGttgttg                 <  2:1880008/100‑1 (MQ=255)
                                                                                   gCGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGttgttg                 <  1:644589/100‑1 (MQ=255)
                                                                                    cGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTgg                <  1:1689272/100‑1 (MQ=255)
                                                                                    cGTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTgg                <  1:1035049/100‑1 (MQ=255)
                                                                                     gTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGt               <  1:4431/100‑1 (MQ=255)
                                                                                     gTCACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGt               <  1:1332448/100‑1 (MQ=255)
                                                                                       cACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGta              <  2:1034711/100‑1 (MQ=255)
                                                                                       cACTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaa             >  2:1062639/1‑100 (MQ=255)
                                                                                        aCTTGTACTGGCCCCACGGTCGTTAAATGGCCCGTAGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaat            <  1:1462646/100‑1 (MQ=255)
                                                                                         cTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaata           >  1:1766153/1‑100 (MQ=255)
                                                                                         cTTTTACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaata           >  1:456319/1‑100 (MQ=255)
                                                                                          ttttACTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaataa          <  1:313473/100‑1 (MQ=255)
                                                                                              aCTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCgg      >  2:1590530/1‑100 (MQ=255)
                                                                                              aCTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCgg      <  2:1073221/100‑1 (MQ=255)
                                                                                              aCTGGCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCgg      <  2:1176269/100‑1 (MQ=255)
                                                                                                   cccccACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAAt  >  1:504302/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
TATGCCCGAGCGCGATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCGGTAAT  >  NC_000913/416286‑416484

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