Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,248,403 C→A 7.5% Q556H (CAG→CAT PP_3724 ← acyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,248,4030CA 7.5% 217.0 / 3.0 69Q556H (CAG→CATPP_3724acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (34/29);  new base A (3/3);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGCACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGC  >  NC_002947/4248258‑4248550
                                                                                                                                                 |                                                                                                                                                   
cgGCACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCa                                                                                                                                                 >  1:111845/1‑150 (MQ=255)
   cACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCa                                                                                                                                              <  1:34263/150‑1 (MQ=255)
     ccGATCACCAAGGCCAGGGGAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGCGAGCGCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGTGCAGATGCGCCAGt                                                                                                                                            <  1:448930/150‑1 (MQ=255)
       gATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGTGCAGATGCGCCAGt                                                                                                                                            <  2:292665/148‑1 (MQ=255)
       gATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGTGCAGATGCGCCAGt                                                                                                                                            >  1:292665/1‑148 (MQ=255)
       gATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGc                                                                                                                                          >  2:1250538/1‑150 (MQ=255)
        aTCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCGCCAGTGCa                                                                                                                                         >  2:1355/1‑150 (MQ=255)
         tCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGTGCAGATGCGCCAGTGCat                                                                                                                                        <  2:1421526/150‑2 (MQ=255)
          cACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGc                                                                                                                                       <  1:1044891/150‑1 (MQ=255)
             cAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGTGCAGATGCGCCAGTGCatcgct                                                                                                                                    >  1:235237/1‑145 (MQ=255)
             cAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTc                                                                                                                                    >  1:1284490/1‑150 (MQ=255)
              aaGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTcc                                                                                                                                   >  2:73323/1‑150 (MQ=255)
               aGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCGCTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCa                                                                                                                                  <  2:247459/150‑1 (MQ=255)
               aGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCa                                                                                                                                  <  1:73323/150‑1 (MQ=255)
                       gCAAATTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCCGGCCCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCgg                                                                                                                          <  1:93826/150‑1 (MQ=255)
                         aaaGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAg                                                                                                                                             >  1:1080133/1‑129 (MQ=255)
                         aaaGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAg                                                                                                                                             <  2:1080133/129‑1 (MQ=255)
                                  ttCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAg                                                                                                               <  2:362478/150‑1 (MQ=255)
                                   tCACGCTTGATGCCCACTATGAGCTTGGCGAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCATGGCGCGGGCATCGCGCAgg                                                                                                              >  1:397667/1‑150 (MQ=255)
                                   tCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAgg                                                                                                              <  1:492027/150‑1 (MQ=255)
                                       gCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAAACTGCGGGGCTGAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCa                                                                                                          >  1:824991/1‑150 (MQ=255)
                                           gATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGc                                                                                                                                                  >  1:1155376/1‑106 (MQ=255)
                                           gATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGc                                                                                                                                                  <  2:1155376/106‑1 (MQ=255)
                                             tGCCCACTATGAGGTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTGCAGGGCGCGGGCATCGCGCAGGTTCAACGCca                                                                                                    <  2:932387/150‑1 (MQ=255)
                                              gCCCACTATGAGCTCGGCCAAGGGCTGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCcac                                                                                                   >  1:897138/1‑150 (MQ=255)
                                               cccACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCcaca                                                                                                  <  1:238824/150‑1 (MQ=255)
                                                  aCTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCCGCACCTGCTCAAACGCAACCTCCGGGGCGTAGACAACGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGcc                                                                                               >  2:719542/1‑149 (MQ=255)
                                                          cTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGTTGGCGGAAAACCGGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTc                                                                                       <  2:397667/150‑1 (MQ=255)
                                                           tCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATc                                                                                                                     >  2:520120/1‑119 (MQ=255)
                                                           tCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATc                                                                                                                     <  1:520120/119‑1 (MQ=255)
                                                           tCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCg                                                                                      >  2:1251947/1‑150 (MQ=255)
                                                           tCGGCCAAGGGCGGCTCGGCCATGAGCTCGAGCAGCACCTGATCAAACGCCACCTGCGGGGCCTAGTCAGCGGTGCACTGGCGCCTCTGCACCAGTGCAGCTTCCAGCGCGCGGGCATCGCGAAGGTTCAACGCCACAGAGCCAGCTTCg                                                                                      >  2:584500/1‑150 (MQ=255)
                                                                    ggcggcTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGc                                                                             <  1:900695/150‑1 (MQ=255)
                                                                     gcggcTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCg                                                                            <  1:1250538/150‑1 (MQ=255)
                                                                           cGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCt                                                                      >  1:1151360/1‑150 (MQ=255)
                                                                             gCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTgc                                                                    >  2:1331051/1‑150 (MQ=255)
                                                                                aTGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAg                                                                                                               >  2:1383278/1‑104 (MQ=255)
                                                                                aTGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAg                                                                                                               <  1:1383278/104‑1 (MQ=255)
                                                                                aTGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTgcgcg                                                                 >  1:1000315/1‑150 (MQ=255)
                                                                                 tGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCAc                                                                                              >  2:851631/1‑120 (MQ=255)
                                                                                 tGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCAc                                                                                              <  1:851631/120‑1 (MQ=255)
                                                                                      tCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGAcc                                                           >  2:1164087/1‑150 (MQ=255)
                                                                                                 cTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTAc                                                <  2:897138/150‑1 (MQ=255)
                                                                                                     tCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCtt                                                                                 >  2:922308/1‑113 (MQ=255)
                                                                                                     tCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCtt                                                                                 <  1:922308/113‑1 (MQ=255)
                                                                                                     tCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAg                                            >  1:831457/1‑150 (MQ=255)
                                                                                                      cAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAgg                                           >  1:446985/1‑150 (MQ=255)
                                                                                                       aaaCGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAggg                                          >  1:601091/1‑150 (MQ=255)
                                                                                                       aaaCGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAggg                                          >  1:532819/1‑150 (MQ=255)
                                                                                                             aaCCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAAcc                                    >  2:1313637/1‑150 (MQ=255)
                                                                                                             aaCCTAAGGGGCGGAGTCACAGCTGCACTACTTCACATGCACCGGGGCAGCGGCTAGACCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGGAAcc                                    >  2:945338/1‑150 (MQ=255)
                                                                                                                  gCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGAcgc                               <  2:831457/150‑1 (MQ=255)
                                                                                                                       gCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACgcgggcg                          <  2:855587/150‑1 (MQ=255)
                                                                                                                             tCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCt                    <  1:1251947/150‑1 (MQ=255)
                                                                                                                             tCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCt                    >  2:1173360/1‑150 (MQ=255)
                                                                                                                             tCAGCGATGCACTGGCGCATCTGCAACAGTGCCGATTCCAGGTCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCt                    <  2:824991/150‑1 (MQ=255)
                                                                                                                               aGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGc                  <  2:1151360/150‑1 (MQ=255)
                                                                                                                                  gATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTgcgc                                                                  <  2:1369380/99‑1 (MQ=255)
                                                                                                                                  gATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTgcgc                                                                  >  1:1369380/1‑99 (MQ=255)
                                                                                                                                       aCTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGgcgc          >  1:549553/1‑150 (MQ=255)
                                                                                                                                        cTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAggg                                          >  2:901005/1‑117 (MQ=255)
                                                                                                                                        cTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAggg                                          <  1:901005/117‑1 (MQ=255)
                                                                                                                                         tGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGt                                        >  1:146187/1‑118 (MQ=255)
                                                                                                                                         tGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGt                                        <  2:146187/118‑1 (MQ=255)
                                                                                                                                             gcATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCgg    >  2:472314/1‑150 (MQ=255)
                                                                                                                                             gcATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCgg    <  1:1331051/150‑1 (MQ=255)
                                                                                                                                             gcATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCgg    >  2:1099095/1‑150 (MQ=255)
                                                                                                                                             gcATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCgg    >  2:1099007/1‑150 (MQ=255)
                                                                                                                                               aTCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGc  <  1:1355/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGGCACCGATCACCAAGGCCAGGGCAAAGTTGTTTTCACGCTTGATGCCCACTATGAGCTCGGCCAAGGGCGGCTCGGCCATGGGCTCGAGCAGCACCTGATCAAACGCAACCTGCGGGGCGTAGTCAGCGATGCACTGGCGCATCTGCACCAGTGCAGCTTCCAGGGCGCGGGCATCGCGCAGGTTCAACGCCACAGCACCAGCTTCGGTCTTGTGCGGCAGCTGCGCGCTGACCGCCTTCAGTACCAAGGGGTAACCGACGCGGGCGGCGTCTGCCACGGCGCGTTCGGGC  >  NC_002947/4248258‑4248550

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: