Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 848,510 T→C 8.3% K21K (AAA→AAG PP_0730 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947848,5100TC 8.3% 245.3 / 8.6 85K21K (AAA→AAGPP_0730hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (35/42);  new base C (4/3);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 6.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.71e-01

CTGAATTCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCATA  >  NC_002947/848361‑848659
                                                                                                                                                      |                                                                                                                                                     
cTGAATTCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGt                                                                                                                                                       <  1:440341/150‑1 (MQ=255)
   aaTTCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTg                                                                                                                                                    <  2:994316/150‑1 (MQ=255)
    aTTCGCCGCCTTCCAGATTACGGATACGCCAATAGGGGTCCTACGCCGTTCGATAAGGCTCTATACGGGATGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTgg                                                                                                                                                   <  2:529652/150‑1 (MQ=255)
     ttCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGc                                                                                                                                                  >  1:268068/1‑150 (MQ=255)
       cgccgcCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGttt                                                                                                                                                     <  1:1491550/145‑1 (MQ=255)
       cgccgcCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGttt                                                                                                                                                     >  2:1491550/1‑145 (MQ=255)
        gccgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGc                                                                                                                                               >  2:651964/1‑150 (MQ=255)
        gccgccTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGc                                                                                                                                               >  1:1365482/1‑150 (MQ=255)
               tCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGCCGTGGCAGCCAAGCCGTTTGGCTGCCTCaaaa                                                                                                                                        <  1:1281334/150‑1 (MQ=255)
                 cAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGc                                                                                                                                      >  2:578901/1‑150 (MQ=255)
                   gATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCca                                                                                                                                    >  1:17459/1‑150 (MQ=255)
                       aTCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCcagcag                                                                                                                                <  1:578901/150‑1 (MQ=255)
                         cGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGc                                                                                                                              <  1:167392/150‑1 (MQ=255)
                                 aaaTAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACg                                                                                                                      >  2:26718/1‑150 (MQ=255)
                                          tCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGACCCACGCGTCGCGCCGATGCATGAAGGTGCCGTGGCAGCCAAGCCGGTTGGCTGCGTCAAAAGCCAGCAGGCCGATTACGGTGCTGTcg                                                                                                             >  1:906230/1‑150 (MQ=255)
                                           ccACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTcga                                                                                                            <  1:651964/150‑1 (MQ=255)
                                           ccACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGACCATGACCGGGCAGGCCGACATGCAACTGCCGGTCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTcga                                                                                                            <  1:285097/150‑1 (MQ=255)
                                            cacaTATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTcgac                                                                                                           >  1:1195289/1‑150 (MQ=255)
                                               atatCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACg                                                                                                                      <  2:1364782/136‑1 (MQ=255)
                                               atatCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACg                                                                                                                      >  1:1364782/1‑136 (MQ=255)
                                                    gCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGGAGGCCGCTTACGGTGCTGTCGACGAAAACCa                                                                                                   <  2:17459/150‑1 (MQ=255)
                                                      cAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCAc                                                                                                  <  1:199930/150‑1 (MQ=255)
                                                      cAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCAc                                                                                                  <  1:244068/150‑1 (MQ=255)
                                                        ggTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGc                                                                                                                                               >  2:2390132/1‑103 (MQ=255)
                                                        ggTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGc                                                                                                                                               <  1:2390132/103‑1 (MQ=255)
                                                        ggTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACgt                                                                                                <  1:26718/150‑1 (MQ=255)
                                                          tAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGAGGCAGCCAAACGGCTTGGCTGcca                                                                                                                                             <  2:98987/103‑2 (MQ=38)
                                                          tAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGAGGCAGCCAAACGGCTTGGCTGcaa                                                                                                                                             >  1:98987/1‑101 (MQ=255)
                                                              gCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCa                                                                                          <  2:1365482/150‑1 (MQ=255)
                                                              gCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCa                                                                                          >  1:1550972/1‑150 (MQ=255)
                                                              gCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCa                                                                                          >  2:941817/1‑150 (MQ=255)
                                                                tCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCAtt                                                                                        >  1:1655838/1‑150 (MQ=255)
                                                                tCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCAtt                                                                                        >  2:2184963/1‑150 (MQ=255)
                                                                    gaCGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGaa                                                                                    <  2:1550972/150‑1 (MQ=255)
                                                                        ggAGGTATCGTCCTTGATCGAGATCGTCAGGAACGAGGTGTTGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTc                                                                                <  1:22865/150‑1 (MQ=255)
                                                                         gAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCTCTTCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCg                                                                               <  1:1934167/150‑1 (MQ=255)
                                                                             tgtCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCt                                                                                                                                                 >  2:730977/1‑80 (MQ=255)
                                                                             tgtCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCt                                                                                                                                                 <  1:730977/80‑1 (MQ=255)
                                                                             tgtCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCtt                                                                                                                         >  2:72576/1‑104 (MQ=255)
                                                                             tgtCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCtt                                                                                                                         <  1:72576/104‑1 (MQ=255)
                                                                               tcgtcCTTTATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCca                                                                                                                                    <  2:325784/91‑1 (MQ=255)
                                                                               tcgtcCTTTATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCca                                                                                                                                    >  1:325784/1‑91 (MQ=255)
                                                                               tcgtcCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCAtt                                                                                        <  1:301009/135‑1 (MQ=255)
                                                                               tcgtcCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCAtt                                                                                        >  2:301009/1‑135 (MQ=255)
                                                                                  tcCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGGTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAg                                                                      >  2:741925/1‑150 (MQ=255)
                                                                                  tcCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACgt                                                                                                <  2:188212/124‑1 (MQ=255)
                                                                                  tcCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACgt                                                                                                >  1:188212/1‑124 (MQ=255)
                                                                                  tcCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCAtt                                                                                        >  2:2304188/1‑132 (MQ=255)
                                                                                  tcCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCAtt                                                                                        <  1:2304188/132‑1 (MQ=255)
                                                                                  tcCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTAAAAAGCCAGCAGGCCGCTTTCGGTGCTGTGGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAg                                                                      >  2:272141/1‑150 (MQ=255)
                                                                                    cTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAggg                                                                    >  1:73416/1‑150 (MQ=255)
                                                                                     ttGATAGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGCCGTGGCAGCCAAGCCTTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTGCATTGGAACCTCGGATTACCAGGGa                                                                   <  2:966644/150‑1 (MQ=255)
                                                                                      tGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGATGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGAt                                                                  <  1:595254/150‑1 (MQ=255)
                                                                                          cGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGaaa                                                              <  1:2184963/150‑1 (MQ=255)
                                                                                           gagaTCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAgg                                                                                                                               >  1:741632/1‑84 (MQ=255)
                                                                                           gagaTCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAgg                                                                                                                               <  2:741632/84‑1 (MQ=255)
                                                                                               tCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTaaa                                                         <  2:1655838/150‑1 (MQ=255)
                                                                                                cGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAAt                                                        <  2:1195289/150‑1 (MQ=255)
                                                                                                  tCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAg                                                      >  2:2309832/1‑150 (MQ=255)
                                                                                                    aGGAACGAGGCGGCGCGCTGATGAAGGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTg                                                    >  2:31411/1‑150 (MQ=255)
                                                                                                     ggAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGCTGAAAGTAAATAGTGc                                                   >  2:2243328/1‑150 (MQ=255)
                                                                                                      gAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTAc                                                                         >  1:1326500/1‑127 (MQ=255)
                                                                                                      gAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTAc                                                                         <  2:1326500/127‑1 (MQ=255)
                                                                                                      gAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCt                                                  <  1:941817/150‑1 (MQ=255)
                                                                                                        aCGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGaaa                                                              <  1:1636780/136‑1 (MQ=255)
                                                                                                        aCGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGaaa                                                              >  2:1636780/1‑136 (MQ=255)
                                                                                                          gAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGc                                              <  1:2309832/150‑1 (MQ=255)
                                                                                                             gCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCg                                           >  2:721724/1‑150 (MQ=255)
                                                                                                                       gatgaATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCt                                                                                 >  1:977294/1‑102 (MQ=255)
                                                                                                                       gatgaATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCt                                                                                 <  2:977294/102‑1 (MQ=255)
                                                                                                                           aatgaaGGTGACGTGGCAGCCAAGCCGTTTGTCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCa                                                                       >  2:1486281/1‑108 (MQ=255)
                                                                                                                           aatgaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCa                                                                       <  1:1486281/108‑1 (MQ=255)
                                                                                                                           aatgaaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaa                             <  1:2275124/150‑1 (MQ=255)
                                                                                                                               aaGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaagcaa                         <  2:73416/150‑1 (MQ=255)
                                                                                                                                 ggTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTaaa                                                         <  1:1719552/116‑1 (MQ=255)
                                                                                                                                 ggTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTaaa                                                         >  2:1719552/1‑116 (MQ=255)
                                                                                                                                     aCGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTc                                                                                           >  1:1049616/1‑78 (MQ=255)
                                                                                                                                     aCGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTc                                                                                           <  2:1049616/78‑1 (MQ=255)
                                                                                                                                        tGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGa                                                                   >  1:1818925/1‑99 (MQ=255)
                                                                                                                                        tGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGa                                                                   <  2:1818925/99‑1 (MQ=255)
                                                                                                                                          gCAGCCAAACGGCTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAg                                                                      <  1:2264728/94‑1 (MQ=255)
                                                                                                                                          gCAGCCAAACGGCTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAg                                                                      >  2:2264728/1‑94 (MQ=255)
                                                                                                                                          gCAGCCAAACGGCTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaag                            <  1:2002008/136‑1 (MQ=255)
                                                                                                                                          gCAGCCAAACGGCTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaag                            >  2:2002008/1‑136 (MQ=255)
                                                                                                                                          gCAGCCAAACGGCTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGa              >  1:1050101/1‑150 (MQ=255)
                                                                                                                                            aGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaagcaa                         >  2:552361/1‑137 (MQ=255)
                                                                                                                                            aGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaagcaa                         <  1:552361/137‑1 (MQ=255)
                                                                                                                                                aaGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTgc        <  1:2049240/150‑1 (MQ=255)
                                                                                                                                                      tttGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaa                             >  1:2195710/1‑123 (MQ=255)
                                                                                                                                                      tttGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGaa                             <  2:2195710/123‑1 (MQ=255)
                                                                                                                                                      tttGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCATa  <  1:741925/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CTGAATTCGCCGCCTTCCAGATTATCGATACGCAAATAAGCATCCACATATCGCAGGTAAGGCTCGAGACGGGAGGTGTCGTCCTTGATCGAGATCGTCAGGAACGAGGCGTCGCGCTGATGAATGAAGGTGACGTGGCAGCCAAGCCGTTTGGCTGCCTCAAAAGCCAGCAGGCCGCTTACGGTGCTGTCGACGAAAACCACGTGTTTCATTGGAACCTCGGATTACCAGGGATGAAAGTAAATAGTGCTGAGCCCGAACAGCATTGCGAAGCAAGCGACCGTCGAAGGTGCGCCATA  >  NC_002947/848361‑848659

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