Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,638,924 T→A 5.3% W532R (TGG→AGG)  PP_4103 → low‑affinity inorganic phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,638,9240TA 5.3% 238.0 / 3.7 75W532R (TGG→AGG) PP_4103low‑affinity inorganic phosphate transporter
Reads supporting (aligned to +/- strand):  ref base T (36/35);  new base A (2/2);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.81e-01

TGTCGACCACGCATGTGCTGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCGGGTGAACCCGCTCCTACGGGC  >  NC_002947/4638778‑4639073
                                                                                                                                                  |                                                                                                                                                      
tgtCGACCACGCATGTGCTGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGgc                                                                                                                                                     <  2:1419669/150‑1 (MQ=255)
   cGACCACGCATGTGCTGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGgctgc                                                                                                                                                  <  2:84122/150‑1 (MQ=255)
    gACCACGCATGTGCTGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGgctgct                                                                                                                                                 >  2:277363/1‑150 (MQ=255)
     aCCACGCATGTGCTGAGTTCGGGGGTGGCCGGCGCCAAGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCCGCGTCGATTATGTTTGCTTGGAGCCTGTACTGGctgctg                                                                                                                                                >  2:185551/1‑150 (MQ=255)
      ccACGCATGTGCTGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCGATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCAGTACAGgctgccag                                                                                                                                               <  1:1215002/150‑4 (MQ=255)
       cACGCATGTGCTGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGaa                                                                                                                                             >  1:1205186/1‑150 (MQ=255)
        aCGCATGTGCTGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAAc                                                                                                                                            <  1:236769/150‑1 (MQ=255)
            aTGTGCTGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACcaca                                                                                                                                        <  1:369027/150‑1 (MQ=255)
               tgCTGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACAtct                                                                                                                                     >  1:1110446/1‑150 (MQ=255)
                 cTGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACAtcttc                                                                                                                                   >  2:39221/1‑150 (MQ=255)
                 cTGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACAtcttc                                                                                                                                   <  2:1025091/150‑1 (MQ=255)
                  tGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACAtcttct                                                                                                                                  >  1:1289986/1‑150 (MQ=255)
                  tGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACAtcttct                                                                                                                                  <  2:862764/150‑1 (MQ=255)
                   gAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACcacagatcgga                                                                                                                                 >  2:536372/1‑143 (MQ=255)
                      ttCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGgc                                                                                                                                                     <  1:811031/128‑1 (MQ=255)
                      ttCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGgc                                                                                                                                                     >  2:811031/1‑128 (MQ=255)
                        cGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAAccc                                                                                                                            <  1:640737/150‑1 (MQ=255)
                         gggggTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGATGCTGAACCACATCTTCTGAACCCg                                                                                                                           >  1:664931/1‑150 (MQ=255)
                                    gcgcCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCtt                                                                                                                <  1:2557/150‑1 (MQ=255)
                                        cATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGt                                                                                                            >  1:282078/1‑150 (MQ=255)
                                           ggtggCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTCgg                                                                                                         >  1:765830/1‑150 (MQ=255)
                                                ccAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGCAGCGAAGTCACCTAACCCGTACGCCGCCTTCTGTAGGAGCgg                                                                                                    >  2:825294/1‑150 (MQ=255)
                                                 cAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCggg                                                                                                   >  1:973088/1‑150 (MQ=255)
                                                 cAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCggg                                                                                                   >  2:909941/1‑150 (MQ=255)
                                                 cAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCggg                                                                                                   <  2:1135498/150‑1 (MQ=255)
                                                 cAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGAGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCggg                                                                                                   >  1:929783/1‑150 (MQ=255)
                                                 cAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCATATGGCCGGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCggg                                                                                                   >  1:685243/1‑150 (MQ=255)
                                                      ggTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTAc                                                                                              <  1:277363/150‑1 (MQ=255)
                                                       gTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCCGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTGCTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTAcc                                                                                             >  1:526546/1‑150 (MQ=255)
                                                        tCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTAccc                                                                                            >  2:1167110/1‑150 (MQ=255)
                                                        tCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTAccc                                                                                            >  2:204060/1‑150 (MQ=255)
                                                        tCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTAccc                                                                                            <  1:204060/150‑1 (MQ=255)
                                                        tCGGGGATTCAGAAACGCACCTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTAccc                                                                                            <  2:685243/150‑1 (MQ=255)
                                                          ggggATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCcgc                                                                                          <  1:909941/150‑1 (MQ=255)
                                                          ggggATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCcgc                                                                                          >  2:1005165/1‑150 (MQ=255)
                                                             gATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGaa                                                                                       >  2:575017/1‑150 (MQ=255)
                                                                  aGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATa                                                                                     >  2:817802/1‑147 (MQ=255)
                                                                  aGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATa                                                                                     <  1:817802/147‑1 (MQ=255)
                                                                  aGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTACTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACgg                                                                                  >  2:985798/1‑150 (MQ=255)
                                                                        gCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAgcgg                                                                            <  1:1326041/150‑1 (MQ=255)
                                                                        gCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAgcgg                                                                            >  2:426190/1‑150 (MQ=255)
                                                                           cTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGcggcgg                                                                         <  2:1110446/150‑1 (MQ=255)
                                                                           cTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGcggcgg                                                                         >  2:1386796/1‑150 (MQ=255)
                                                                                          tgATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAAcggcg                                                          <  1:39221/150‑1 (MQ=255)
                                                                                            aTGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCggg                                                        >  1:236333/1‑150 (MQ=255)
                                                                                                ccTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCCGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTggc                                                    <  2:973088/150‑2 (MQ=255)
                                                                                                 cTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTa                                                                                                           <  2:138813/94‑1 (MQ=255)
                                                                                                 cTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTa                                                                                                           >  1:138813/1‑94 (MQ=255)
                                                                                                 cTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGgag                                                   >  2:1200246/1‑150 (MQ=255)
                                                                                                 cTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGCACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGgag                                                   >  2:1136770/1‑150 (MQ=255)
                                                                                                     gTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTg                                               <  2:90079/150‑1 (MQ=255)
                                                                                                     gTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTg                                               <  2:20748/150‑1 (MQ=255)
                                                                                                      tGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTgg                                              >  1:299398/1‑150 (MQ=255)
                                                                                                             ccTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTgg                                              >  2:952210/1‑143 (MQ=255)
                                                                                                             ccTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTgg                                              <  1:952210/143‑1 (MQ=255)
                                                                                                             ccTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCcagcag                                       >  1:501328/1‑150 (MQ=255)
                                                                                                                gCCAGCGGCAATGCTGCTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCTAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAggcc                                    <  1:432789/150‑1 (MQ=255)
                                                                                                                  cAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAggcccg                                  <  1:619225/150‑1 (MQ=255)
                                                                                                                     cGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAggcccggcc                               <  1:185551/150‑1 (MQ=255)
                                                                                                                          aTGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAAc                                                              >  2:921064/1‑114 (MQ=255)
                                                                                                                          aTGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAAc                                                              <  1:921064/114‑1 (MQ=255)
                                                                                                                             cTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACgg                                                                   >  2:1400549/1‑106 (MQ=255)
                                                                                                                             cTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACgg                                                                   <  1:1400549/106‑1 (MQ=255)
                                                                                                                              tGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCgg                      <  2:299398/150‑1 (MQ=255)
                                                                                                                              tGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCgg                      <  1:1200246/150‑1 (MQ=255)
                                                                                                                              tGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCgg                      <  2:1289986/150‑1 (MQ=255)
                                                                                                                                    tcagcAGCCAGTACAGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGgaga                                                  <  2:660801/113‑1 (MQ=255)
                                                                                                                                    tcagcAGCCAGTACAGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGgaga                                                  >  1:660801/4‑116 (MQ=255)
                                                                                                                                    ctggcAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCGGGTGAAc                <  1:575017/150‑1 (MQ=255)
                                                                                                                                     tggcAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAggcccgg                                 >  1:297259/1‑132 (MQ=255)
                                                                                                                                     tggcAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAggcccgg                                 <  2:297259/132‑1 (MQ=255)
                                                                                                                                     tggcAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCGGGTGAAcc               >  1:687871/1‑150 (MQ=255)
                                                                                                                                     tctcAGCCTGTACTGGCTGCTTATCCACATCTTCTCAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCGGGTGAAcc               <  2:664931/147‑1 (MQ=255)
                                                                                                                                       gcAGCCAGTACAGGCTGCTGAACCCCATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCGGGTGAACCCg             >  1:1260797/1‑150 (MQ=255)
                                                                                                                                                 cTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCCAATACTGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCGCGGCCCTTTCGCGGGTGAACCCGCTCCTACggg   <  1:114599/150‑1 (MQ=255)
                                                                                                                                                  tGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCcagca                                        >  1:337520/1‑112 (MQ=255)
                                                                                                                                                  tGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCcagca                                        <  2:337520/112‑1 (MQ=255)
                                                                                                                                                  tGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCGGGTGAACCCGCTCCTACGGGc  <  2:198543/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
TGTCGACCACGCATGTGCTGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCGGGTGAACCCGCTCCTACGGGC  >  NC_002947/4638778‑4639073

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: