Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,638,919 T→A 5.3% L530Q (CTG→CAG)  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,9190TA 5.3% 243.8 / 3.7 76L530Q (CTG→CAG) PP_4103low‑affinity inorganic phosphate transporter
Reads supporting (aligned to +/- strand):  ref base T (38/34);  new base A (2/2);  total (40/36)
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.88e-01

CTGCCGGTGTCGACCACGCATGTGCTGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCGGGTGAACCCG  >  NC_002947/4638771‑4639062
                                                                                                                                                     |                                                                                                                                                
cTGCCGGTGTCGACCACGCATGTGCTGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGAATCACACTCCTGCGGGCGGCAATGCTGTTGGCTGGCAGCCTg                                                                                                                                                 >  2:114599/1‑150 (MQ=255)
   ccGGTGTCGACCACGCATGTGCTGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTAc                                                                                                                                              >  1:1135498/1‑150 (MQ=255)
       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                                                                                                                                          >  2:811031/1‑128 (MQ=255)
                             ttCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGgc                                                                                                                                          <  1:811031/128‑1 (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                                                                                        >  2:909941/1‑150 (MQ=255)
                                                        cAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCggg                                                                                        <  2:1135498/150‑1 (MQ=255)
                                                        cAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCggg                                                                                        >  1:973088/1‑150 (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                                                                                 <  1:204060/150‑1 (MQ=255)
                                                               tCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTAccc                                                                                 >  2:204060/1‑150 (MQ=255)
                                                               tCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTAccc                                                                                 >  2:1167110/1‑150 (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                                                                          <  1:817802/147‑1 (MQ=255)
                                                                         aGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATa                                                                          >  2:817802/1‑147 (MQ=255)
                                                                         aGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTACTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACgg                                                                       >  2:985798/1‑150 (MQ=255)
                                                                               gCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAgcgg                                                                 >  2:426190/1‑150 (MQ=255)
                                                                               gCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAgcgg                                                                 <  1:1326041/150‑1 (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:20748/150‑1 (MQ=255)
                                                                                                            gTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTg                                    <  2:90079/150‑1 (MQ=255)
                                                                                                             tGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTgg                                   >  1:299398/1‑150 (MQ=255)
                                                                                                                     ccTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTgg                                   <  1:952210/143‑1 (MQ=255)
                                                                                                                     ccTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTgg                                   >  2:952210/1‑143 (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                                                   <  1:921064/114‑1 (MQ=255)
                                                                                                                                  aTGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAAc                                                   >  2:921064/1‑114 (MQ=255)
                                                                                                                                     cTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACgg                                                        <  1:1400549/106‑1 (MQ=255)
                                                                                                                                     cTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACgg                                                        >  2:1400549/1‑106 (MQ=255)
                                                                                                                                      tGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCgg           <  2:299398/150‑1 (MQ=255)
                                                                                                                                      tGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCgg           <  2:1289986/150‑1 (MQ=255)
                                                                                                                                      tGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCgg           <  1:1200246/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                      <  2:297259/132‑1 (MQ=255)
                                                                                                                                             tggcAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAggcccgg                      >  1:297259/1‑132 (MQ=255)
                                                                                                                                             tggcAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCGGGTGAAcc    >  1:687871/1‑150 (MQ=255)
                                                                                                                                             tctcAGCCTGTACTGGCTGCTTATCCACATCTTCTCAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCGGGTGAAcc    <  2:664931/147‑1 (MQ=255)
                                                                                                                                               gcAGCCAGTACAGGCTGCTGAACCCCATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCGGGTGAACCCg  >  1:1260797/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
CTGCCGGTGTCGACCACGCATGTGCTGAGTTCGGGGGTGGCCGGCGCCATGGTGGCCAATGGGTCGGGGATTCAGAAACGCACTTTGTTCAATCTGCTGATGGCCTGGGTGCTGACCCTGCCAGCGGCAATGCTGTTGGCTGGCAGCCTGTACTGGCTGCTGAACCACATCTTCTGAACCCGTACGCCGCCTTCTGTAGGAGCGGGTTTACCCGCGAATACGGTAGCGGCGGCAACGGTGAACGGCGGGTGGAGATTGGCCAGCAGGCCCGGCCCTTTCGCGGGTGAACCCG  >  NC_002947/4638771‑4639062

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: