Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,049,182 T→G 11.1% K141Q (AAG→CAG)  PP_0907 ← RND transporter membrane fusion protein‑like protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,049,1820TG11.1% 254.5 / 25.2 99K141Q (AAG→CAG) PP_0907RND transporter membrane fusion protein‑like protein
Reads supporting (aligned to +/- strand):  ref base T (40/48);  new base G (5/6);  total (45/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.64e-01

GTGCATTGATGCGTTTCTTGATCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCTTCCACCTCCA  >  NC_002947/1049035‑1049319
                                                                                                                                                   |                                                                                                                                         
gTGCATTGATGCGTTTCTTGATCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTgg                                                                                                                                         >  2:2197683/1‑150 (MQ=255)
 ggCATTGATGCGGTTCTTGGTCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTAATGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGc                                                                                                                                        <  2:394733/149‑1 (MQ=255)
     ttGATGCGTTTCTTGATCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAgg                                                                                                                                    >  2:952009/1‑150 (MQ=255)
         tGCGTTTCTTGATCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCa                                                                                                                                >  1:1940847/1‑150 (MQ=255)
            gTTTCTTGATCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTGGGCCTGGCCAAGGt                                                                                                                             <  2:1328095/150‑1 (MQ=255)
             tttCTTGATCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTc                                                                                                                            <  1:2075179/150‑1 (MQ=255)
             tttCTTGATCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTc                                                                                                                            >  1:603667/1‑150 (MQ=255)
              ttCTTGATCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCg                                                                                                                           <  1:250106/150‑1 (MQ=255)
                   gATCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAg                                                                                                                      <  1:1057537/150‑1 (MQ=255)
                    aTCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTGGGCCTGGCCAAGGTCgactt                                                                                                                       >  2:2186483/1‑144 (MQ=255)
                    aTCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTGGGCCTGGCCAAGGTCgactt                                                                                                                       <  1:2186483/148‑5 (MQ=255)
                       ttctGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCg                                                                                                                  >  2:2122142/1‑150 (MQ=255)
                       ttctGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGCCTGAGACGCCCGTAGTTGATCCTGGTCAGGCCCAGGGCGGCCAGGGCg                                                                                                                  >  2:2069954/1‑150 (MQ=255)
                          tGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTa                                                                                                               <  1:2069954/150‑1 (MQ=255)
                          tGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTa                                                                                                               >  2:445034/1‑150 (MQ=255)
                          tGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTa                                                                                                               <  1:1925346/150‑1 (MQ=255)
                          tGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCACAGGTCGGCCAGGGCGGTa                                                                                                               >  2:984691/1‑150 (MQ=255)
                           gcgcCTTAAGCTGCTCGAGCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAAGCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTAc                                                                                                              <  2:1422983/150‑1 (MQ=255)
                            cgcCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTAcc                                                                                                             >  1:2284858/1‑150 (MQ=255)
                             gcCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCg                                                                                                            >  1:1206959/1‑150 (MQ=255)
                               cTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAg                                                                                                          >  1:1066165/1‑150 (MQ=255)
                                ttAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGc                                                                                                         >  1:491879/1‑150 (MQ=255)
                                 tAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCa                                                                                                        >  1:1242627/1‑150 (MQ=255)
                                 tAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCa                                                                                                        <  2:1066165/150‑1 (MQ=255)
                                   aGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTGGGCCTGGCCAAGGTc                                                                                                                            <  2:1384843/128‑1 (MQ=255)
                                   aGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTGGGCCTGGCCAAGGTc                                                                                                                            >  1:1384843/1‑128 (MQ=255)
                                    gctgctCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGa                                                                                                           <  2:907841/144‑1 (MQ=255)
                                    gctgctCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGa                                                                                                           >  1:907841/1‑144 (MQ=255)
                                    gctgctCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAggg                                                                                                     <  2:213355/150‑1 (MQ=255)
                                     ctgctCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGc                                                                                                    <  1:1519479/150‑1 (MQ=255)
                                        ctCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGt                                                                                                 >  2:718337/1‑150 (MQ=255)
                                         tCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCCGCTGGCGGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGtt                                                                                                >  1:761986/1‑150 (MQ=255)
                                         tCGACCACCGCCTGGTTGCGACGGTAATGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGtt                                                                                                >  2:2082656/1‑150 (MQ=255)
                                           taccaccGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTc                                                                                              <  1:1715342/149‑1 (MQ=255)
                                           gACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTc                                                                                              >  1:522670/1‑150 (MQ=255)
                                            accaccGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTcc                                                                                             <  1:718337/150‑1 (MQ=255)
                                             ccaccGCTGGGTTGCGCCGGTACTGGGAGGAAAGGCTGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGTCCGAAGTCTACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCt                                                                                            <  1:112767/150‑1 (MQ=255)
                                             ccaccGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCACGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCt                                                                                            >  1:1036757/1‑150 (MQ=255)
                                               accGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGa                                                                                          <  2:603667/150‑1 (MQ=255)
                                                ccGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAggg                                                                                                     <  1:1765235/138‑1 (MQ=255)
                                                ccGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAggg                                                                                                     >  2:1765235/1‑138 (MQ=255)
                                                 cGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAggg                                                                                                                    <  1:1959897/122‑1 (MQ=255)
                                                 cGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAggg                                                                                                                    >  2:1959897/1‑122 (MQ=255)
                                                 cGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGAtc                                                                                        <  1:445034/150‑1 (MQ=255)
                                                       gTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGtt                                                                                  <  2:249041/150‑1 (MQ=255)
                                                                gTAGTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCaa                                                                         <  1:748207/150‑1 (MQ=255)
                                                                  aCTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAaga                                                                       <  2:1206959/150‑1 (MQ=255)
                                                                    tGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGc                                                                                                                   <  1:1572225/104‑1 (MQ=255)
                                                                    tGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGc                                                                                                                   >  2:1572225/1‑104 (MQ=255)
                                                                     gggAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCa                                                                    <  1:1088275/150‑1 (MQ=255)
                                                                     gggAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCa                                                                    <  1:2082656/150‑1 (MQ=255)
                                                                      ggaggaAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGa                                                                                          <  2:1198875/127‑1 (MQ=255)
                                                                      ggaggaAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGa                                                                                          >  1:1198875/1‑127 (MQ=255)
                                                                      ggaggaAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAagagacg                                                                   <  2:1242627/150‑4 (MQ=255)
                                                                      ggaggaAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAg                                                                   >  1:587408/1‑150 (MQ=255)
                                                                       gaggaAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGt                                                                  <  1:306790/150‑1 (MQ=255)
                                                                                cGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGCGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGcc                                                         <  1:1315350/150‑1 (MQ=255)
                                                                                  gTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGATGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCtt                                                       <  2:2284858/150‑1 (MQ=255)
                                                                                    cGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGAATGACGTCCACGAGGTGGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGa                                                     <  2:674519/150‑1 (MQ=255)
                                                                                        cTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTg                                                 <  2:2260321/150‑1 (MQ=255)
                                                                                                 ggAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGa                                        <  1:952009/150‑1 (MQ=255)
                                                                                                  gatctGGCCGAATCGCCCACCAGCTAAGGTCCTCGGCAGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCTTTTAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAg                                       <  1:1258/146‑1 (MQ=255)
                                                                                                  gagaTGGCCGAATCGCCCACCAGCTGGCTGCCGAGGGTGATGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAg                                       <  2:1999624/150‑1 (MQ=255)
                                                                                                       ggCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGaa                                  <  2:1940847/150‑1 (MQ=255)
                                                                                                        gCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGa                                                                                                           <  2:69496/76‑1 (MQ=255)
                                                                                                        gCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGa                                                                                                           >  1:69496/1‑76 (MQ=255)
                                                                                                        gCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAg                                 >  2:733220/1‑150 (MQ=255)
                                                                                                              tCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCCGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAgg                           <  1:1728798/150‑1 (MQ=255)
                                                                                                              tCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAgg                           >  1:1984819/1‑150 (MQ=255)
                                                                                                                gCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTg                                                                <  1:167779/111‑1 (MQ=255)
                                                                                                                gCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTg                                                                >  2:167779/1‑111 (MQ=255)
                                                                                                                 cccACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGc                                                                           >  2:1166316/1‑99 (MQ=255)
                                                                                                                 cccACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGc                                                                           <  1:1166316/99‑1 (MQ=255)
                                                                                                                 cccACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGAt                        <  1:2122142/150‑1 (MQ=255)
                                                                                                                 cccACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCCGTTGGCCTGAGTCGAAGTGCAGGGAt                        >  2:107912/1‑150 (MQ=255)
                                                                                                                  ccaccaGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGAtt                       <  1:1287106/150‑1 (MQ=255)
                                                                                                                  ccaccaGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGAtt                       >  2:552720/1‑150 (MQ=255)
                                                                                                                   caccaGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGAttt                      <  1:2197683/150‑1 (MQ=255)
                                                                                                                   caccaGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGAttt                      >  1:1562956/1‑150 (MQ=255)
                                                                                                                   caccaGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGAttt                      <  2:1036757/150‑1 (MQ=255)
                                                                                                                    accaGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGt                                <  2:2147455/139‑1 (MQ=255)
                                                                                                                    accaGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGt                                >  1:2147455/1‑139 (MQ=255)
                                                                                                                     ccaGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGa                    <  1:552720/150‑1 (MQ=255)
                                                                                                                     ccaGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGa                    <  1:733220/150‑1 (MQ=255)
                                                                                                                         ctggctgCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGc                                                                     <  2:2401214/97‑1 (MQ=255)
                                                                                                                         ctggctgCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGc                                                                     >  1:2401214/1‑97 (MQ=255)
                                                                                                                           ggctgCCACGGCCCAGGGCGACCTTGGCCAGGCCCAGGTGGGCGCTGGCGGTAGCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGc              >  2:341761/1‑150 (MQ=255)
                                                                                                                            gctgCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGcc             >  1:324551/1‑150 (MQ=255)
                                                                                                                                ccACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATggtggtc         >  2:1571264/1‑144 (MQ=255)
                                                                                                                                 cACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCTTCCa        >  2:742617/1‑150 (MQ=255)
                                                                                                                                  aCGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGtt                                                                 <  2:1998169/92‑1 (MQ=255)
                                                                                                                                  aCGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGtt                                                                 >  1:1998169/1‑92 (MQ=255)
                                                                                                                                       cGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCTTCCACCTCCa  >  2:1365486/1‑150 (MQ=255)
                                                                                                                                             cGACCTGGGCCTGGCCAAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCt                                          <  2:713334/104‑1 (MQ=255)
                                                                                                                                             cGACCTGGGCCTGGCCAAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCt                                          >  1:713334/1‑104 (MQ=255)
                                                                                                                                             cGACCTGGGCCTGGCCAAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGa                                        >  1:187613/1‑106 (MQ=255)
                                                                                                                                             cGACCTGGGCCTGGCCAAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGa                                        <  2:187613/106‑1 (MQ=255)
                                                                                                                                               aCCTGGGCCTGGCCAAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGa                                        >  1:1630680/1‑104 (MQ=255)
                                                                                                                                               aCCTGGGCCTGGCCAAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGa                                        <  2:1630680/104‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
GTGCATTGATGCGTTTCTTGATCTTCTGCGCCTTAAGCTGCTCGACCACCGCCTGGTTGCGCCGGTACTGGGAGGAAAGGCGGTCGAACTCGCCGCGGGAGATGGCCGAATCGCCCACCAGCTGGCTGCCACGGCCGAAGTCGACCTTGGCCAGGCCCAGGTCGGCCAGGGCGGTACCGAGCAGGGCTGTTTCCTGATCGTCGTTGAGCTGCAAGAGCAGTTGCCCGGCCTTGACCTGTTGGCCTGAGTCGAAGTGCAGGGATTTGACGATGCCTTCCACCTCCA  >  NC_002947/1049035‑1049319

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: