Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,057,001:1 +A 7.1% coding (390/990 nt) PP_2669 ← outer membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,057,0011.A 7.1% 322.2 / 19.0 84coding (390/990 nt)PP_2669outer membrane protein
Reads supporting (aligned to +/- strand):  ref base . (40/38);  new base A (2/4);  total (42/42)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

CGGCCAGGGTCTGGGTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTCC  >  NC_002947/3056853‑3057148
                                                                                                                                                     |                                                                                                                                                   
cGGCCAGGGTCTGGGTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGGCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGgtga                                                                                                                                                     <  1:749552/150‑2 (MQ=255)
  gCCAGGGTCTGGGTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGCGAACCTTGTCGGTGACGGTG‑TCg                                                                                                                                                  >  1:585835/1‑150 (MQ=255)
       ggTCTGGGTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCAc                                                                                                                                             >  1:841839/1‑150 (MQ=255)
        gTCTGGGTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGGCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCAcc                                                                                                                                            <  1:596562/150‑1 (MQ=255)
            gggTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTGGATGTCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTGATCGA‑CACCGTCa                                                                                                                                        <  2:833990/150‑1 (MQ=255)
             ggTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTGATCGA‑CACCGTCAc                                                                                                                                       <  2:1758211/150‑1 (MQ=255)
               tGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGACCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCa                                                                                                                                     >  1:1900540/1‑150 (MQ=255)
                  tGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAgcg                                                                                                                                  >  2:1566332/1‑150 (MQ=255)
                    gtgtCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCa                                                                                                                                >  2:1355509/1‑150 (MQ=255)
                    gtgtCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTGCGATCTGCCCGAGCACCTTGTCGGTGATGGTG‑TCGAGCACCGTCACCAGCGCa                                                                                                                                >  2:907443/1‑150 (MQ=255)
                     tgtCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCAt                                                                                                                               <  2:269686/150‑1 (MQ=255)
                         gATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGCCCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTc                                                                                                                           <  2:137861/150‑1 (MQ=255)
                         gATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCCACAGCGCATCGTc                                                                                                                           >  1:1433457/1‑150 (MQ=255)
                         gATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTc                                                                                                                           >  1:1921104/1‑150 (MQ=255)
                            ccAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTc                                                                                                                        >  2:1395556/1‑150 (MQ=255)
                            ccAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTc                                                                                                                        <  1:47313/150‑1 (MQ=255)
                             cAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCg                                                                                                                       >  1:427686/1‑150 (MQ=255)
                                tGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGttt                                                                                                                    <  1:70150/150‑1 (MQ=255)
                                   agcaTGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGat                                                                                                                 <  1:549235/150‑1 (MQ=255)
                                    gcaGGTTGGTGGGTTCGCTGGTGTTGACCGACCCCTTGCCTGGTGTGCGGTCCGCCATGCCTTCGGGTTCCTGGCGGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGata                                                                                                                <  1:1374263/150‑1 (MQ=255)
                                        gTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGatacata                                                                                                            >  1:1899090/1‑150 (MQ=255)
                                          tggtggTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGAtacataca                                                                                                          >  2:734428/1‑150 (MQ=255)
                                               gtTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATAcagccag                                                                                                     >  2:1270782/1‑150 (MQ=255)
                                                tttCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATc                                                                                                                              >  2:46432/1‑124 (MQ=255)
                                                tttCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATc                                                                                                                              <  1:46432/124‑1 (MQ=255)
                                                      tGGTGTTGACCGCCCACTTGCCGACGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGGGCACCTTGTCGGTGACGGTG‑TCGATCTCCGTAACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGt                                                                                              >  1:911809/1‑150 (MQ=255)
                                                          gttgcgcGCCCACTTGGCGTCGGGGCTGACCGCCAGGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGGCCTCGTTTGATACATACAGCCAGCGGTCGTTggg                                                                                          <  2:293179/144‑1 (MQ=255)
                                                          gtTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTggg                                                                                          <  1:1355509/150‑1 (MQ=255)
                                                          gtTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTggg                                                                                          <  2:427686/150‑1 (MQ=255)
                                                           ttGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTgggg                                                                                         >  1:1541833/1‑150 (MQ=255)
                                                                cGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACGCCCAGCGGTCGTTGGGGTGCAg                                                                                    <  1:537752/150‑1 (MQ=255)
                                                                cGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAg                                                                                    >  2:1246235/1‑150 (MQ=255)
                                                                cGCCCACTGGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAg                                                                                    >  1:643387/1‑150 (MQ=255)
                                                                 gCCCACTTGCCGTCGGTGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCGCCTTGTCGGAGACGGTG‑TCGATCAACTTCACCAGCGCATCGTCCTCGTTTTATACATACAGCCAGCGGTCGTTGGGGTGCAgg                                                                                   >  1:1455117/1‑150 (MQ=255)
                                                                 gCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAgg                                                                                   <  1:1566332/150‑1 (MQ=255)
                                                                   ccACTTGCCGTCGGGGCTGACCGCCATGCCTTAGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATTGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGc                                                                                 >  1:1899984/1‑150 (MQ=255)
                                                                    cACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCg                                                                                >  2:1724409/1‑150 (MQ=255)
                                                                     aCTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGa                                                                               >  1:998717/1‑150 (MQ=255)
                                                                      cTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTGATCGA‑CACCGTCa                                                                                                                                        <  1:1800358/92‑1 (MQ=255)
                                                                      cTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTGATCGA‑CACCGTCa                                                                                                                                        >  2:1800358/1‑92 (MQ=255)
                                                                            gTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTc                                                                        >  1:219141/1‑150 (MQ=255)
                                                                                gggCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCgggg                                                                    <  2:1285746/150‑1 (MQ=255)
                                                                                 ggCTGACCGCCACGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGGGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGt                                                                   <  1:568332/150‑1 (MQ=255)
                                                                                   cTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTcc                                                                 <  2:1906532/150‑1 (MQ=255)
                                                                                   cTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTcc                                                                 >  1:165135/1‑150 (MQ=255)
                                                                                   cTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTcc                                                                 <  1:615362/150‑1 (MQ=255)
                                                                                       ccGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTCTGATACATACAGCCAGCGGTCGTTGGGCTGCAGGGCGAACTGCTCGGGGTCCTTGc                                                             <  2:998717/150‑1 (MQ=255)
                                                                                        cGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTAGGCTGCAGTGCGAACTGCTCGGGGTCCTTGcc                                                            <  2:841839/150‑1 (MQ=255)
                                                                                            aTGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCCTCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGaa                                                        <  2:219141/150‑1 (MQ=255)
                                                                                            aTGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCGTGCCGGaa                                                        <  2:78481/150‑1 (MQ=255)
                                                                                              gCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAgg                                                      >  1:1777836/1‑150 (MQ=255)
                                                                                               ccTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGc                                                     <  1:1518375/150‑1 (MQ=255)
                                                                                               ccTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGc                                                     <  2:920024/150‑1 (MQ=255)
                                                                                                 ttCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCggg                                                                     >  2:1874877/1‑132 (MQ=255)
                                                                                                 ttCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCggg                                                                     <  1:1874877/132‑1 (MQ=255)
                                                                                                 ttCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCaa                                                   >  1:1775673/1‑150 (MQ=255)
                                                                                                  tCGGGTTCGATGCCGACTTCGATCTCCCCGAGTATCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAAc                                                  <  2:1455117/150‑1 (MQ=255)
                                                                                                       ttCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCTCCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCtt                                             <  1:881462/150‑1 (MQ=255)
                                                                                                       ttCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCtt                                             <  1:1395556/150‑1 (MQ=255)
                                                                                                         cGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCtggg                                           <  2:656176/150‑4 (MQ=255)
                                                                                                         cGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTtga                                           <  2:1899090/150‑1 (MQ=255)
                                                                                                           aTGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGc                                                     >  1:1043984/1‑138 (MQ=255)
                                                                                                           aTGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGc                                                     <  2:1043984/138‑1 (MQ=255)
                                                                                                           aTGCCGACTTCGATCTGCCCGAGCACCTTGTCAGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTtgatg                                         >  2:1070978/1‑150 (MQ=255)
                                                                                                            tGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGaa                                                                              >  2:956970/1‑112 (MQ=255)
                                                                                                            tGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGaa                                                                              <  1:956970/112‑1 (MQ=255)
                                                                                                               cGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGcc                                                            >  2:1533567/1‑127 (MQ=255)
                                                                                                               cGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGcc                                                            <  1:1533567/127‑1 (MQ=255)
                                                                                                                 aCTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTg                                   <  1:180398/150‑1 (MQ=255)
                                                                                                                            cccGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCAc                        >  1:1521562/1‑150 (MQ=255)
                                                                                                                               gAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAg                                                                                    <  1:667010/87‑1 (MQ=255)
                                                                                                                               gAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAg                                                                                    >  2:667010/1‑87 (MQ=255)
                                                                                                                               gAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTc                     >  2:837773/1‑150 (MQ=255)
                                                                                                                                            gTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGatacata                                                                                                            >  1:1350287/1‑50 (MQ=255)
                                                                                                                                            gTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGatacata                                                                                                            <  2:1350287/50‑1 (MQ=255)
                                                                                                                                            gTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTACCGGAAGGCAACTCCTGGATGATCTGGCGTGTGGCCACGTCCATCACCTGCAcc        >  2:1754254/1‑150 (MQ=255)
                                                                                                                                             tGACGGTGATCGA‑CACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTg                                                                                            >  1:547695/1‑65 (MQ=37)
                                                                                                                                             tGACGGTGATCGA‑CACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTg                                                                                            <  2:547695/65‑1 (MQ=37)
                                                                                                                                             tGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCATCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCAccc       >  1:1251471/1‑150 (MQ=255)
                                                                                                                                             tGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCAccc       >  2:1643426/1‑150 (MQ=255)
                                                                                                                                              gACGGTG‑TCGATCACCGTCACCAGCGCATCGGCCGCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCg      <  1:679776/150‑1 (MQ=255)
                                                                                                                                               aCGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCgg     <  2:1541833/150‑1 (MQ=255)
                                                                                                                                                 ggTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTc   <  1:875159/150‑1 (MQ=255)
                                                                                                                                                 ggTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTc   >  2:182870/1‑150 (MQ=255)
                                                                                                                                                 ggTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACATGCACCCGGTc   >  1:1022499/1‑150 (MQ=255)
                                                                                                                                                  gtg‑tCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTTCTTTACGATATGCTTCTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTcc  <  2:1699262/150‑1 (MQ=255)
                                                                                                                                                  gtg‑tCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTcc  <  1:182870/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CGGCCAGGGTCTGGGTGCTGGTGTCGATCCAGTGCAGCATGTTGGTGGTTTCGCTGGTGTTGACCGCCCACTTGCCGTCGGGGCTGACCGCCATGCCTTCGGGTTCGATGCCGACTTCGATCTGCCCGAGCACCTTGTCGGTGACGGTG‑TCGATCACCGTCACCAGCGCATCGTCCTCGTTTGATACATACAGCCAGCGGTCGTTGGGGTGCAGGGCGAACTGCTCGGGGTCCTTGCCGGAAGGCAACTCCTTGATGATCTTGCGGGTGGCCACGTCCATCACCTGCACCCGGTCC  >  NC_002947/3056853‑3057148

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