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)
          cgcgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCg                                                                                          <  2:1045937/100‑1 (MQ=255)
            cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAATGCGTCTAGCGTGCCAATATAAATGTTACGCATGGTGTTACTTTTACTGG‑CCCCTCGGTCGtt                                                                                        >  2:1616871/1‑100 (MQ=255)
            cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGtt                                                                                        <  1:219123/100‑1 (MQ=255)
            cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGtt                                                                                        <  1:896851/100‑1 (MQ=255)
            cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGtt                                                                                        >  2:1916307/1‑100 (MQ=255)
            cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGtt                                                                                        >  2:125500/1‑100 (MQ=255)
            cgATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGtt                                                                                        <  1:1565339/100‑1 (MQ=255)
             gATGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTa                                                                                       >  2:127500/1‑100 (MQ=255)
              aTGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTaa                                                                                      <  2:674824/100‑1 (MQ=255)
              aTGTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTaa                                                                                      >  1:694989/1‑100 (MQ=255)
                gTAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGAGTCACTTTTACTGG‑CCCCACGGTCGTTAAAt                                                                                    >  2:1697802/1‑100 (MQ=255)
                 tAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATg                                                                                   <  2:640304/100‑1 (MQ=255)
                 tAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATg                                                                                   <  2:1744091/100‑1 (MQ=255)
                 tAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATg                                                                                   >  2:1720173/1‑100 (MQ=255)
                 tAGTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCCTTAAATg                                                                                   <  2:1022124/100‑1 (MQ=255)
                   gTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTc                                                                                 >  1:204411/1‑100 (MQ=255)
                   gTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTc                                                                                 >  2:1689282/1‑100 (MQ=255)
                   gTCGGCTGGTGGAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTc                                                                                 >  2:1372285/1‑100 (MQ=255)
                         tggtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTgg                                                                           <  2:777083/100‑1 (MQ=255)
                          ggtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGc                                                                          <  1:1663025/100‑1 (MQ=255)
                          ggtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGc                                                                          <  2:1513849/100‑1 (MQ=255)
                          ggtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGc                                                                          <  1:1720162/100‑1 (MQ=255)
                           gtggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCg                                                                         <  2:1703405/100‑1 (MQ=255)
                             ggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCgat                                                                       <  2:1096152/100‑1 (MQ=255)
                             ggAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCgat                                                                       >  2:1392473/1‑100 (MQ=255)
                              gAAAGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCg                                                                         >  1:367121/1‑98 (MQ=255)
                               aaaGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCgatga                                                                     >  1:355128/1‑100 (MQ=255)
                               aaaGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCgatga                                                                     >  1:1706755/1‑100 (MQ=255)
                               aaaGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCgatga                                                                     >  2:153503/1‑100 (MQ=255)
                               aaaGTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGGCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCgatg                                                                      <  1:295712/100‑1 (MQ=255)
                                aaGTTTTGTGTCGAAAATGCGTCCATCATGCAAAGAGACATGTCACGCACGGCGACACTTTTACTGG‑CCCCACGGTCGTTGAATGTCCCGTGGCgatg                                                                      <  1:871541/98‑1 (MQ=255)
                                 aGTTTTATGCCGGAAACGCGTCCAACGTGCCAATATAAATGTCACGCACGACGTCACTGTTACTTG‑CCCCACGGGCGGTAAATGTCCCGTGGCTatgatg                                                                   >  1:1621579/1‑100 (MQ=255)
                                  gTTTTGTGCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATgg                                                                  >  2:1833312/1‑100 (MQ=255)
                                      tgtgCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCag                                                              <  1:1916294/100‑1 (MQ=255)
                                        tgCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCagag                                                            >  2:1044444/1‑100 (MQ=255)
                                         gCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACTGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAgg                                                           <  2:1623941/100‑1 (MQ=255)
                                         gCCGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAgg                                                           >  1:1211276/1‑100 (MQ=255)
                                          ccGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGc                                                          >  2:854791/1‑100 (MQ=255)
                                          ccGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGc                                                          >  2:313475/1‑100 (MQ=255)
                                           cGGAAACGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCt                                                         >  2:119005/1‑100 (MQ=255)
                                              aaaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGAt                                                      <  2:1012246/100‑1 (MQ=255)
                                              aaaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGAt                                                      <  2:1781590/100‑1 (MQ=255)
                                               aaCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATc                                                     <  1:127499/100‑1 (MQ=255)
                                               aaCGCGTCCAACGTGCCAATATAAATCTCACACACACCATCACTTTTACT‑AACCCCACAATCGTTAAATATCCCATAGCAATAATAAGCAGAGGCTGATc                                                     <  1:1616861/100‑1 (MQ=255)
                                                aCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATcg                                                    <  2:912139/100‑1 (MQ=255)
                                                aCGCGTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATcg                                                    <  1:1184626/100‑1 (MQ=255)
                                                    gTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATcgccgc                                                >  1:1224900/1‑100 (MQ=255)
                                                    gTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATcgccgc                                                >  2:1810123/1‑100 (MQ=255)
                                                    gTCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATcgccgc                                                <  2:426688/100‑1 (MQ=255)
                                                     tCCAGCGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCcgcg                                               <  2:1760885/100‑1 (MQ=255)
                                                          cGTGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGtt                                          <  1:1831132/100‑1 (MQ=255)
                                                            tGCCAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTg                                        <  1:976813/100‑1 (MQ=255)
                                                             gCCAATATGAATGTCACGCACGGCGTCAGTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGc                                       <  1:481311/100‑1 (MQ=255)
                                                              ccAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCa                                      >  2:10964/1‑100 (MQ=255)
                                                              ccAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCa                                      <  1:905213/100‑1 (MQ=255)
                                                               cAATATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAg                                     <  2:1227936/100‑1 (MQ=255)
                                                                aaTATAAATGTCACGCACGGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAg                                     >  2:155493/1‑100 (MQ=255)
                                                                aaTATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAgg                                    >  2:136075/1‑100 (MQ=255)
                                                                aaTATAAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAgg                                    >  2:257892/1‑100 (MQ=255)
                                                                   ataAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCAt                                 >  1:1394023/1‑100 (MQ=255)
                                                                    taAATGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATc                                >  1:1788246/1‑100 (MQ=255)
                                                                      aaTGTCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCgg                              <  2:355130/100‑1 (MQ=255)
                                                                          tCACGCACGGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATa                          >  2:1713728/1‑100 (MQ=255)
                                                                           cacgcacgGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAg                         >  1:1478233/1‑100 (MQ=255)
                                                                           cacgcacgGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAg                         <  2:1112051/100‑1 (MQ=255)
                                                                            acgcacgGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAgt                        >  1:1428242/1‑100 (MQ=255)
                                                                            acgcacgGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAgt                        >  2:1268159/1‑100 (MQ=255)
                                                                             cgcacgGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAgt                        >  2:569877/1‑100 (MQ=255)
                                                                               cacgGCGTCACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAgtgt                      <  2:1436748/100‑1 (MQ=255)
                                                                               cacgGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGtgtt                     <  1:263925/100‑1 (MQ=255)
                                                                                acgGCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGtgttg                    <  2:1706766/100‑1 (MQ=255)
                                                                                   gCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGttgttg                 <  2:1880008/100‑1 (MQ=255)
                                                                                   gCGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGttgttg                 <  1:644589/100‑1 (MQ=255)
                                                                                    cGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTgg                <  1:1689272/100‑1 (MQ=255)
                                                                                    cGTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTgg                <  1:1035049/100‑1 (MQ=255)
                                                                                     gTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGt               <  1:4431/100‑1 (MQ=255)
                                                                                     gTCACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGt               <  1:1332448/100‑1 (MQ=255)
                                                                                       cACTTTTACTGGCCCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGta              <  2:1034711/100‑1 (MQ=255)
                                                                                       cACTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaa             >  2:1062639/1‑100 (MQ=255)
                                                                                        aCTTGTACTGG‑CCCCACGGTCGTTAAATGGCCCGTAGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaat            <  1:1462646/100‑1 (MQ=255)
                                                                                         cTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaata           >  1:1766153/1‑100 (MQ=255)
                                                                                         cTTTTACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaata           >  1:456319/1‑100 (MQ=255)
                                                                                          ttttACTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGtaataa          <  1:313473/100‑1 (MQ=255)
                                                                                              aCTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCgg      >  2:1590530/1‑100 (MQ=255)
                                                                                              aCTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCgg      <  2:1073221/100‑1 (MQ=255)
                                                                                              aCTGG‑CCCCACGGTCGTTAAATGTCCCGTGGCGATGATGGGCAGAGGCTGATCGCCGCGCAGTTTGCAGGCATCGGCATAGTGTTGTTGGTAATAATCgg      <  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: ‑