Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,895,757 T→A 7.8% L56Q (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,7570TA 7.8% 195.4 / 7.4 64L56Q (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 = 9.09e-01

AATACGCGAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCT  >  NC_002947/4895608‑4895904
                                                                                                                                                     |                                                                                                                                                   
aatACGCGAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCct                                                                                                                                                     >  2:11732/1‑150 (MQ=255)
   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:310164/1‑150 (MQ=255)
                                                                                                    aCGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTc                                                 >  2:1136724/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                                                           >  2:191083/1‑122 (MQ=255)
                                                                                                                      aTCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGc                                                           <  1:191083/122‑1 (MQ=255)
                                                                                                                      aTCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGCTTTa                               >  2:347317/1‑150 (MQ=255)
                                                                                                                      aTCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGCGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTa                               <  2:1297200/150‑1 (MQ=255)
                                                                                                                          aTGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTgg                           <  1:1687089/150‑1 (MQ=255)
                                                                                                                           tGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTAt                                                             >  1:1924705/1‑115 (MQ=255)
                                                                                                                           tGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTAt                                                             <  2:1924705/115‑1 (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                                  >  1:1778201/3‑124 (MQ=255)
                                                                                                                                             taccGGCCAGCAGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGAt                                  <  2:1778201/122‑1 (MQ=255)
                                                                                                                                              gCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTc       >  1:456538/1‑150 (MQ=255)
                                                                                                                                                cGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTcgc     >  2:1084941/1‑150 (MQ=255)
                                                                                                                                                   ccTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCt  <  1:1136724/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
AATACGCGAGCGGCTCTTCGTCATCGTCTTCCAGACCGACACCGTAGCCGGTCGCCGCTTCGACAAGATCCTGTTGCTGATCATCCTGGCCAGCCTGGTCACGGTGATCCTCGACAGCATCGATGAGGTGCACCAGGGCTACGCCGGCCTGCTGGCCGGTATCGAGTGGGGCTTTACCGCGATCTTCCTGGCCGAATACATCACCCGCCTGTACTGTTCGCCCAAGCCGCTGCGCTATGCCTTCAGTTTCTACGGCCTGGTCGATTTACTGGCGATCGTGCCGGGGATCCTCGCGCT  >  NC_002947/4895608‑4895904

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