Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,895,760 T→A 7.8% L57Q (CTG→CAG)  PP_4304 → VIC family cation transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,895,7600TA 7.8% 192.6 / 7.7 64L57Q (CTG→CAG) PP_4304VIC family cation transporter
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base A (3/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.08e-01

ACGCGAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTA  >  NC_002947/4895611‑4895907
                                                                                                                                                     |                                                                                                                                                   
aCGCGAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCctgct                                                                                                                                                     >  2:1329321/1‑150 (MQ=255)
  gcgAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTgg                                                                                                                                                   >  1:1297200/1‑150 (MQ=255)
   cgAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGc                                                                                                                                                  <  1:1011462/150‑1 (MQ=255)
        ggCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTa                                                                                                                                             <  2:333903/150‑1 (MQ=255)
            ctTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGa                                                                                                                                         <  1:1360978/150‑1 (MQ=255)
               cGTCATCGTCTTCCAGACCGACCCCGTAGCCGGTCGACTCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCAGGCCTGCTGGCCGGTATCGAGTg                                                                                                                                      <  1:892020/150‑1 (MQ=255)
               cGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACACCGGCCAGCAGGCCggcgta                                                                                                                                           >  2:1782412/1‑141 (MQ=255)
               cGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACACCGGCCAGCAGGCCggcgta                                                                                                                                           <  1:1782412/145‑5 (MQ=255)
                 tCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTggg                                                                                                                                    <  1:1326627/150‑1 (MQ=255)
                         ttCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTAcc                                                                                                                            <  2:1225171/150‑1 (MQ=255)
                             aGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGa                                                                                                                        <  1:11732/150‑1 (MQ=255)
                                                    cgcTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCAcc                                                                                                 >  1:1104128/1‑150 (MQ=255)
                                                     gcTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGc                                                                                                                                  >  2:1233989/1‑116 (MQ=255)
                                                     gcTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGc                                                                                                                                  <  1:1233989/116‑1 (MQ=255)
                                                     gcTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCcccc                                                                                                >  2:949791/1‑150 (MQ=255)
                                                       ttCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGc                                                                                              <  2:1651234/150‑1 (MQ=255)
                                                         cGACAAGAGCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCt                                                                                            <  2:170604/150‑1 (MQ=255)
                                                              aGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCt                                                                                                                     >  1:1578789/1‑120 (MQ=255)
                                                              aGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCt                                                                                                                     <  2:1578789/120‑1 (MQ=255)
                                                               gATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTg                                                                                      <  1:1390365/150‑1 (MQ=255)
                                                                         cTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGaa                                                                                                          <  1:1001836/120‑1 (MQ=255)
                                                                         cTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGaa                                                                                                          >  2:1001836/1‑120 (MQ=255)
                                                                         cTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTc                                                                                   >  1:1154593/1‑143 (MQ=255)
                                                                         cTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTc                                                                                   <  2:1154593/143‑1 (MQ=255)
                                                                          tGATAATCCTGGCCAACCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGc                                                                           >  2:73377/1‑150 (MQ=255)
                                                                           gATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGCTATCGAGTGGCGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGcc                                                                          >  2:294694/1‑150 (MQ=255)
                                                                            atcatcCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCg                                                                         >  2:1765388/1‑150 (MQ=255)
                                                                              catcCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCt                                                                       >  1:43216/1‑150 (MQ=255)
                                                                               atcCTGGCCAGCTTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAg                                                                                                                                        <  1:1592543/84‑1 (MQ=255)
                                                                               atcCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAg                                                                                                                                        >  2:1592543/1‑84 (MQ=255)
                                                                                cccTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTgc                                                                     <  2:1104128/149‑1 (MQ=255)
                                                                                      ccAGCCTGGTCACGGTGATCATAGACAGCAGCGATAAGTGTCACCAGGGCGACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATg                                                               >  2:1876051/1‑150 (MQ=255)
                                                                                        aGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGcc                                                             <  2:43216/150‑1 (MQ=255)
                                                                                                 aCGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTc                                                    >  2:881404/1‑150 (MQ=255)
                                                                                                 aCGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTc                                                    >  2:1136724/1‑150 (MQ=255)
                                                                                                 aCGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTc                                                    >  2:310164/1‑150 (MQ=255)
                                                                                                      gATCCTCGACAGCATCGATGAGGTGCACCAGAGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCATTACCGCCATCTTCATGGCCTAATACATCAACAGAATGTAGTTTTCGGCCAAACTGTGGCGCTTTGCCTTCATTTTCttctg                                               >  2:702984/1‑146 (MQ=255)
                                                                                                          cTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGTCCTg                                           >  1:85726/1‑150 (MQ=255)
                                                                                                          cTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTg                                           >  2:805826/1‑150 (MQ=255)
                                                                                                           tCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGTCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCa                                                         <  1:1208021/135‑1 (MQ=255)
                                                                                                           tCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCa                                                         >  2:1208021/1‑135 (MQ=255)
                                                                                                            cGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGt                                         <  1:73377/150‑1 (MQ=255)
                                                                                                                   aTCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGc                                                              <  1:191083/122‑1 (MQ=255)
                                                                                                                   aTCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGc                                                              >  2:191083/1‑122 (MQ=255)
                                                                                                                   aTCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGCTTTa                                  >  2:347317/1‑150 (MQ=255)
                                                                                                                   aTCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGCGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTa                                  <  2:1297200/150‑1 (MQ=255)
                                                                                                                       aTGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTgg                              <  1:1687089/150‑1 (MQ=255)
                                                                                                                        tGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTAt                                                                <  2:1924705/115‑1 (MQ=255)
                                                                                                                        tGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTAt                                                                >  1:1924705/1‑115 (MQ=255)
                                                                                                                        tGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGc                             >  1:365646/1‑150 (MQ=255)
                                                                                                                          aGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGa                           <  1:347317/150‑1 (MQ=255)
                                                                                                                          aGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTCGGCGTATACCGCGATCTTCCTGGCCGAATACATCGCCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGa                           <  1:1181132/150‑1 (MQ=255)
                                                                                                                                 ccAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTATGGGCGATCGTGcc                    >  2:1337087/1‑150 (MQ=255)
                                                                                                                                 ccAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGcc                    >  1:1367597/1‑150 (MQ=255)
                                                                                                                                  cAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCg                   >  1:452968/1‑150 (MQ=255)
                                                                                                                                   aGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCgg                  <  2:365646/150‑1 (MQ=255)
                                                                                                                                    gggCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCggg                 <  2:1367597/150‑1 (MQ=255)
                                                                                                                                      gCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGa               >  2:207214/1‑150 (MQ=255)
                                                                                                                                        tACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATc             >  2:1286122/1‑150 (MQ=255)
                                                                                                                                        gataccGGCCAGCAGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATc             >  1:1053715/5‑150 (MQ=255)
                                                                                                                                          taccGGCCAGCAGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGAt                                     <  2:1778201/122‑1 (MQ=255)
                                                                                                                                          taccGGCCAGCAGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGAt                                     >  1:1778201/3‑124 (MQ=255)
                                                                                                                                           gCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTc          >  1:456538/1‑150 (MQ=255)
                                                                                                                                             cGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTcgc        >  2:1084941/1‑150 (MQ=255)
                                                                                                                                                ccTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCt     <  1:1136724/150‑1 (MQ=255)
                                                                                                                                                   gctgGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGta  <  1:1765388/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
ACGCGAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCTGTA  >  NC_002947/4895611‑4895907

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