Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,703,102 G→A 5.1% V95V (GTC→GTT)  era ← membrane‑associated, 16S rRNA‑binding GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,703,1020GA 5.1% 255.5 / 2.4 78V95V (GTC→GTT) eramembrane‑associated, 16S rRNA‑binding GTPase
Reads supporting (aligned to +/- strand):  ref base G (35/39);  new base A (2/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

GCTTGCCAGGAACTGCAGGTGCGGCAGCAGATCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGACGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAA  >  NC_000913/2702955‑2703247
                                                                                                                                                   |                                                                                                                                                 
gCTTGCCAGGAACTGCAGGTGCGGCAGCAGATCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGa                                                                                                                                                  <  1:332194/149‑1 (MQ=255)
       aGGAACTGCAGGTGCGGCAGCAGATCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAAt                                                                                                                                           <  2:24545/149‑1 (MQ=255)
       aGGAACTGCAGGTGCGGCAGCAGATCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACAACAAAAAt                                                                                                                                           >  7:325411/1‑149 (MQ=255)
       aGGAACTGCAGGTGCGGCAGCAGATCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACAACAAAAAt                                                                                                                                           >  3:13228/1‑149 (MQ=255)
        ggAACTGCAGGTGCGGCAGCAGATCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAAt                                                                                                                                           >  1:154104/1‑148 (MQ=255)
          aaCTGCAGGTGCGGCAGCAGATCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCAc                                                                                                                                        <  8:361679/149‑1 (MQ=255)
            cTGCAGTGGCGGCAGCAGAGAGGCTTTCTCCTGCACGTTGCCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTTAGCACCATTTCGGCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCa                                                                                                                                      <  1:109070/149‑1 (MQ=255)
                aGGTGCGGCAGCAGATCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTc                                                                                                                                  <  1:165033/149‑1 (MQ=255)
                  gTGCGGCAGCAGATCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGGTc                                                                                                                                  <  8:378093/147‑1 (MQ=255)
                  gTGCGGCAGCAGATCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTc                                                                                                                                  >  7:378093/1‑147 (MQ=255)
                   tGCGGCAGCAGATCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAAc                                                                                                                               <  7:203113/149‑1 (MQ=255)
                        cagcagATCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGCCAAAAATCACCAGCTCAACATCGc                                                                                                                          >  7:75237/1‑149 (MQ=255)
                        cagcagATCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACAACAAAAATCACCAGCTCAACATCGc                                                                                                                          <  4:13228/149‑1 (MQ=255)
                        cagcagATCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACAACAAAAATCACCAGCTCAACATCGc                                                                                                                          <  8:325411/149‑1 (MQ=255)
                             gATCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATa                                                                                                                     <  7:303757/149‑1 (MQ=255)
                              aTCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCg                                                                                                                           <  2:227879/142‑1 (MQ=255)
                              aTCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCg                                                                                                                           >  1:227879/1‑142 (MQ=255)
                              aTCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCaa                                                                                                                       >  2:191960/1‑146 (MQ=255)
                              aTCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATag                                                                                                                    >  2:75460/1‑149 (MQ=255)
                              aTCGGCTTTCTCCTGCACATTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGCGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCaa                                                                                                                       >  2:136764/1‑146 (MQ=255)
                               tCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCCCCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATa                                                                                                                     >  8:301102/1‑147 (MQ=255)
                               tCGGCTTGCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTGTGCATTCGCGCAGTTTGTGTAGCACCATTTCGGCGTCCGGCGGCCAGCGGGGGCCTGCACCGACAAAAATCACCTGCTCAACATCGCCAATagc                                                                                                                   <  5:262507/149‑2 (MQ=255)
                                cGGCTTTCTCCTGCACTTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCTTCGTCCGGCGTCCAGCTGGTTCCTTCAACGACAAAAATCACCAGCTCAACATCTCCAATagag                                                                                                                  >  2:265821/1‑149 (MQ=255)
                                   cTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAgctg                                                                                                               >  6:151269/1‑149 (MQ=255)
                                       ctcCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTc                                                                                                            <  8:135340/148‑1 (MQ=255)
                                           tGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTcgccgc                                                                                                       >  7:236084/1‑149 (MQ=255)
                                             cACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAGCATCGCCAATAGAGCTGCTCGCCGCtt                                                                                                     >  8:86921/1‑149 (MQ=255)
                                             cACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCtt                                                                                                     >  6:200325/1‑149 (MQ=255)
                                                 ttGTCCACTTTGTTCACCGCGAGGATTAACGGCGCTTTGCCTTCGCGCATTTTGTTGAGCACAATTTCGTCGTCCGGCGTCCAGCGGGTGACTTCAGAGACAAAAATCACAAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGtt                                                                                                 >  5:40650/1‑149 (MQ=255)
                                                   gTCCACTTTGTTCACCGCGATGATTACCGGCGCTTTGCCTTCTCGCAGTTTGTTTAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTtc                                                                                                <  1:265821/148‑1 (MQ=255)
                                                   gTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTtc                                                                                                <  2:370841/148‑1 (MQ=255)
                                                    tCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTtcat                                                                                              <  8:236084/149‑1 (MQ=255)
                                                      cACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAgc                                                                                                                 <  5:84904/128‑1 (MQ=255)
                                                      cACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAgc                                                                                                                 >  6:84904/1‑128 (MQ=255)
                                                      cACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTtcatc                                                                                             <  2:9609/148‑1 (MQ=255)
                                                      cACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTtcat                                                                                              >  5:75916/1‑147 (MQ=255)
                                                       aCTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAg                                                                                           <  8:98545/149‑1 (MQ=255)
                                                            gTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCa                                                                                                                                      >  8:230219/1‑101 (MQ=255)
                                                            gTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCa                                                                                                                                      <  7:230219/101‑1 (MQ=255)
                                                             ttCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGt                                                                                      <  2:154104/148‑1 (MQ=255)
                                                                 ccGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATg                                                                                 <  7:194202/149‑1 (MQ=255)
                                                                  cgcgAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACaaaaa                                                                                                                                            <  2:55086/89‑1 (MQ=255)
                                                                  cgcgAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAACTCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATgg                                                                                >  5:296094/1‑149 (MQ=255)
                                                                         aTTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGttt                                                                         <  1:136764/149‑1 (MQ=255)
                                                                          ttACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGtttt                                                                        >  2:263232/1‑149 (MQ=255)
                                                                             ccGGCGCTTTGCCTTCGCTCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCATTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTttct                                                                     <  6:40650/149‑1 (MQ=255)
                                                                             ccGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCTGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTttct                                                                     >  6:133411/1‑149 (MQ=255)
                                                                                   cTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCata                                                               <  6:66042/149‑1 (MQ=255)
                                                                                    tttGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCatat                                                              >  3:314380/1‑149 (MQ=255)
                                                                                      tGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATg                                                             <  5:151269/148‑1 (MQ=255)
                                                                                        ccTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAg                                                          <  8:26758/149‑1 (MQ=255)
                                                                                         cTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTcgccg                                                                                                        <  8:59572/102‑1 (MQ=255)
                                                                                         cTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAgg                                                         >  4:258391/1‑149 (MQ=255)
                                                                                          ttCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGc                                                        >  1:234510/1‑149 (MQ=255)
                                                                                          ttCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGc                                                        >  3:182737/1‑149 (MQ=255)
                                                                                               gcAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGgt                                                   <  5:200325/149‑1 (MQ=255)
                                                                                                   tttGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTAt                                               <  6:357167/149‑1 (MQ=255)
                                                                                                       ttgAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGGTCAACATCGCGAATAGAGCTGCTCGCCGGTTTGTTCATGAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGAc                                           >  3:9242/1‑149 (MQ=255)
                                                                                                       ttgAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGAc                                           <  1:263232/149‑1 (MQ=255)
                                                                                                         gAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGACGt                                         >  7:3061/1‑149 (MQ=255)
                                                                                                            cACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGACg                                          >  7:228175/1‑145 (MQ=255)
                                                                                                             aCCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGACGTAGAt                                     <  7:41551/149‑1 (MQ=255)
                                                                                                             aCCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGGGTATCGACGTAGAt                                     <  3:258391/149‑1 (MQ=255)
                                                                                                               cATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGgt                                                   >  4:77886/1‑133 (MQ=255)
                                                                                                               cATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGgt                                                   <  3:77886/133‑1 (MQ=255)
                                                                                                                   tcgtcgTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGACGTAGATCGCCt                                <  1:191960/148‑1 (MQ=255)
                                                                                                                     gtcgtcCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGGTTTTCTTCCATATGCAGGCCCGGTGTATCGACGTAGATCGCCTGa                              <  4:26417/148‑1 (MQ=255)
                                                                                                                        gtcCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGACGTAGATCGCCTGATAcg                          >  5:352345/1‑149 (MQ=255)
                                                                                                                        gtcCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGACGTAGATCGCCTGATAcg                          >  5:346694/1‑149 (MQ=255)
                                                                                                                                gccAGCGGGTGCCTTCAACGACAAAAATCCCCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTGTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGACGTAGATCGCCTGATACGCGCCTTCa                  <  7:381489/148‑1 (MQ=255)
                                                                                                                                  cAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCGGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGACGTAGATCGCCTGATACGCGCCTTCAg                 >  2:400808/1‑148 (MQ=255)
                                                                                                                                     cGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGACGTAGATCGCCTGATACGCGCCTTCAGTATg             <  8:182225/149‑1 (MQ=255)
                                                                                                                                       ggTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGACGTAGATCGCCTGATACGCGCCTTCAGTATGGa           <  2:234510/149‑1 (MQ=255)
                                                                                                                                            cTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGACGTAGATCGCCTGATACGCGCCTTCAGTATGGATccc       >  5:125483/1‑148 (MQ=255)
                                                                                                                                             ttCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGACGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCca     >  3:338219/1‑149 (MQ=255)
                                                                                                                                             ttCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGACGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCca     >  1:13668/1‑149 (MQ=255)
                                                                                                                                               cAAAGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCAGGTGTATCGACGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCcac    >  5:140817/1‑148 (MQ=255)
                                                                                                                                                aaCGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGACGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACaa  <  1:75460/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GCTTGCCAGGAACTGCAGGTGCGGCAGCAGATCGGCTTTCTCCTGCACGTTGTCCACTTTGTTCACCGCGAGGATTACCGGCGCTTTGCCTTCGCGCAGTTTGTTGAGCACCATTTCGTCGTCCGGCGTCCAGCGGGTGCCTTCAACGACAAAAATCACCAGCTCAACATCGCCAATAGAGCTGCTCGCCGCTTTGTTCATCAGGCGGTTAATGGCGCGTTTTTCTTCCATATGCAGGCCCGGTGTATCGACGTAGATCGCCTGATACGCGCCTTCAGTATGGATCCCCACAA  >  NC_000913/2702955‑2703247

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