Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,918,010 G→A 13.7% V122I (GTC→ATC) ‡ PP_4328 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,918,0100GA13.7% 157.6 / 15.4 66V122I (GTC→ATC) ‡PP_4328hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (31/25);  new base A (4/5);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

CTTCAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCG  >  NC_002947/4917863‑4918147
                                                                                                                                                   |                                                                                                                                         
cTTCAGGCCAACGCCAGCCAGCCGTTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTc                                                                                                                                         <  1:975105/150‑1 (MQ=255)
      gCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCAcc                                                                                                                                   >  1:90915/1‑150 (MQ=255)
       ccAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCt                                                                                                                                  <  1:1368516/150‑1 (MQ=255)
        cAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCtt                                                                                                                                 >  1:394863/1‑150 (MQ=255)
         aaCGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtg                                                                                                                                <  2:394863/150‑1 (MQ=255)
           cGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtgct                                                                                                                              <  2:1288478/150‑1 (MQ=255)
            gccagccaGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTtgcct                                                                                                                             <  1:769118/150‑3 (MQ=255)
               agccagccGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACACTCACCCAGCTGGACACCATCAATGTGCCGGTCGGCACCTTGCTgcaa                                                                                                                          <  2:430991/150‑1 (MQ=255)
                       gCTGCCCTGGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTAAGAACCTTGCTGCAAGGCAAGGt                                                                                                                  >  2:423348/1‑150 (MQ=255)
                            cccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGa                                                                                                             >  1:548857/1‑150 (MQ=255)
                            cccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGa                                                                                                             <  1:582363/150‑1 (MQ=255)
                             ccAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGTTCGGCACCTTGCTGCAAGGCAAGGTCCTGac                                                                                                            <  2:1571328/150‑1 (MQ=255)
                                gggCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGaccaa                                                                                                         >  2:1505887/1‑150 (MQ=255)
                                 ggCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCTGACGGCTACCCTCACTCATCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGaccaac                                                                                                        <  1:197881/150‑1 (MQ=255)
                                        aGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAgctgc                                                                                                 >  1:1408729/1‑150 (MQ=255)
                                                gTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATTCTGCAACAGGCCAGCGCAAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTCCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAg                                                                                         <  2:549685/150‑1 (MQ=255)
                                                 tCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAgg                                                                                        >  2:532980/1‑150 (MQ=255)
                                                       aGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACACg                                                                                  >  1:490028/1‑150 (MQ=255)
                                                        gACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGACAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCgg                                                                                 >  1:388033/1‑150 (MQ=255)
                                                           cGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAgc                                                                                                    <  2:1070421/128‑1 (MQ=255)
                                                           cGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAgc                                                                                                    >  1:1070421/1‑128 (MQ=255)
                                                            gAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtgct                                                                                                                              <  1:567008/101‑1 (MQ=255)
                                                            gAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtgct                                                                                                                              >  2:567008/1‑101 (MQ=255)
                                                               aGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCCCACAGCTAGAAACCAGCAAAATGAGGGTCGGCACCTTGCAGCAACGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCg                                                                          >  2:595771/1‑150 (MQ=255)
                                                                      gCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAg                                                                                         <  2:291004/128‑1 (MQ=255)
                                                                      gCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAg                                                                                         >  1:291004/1‑128 (MQ=255)
                                                                        cTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTAc                                                                 >  2:303908/1‑150 (MQ=255)
                                                                          ggCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCGAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCg                                                               <  1:1546645/150‑1 (MQ=255)
                                                                           gCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCgg                                                              >  1:907228/1‑150 (MQ=255)
                                                                                 aTGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTg                                                        >  1:829193/1‑150 (MQ=255)
                                                                                 aTGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTg                                                        <  2:907228/150‑1 (MQ=255)
                                                                                    ctgcAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTgcaa                                                                                                                          <  1:954171/81‑1 (MQ=255)
                                                                                    ctgcAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTgcaa                                                                                                                          >  2:954171/1‑81 (MQ=255)
                                                                                        aaCAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCgg                                                              <  2:1017213/137‑1 (MQ=255)
                                                                                        aaCAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCgg                                                              >  1:1017213/1‑137 (MQ=255)
                                                                                        aaCAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtgct                                                                                                                              >  1:370305/1‑73 (MQ=38)
                                                                                        aaCAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGACCGGCACCTtgct                                                                                                                              <  2:370305/73‑1 (MQ=38)
                                                                                           aGGCCAGCGACAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCTGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGc                                              <  2:388033/150‑1 (MQ=255)
                                                                                              ccagcgccagcCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCa                                           >  2:1137166/1‑150 (MQ=255)
                                                                                               cagcgccagcCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCaa                                          >  1:725111/1‑150 (MQ=255)
                                                                                                 gcgccagcCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACa                                        <  2:548857/150‑1 (MQ=255)
                                                                                                   gccagccaGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGc                                      >  1:533149/1‑150 (MQ=255)
                                                                                                        ccagGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTg                                                        <  2:1410998/127‑1 (MQ=255)
                                                                                                        ccagGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTg                                                        >  1:1410998/1‑127 (MQ=255)
                                                                                                           ggTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTgcaaggcaag                                                                                                                    <  1:750891/64‑1 (MQ=255)
                                                                                                           ggTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTgcaaggcaag                                                                                                                    >  2:750891/1‑64 (MQ=255)
                                                                                                                    ccTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCcac                     >  1:310067/1‑150 (MQ=255)
                                                                                                                     cTCGCCCAGCTGGACACCGGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCcaca                    <  2:1595683/150‑1 (MQ=255)
                                                                                                                      tCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGcc                                                                       >  2:1044495/1‑98 (MQ=255)
                                                                                                                      tCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGcc                                                                       <  1:1044495/98‑1 (MQ=255)
                                                                                                                          ccAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGcc                             <  1:830830/136‑1 (MQ=255)
                                                                                                                          ccAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGcc                             >  2:830830/1‑136 (MQ=255)
                                                                                                                              ccggACACCAGCAAGGTGCCCGTCGGCACCTTGCTGCTAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATc           <  1:423348/148‑1 (MQ=255)
                                                                                                                                  acacCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACa       >  1:1254847/1‑150 (MQ=255)
                                                                                                                                      cAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAAc                                   >  2:913657/1‑118 (MQ=255)
                                                                                                                                      cAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAAc                                   <  1:913657/118‑1 (MQ=255)
                                                                                                                                      cAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgccc   >  2:1239192/1‑150 (MQ=255)
                                                                                                                                      cAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgccc   <  2:725111/150‑1 (MQ=255)
                                                                                                                                      cAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgccc   >  1:886132/1‑150 (MQ=255)
                                                                                                                                       aGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAgcccg  >  1:1478336/1‑150 (MQ=255)
                                                                                                                                        gCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCcaca                    >  1:1504550/1‑131 (MQ=255)
                                                                                                                                        gCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCcaca                    <  2:1504550/131‑1 (MQ=255)
                                                                                                                                        gCAAGGTGCCGACCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGa                >  2:1073087/1‑135 (MQ=255)
                                                                                                                                        gCAAGGTGCCGACCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGa                <  1:1073087/135‑1 (MQ=255)
                                                                                                                                            ggTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAg                                >  2:174059/1‑115 (MQ=255)
                                                                                                                                            ggTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAg                                <  1:174059/115‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
CTTCAGGCCAACGCCAGCCAGCCGCTGCCCCAGGGCAGCCAGATCACCGTCAGCCAGACCGAAAGCAATCGCCTGGCGATCATGCTGCAACAGGCCAGCGCCAGCCAGGTCGCCACCCTCACCCAGCTGGACACCAGCAAGGTGCCGGTCGGCACCTTGCTGCAAGGCAAGGTCCTGACCAACCAGCTGCTGGCCCAGGCACCCGGCCAACCGGCCGTGTACCGGGCGCTGGTCAGCCTGCTCAACAGCAACCAGGCCGGTGCCACACTGAGCATCGACAGCCCG  >  NC_002947/4917863‑4918147

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: