Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 980,749 T→G 6.2% L64R (CTG→CGG)  hscB → DnaJ‑like molecular chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947980,7490TG 6.2% 309.6 / 9.4 97L64R (CTG→CGG) hscBDnaJ‑like molecular chaperone
Reads supporting (aligned to +/- strand):  ref base T (49/42);  new base G (3/3);  total (52/45)
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.22e-01

TCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTT  >  NC_002947/980609‑980897
                                                                                                                                             |                                                                                                                                                      
tCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTg                                                                                                                                                >  7:261441/1‑149 (MQ=255)
    gATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCGGCGCACTgcgca                                                                                                                                            <  5:239931/149‑2 (MQ=255)
      tctcGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGcgccgc                                                                                                                                          >  7:11728/1‑149 (MQ=255)
      tctcGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGcgccgc                                                                                                                                          >  1:324392/1‑149 (MQ=255)
      tctcGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGcgccgc                                                                                                                                          >  8:17716/1‑149 (MQ=255)
          gACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGCTCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGcgccgc                                                                                                                                          <  5:345180/145‑1 (MQ=255)
          gACAAGCTGGCCACACGCTAACGCCAGCAGGCCCGCGAAGTCCATCCTAATCGCTATGCCGACGCCTCCGAGCGTGAAAAGAGTATTGATCTGCAAAAATACCCGCCCCTCAACCACGCTTACCAGA‑CCCTCCGCAGTGCGCcgcgaag                                                                                                                                      >  4:18181/1‑146 (MQ=255)
               gCTGGCCACTCGCTATCGCGAGCTGGCCCGCCAAGTCCATACTGATCGATTTGCCGACGCCACCCATCGTCAGCACCGTATTGCGCTGGAAAAATACGCGGCCCTCAACCACGCTTACAAGA‑CGCTGCGCAGTGCGCCGCGTcgcccgc                                                                                                                                 >  2:244102/1‑149 (MQ=255)
                    ccACTCGCTATCGCGAGCTGGGCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTAc                                                                                                                            <  5:200260/149‑1 (MQ=255)
                    ccACTCGCTATCGCGAGCTGGCCCGCTAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTAc                                                                                                                            <  3:18181/149‑1 (MQ=255)
                     cACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTAcc                                                                                                                           >  2:196914/1‑149 (MQ=255)
                        tCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgc                                                                                                                        >  3:98558/1‑149 (MQ=255)
                          gCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgctg                                                                                                                      <  3:308107/149‑1 (MQ=255)
                          gCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACctgctg                                                                                                                      >  3:184643/1‑149 (MQ=255)
                                 cgAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCAt                                                                                                                 >  4:45099/1‑147 (MQ=255)
                                 cgAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCa                                                                                                               <  7:215850/149‑1 (MQ=255)
                                  gAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAg                                                                                                              <  4:184643/149‑1 (MQ=255)
                                      tGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGc                                                                                                          >  7:301636/1‑149 (MQ=255)
                                      tGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGC‑TACCAGA‑CGCGGCGCACTGCGCAGCGTc                                                                                                                                       >  5:86076/1‑119 (MQ=255)
                                      tGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGC‑TACCAGA‑CGCGGCGCACTGCGCAGCGTc                                                                                                                                       <  6:86076/119‑1 (MQ=255)
                                        gCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCa                                                                                                        <  4:212561/149‑1 (MQ=255)
                                           cgcgAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGACCGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACg                                                                                                      <  1:259150/149‑1 (MQ=255)
                                                     aTCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTAc                                                                                                                            <  5:222627/116‑1 (MQ=255)
                                                     aTCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTAc                                                                                                                            >  6:222627/1‑116 (MQ=255)
                                                     aTCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCag                                                                                         >  3:220624/1‑149 (MQ=255)
                                                      tCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCag                                                                                         >  7:276568/1‑148 (MQ=255)
                                                        cTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCaggag                                                                                      <  7:17716/149‑1 (MQ=255)
                                                          gATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGt                                                                                    >  4:184932/1‑149 (MQ=255)
                                                             cGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCa                                                                                  <  8:11728/148‑1 (MQ=255)
                                                               cTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACgg                                                                               >  4:202633/1‑149 (MQ=255)
                                                                tttGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAAGTCAGGGt                                                                              >  1:304701/1‑149 (MQ=255)
                                                                  tGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTc                                                                             <  2:324392/148‑1 (MQ=255)
                                                                     cGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACg                                                                         >  7:242030/1‑149 (MQ=255)
                                                                        cGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGAcc                                                                      <  3:184932/149‑1 (MQ=255)
                                                                         gCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGAccc                                                                     >  2:271121/1‑149 (MQ=255)
                                                                         gCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGAcc                                                                      >  6:162689/1‑148 (MQ=255)
                                                                          ccTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCt                                                                    >  1:68251/1‑149 (MQ=255)
                                                                            tCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCaggag                                                                                      >  6:167885/1‑129 (MQ=255)
                                                                            tCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCaggag                                                                                      <  5:167885/129‑1 (MQ=255)
                                                                              cGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGa                                                                  <  8:261441/147‑1 (MQ=255)
                                                                              cGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACt                                                                >  6:254415/1‑149 (MQ=255)
                                                                               gAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTc                                                                             >  4:204113/1‑135 (MQ=255)
                                                                               gAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTc                                                                             <  3:204113/135‑1 (MQ=255)
                                                                                aGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTc                                                              <  3:45099/149‑1 (MQ=255)
                                                                                   gTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCctg                                                           >  2:115056/1‑149 (MQ=255)
                                                                                    ttatcagcGTATTTCGCTGTAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAAGTCAGGGTCGTGGACCGTGACCTCctg                                                           <  2:304701/144‑1 (MQ=255)
                                                                                    tGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAAGTCAGGGTCGTGGACCGTGACCTCctgg                                                          <  4:98558/149‑2 (MQ=255)
                                                                                     gAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTcg                                                                                                                                      <  2:270637/74‑1 (MQ=255)
                                                                                     gAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCctgct                                                         >  7:289428/1‑149 (MQ=255)
                                                                                      agcagcGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCctgct                                                         <  8:289428/148‑1 (MQ=255)
                                                                                        cagcGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCctgctg                                                        >  3:269799/1‑147 (MQ=255)
                                                                                        cagcGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTgca                                                      <  1:82632/149‑1 (MQ=255)
                                                                                         agcGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTgcag                                                     <  4:269799/149‑1 (MQ=255)
                                                                                           cGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTgcagca                                                   <  5:254415/149‑1 (MQ=255)
                                                                                                  cgcTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTgg                                                                                                                     <  1:127845/78‑1 (MQ=255)
                                                                                                  cgcTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTgg                                                                                                                     >  2:127845/1‑78 (MQ=255)
                                                                                                  cgcTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAg                                            >  7:127558/1‑149 (MQ=255)
                                                                                                   gcTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGt                                           >  7:12515/1‑149 (MQ=255)
                                                                                                       gAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGGAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGc                                        <  1:115056/148‑1 (MQ=255)
                                                                                                       gAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGCCTTCCTGCTGCAGCAGATGCAGTGGCg                                       >  7:311374/1‑149 (MQ=255)
                                                                                                       gAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGc                                        <  5:191582/148‑1 (MQ=255)
                                                                                                       gAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCg                                       >  3:264918/1‑149 (MQ=255)
                                                                                                       gAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTAGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTgcagc                                                    >  2:46899/1‑136 (MQ=255)
                                                                                                       gAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTAGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTgcagc                                                    <  1:46899/136‑1 (MQ=255)
                                                                                                       gAAAAATCCGCCGCCCTCAACGAAGCTTACCAGC‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACGTAT‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGc                                        <  1:267168/149‑1 (MQ=255)
                                                                                                        aaaaaTCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTAGCTGCTGGCCATCAGCGGCGACG‑AA‑GTGCCCCAGGAGGTGACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGt                                      <  1:244102/149‑1 (MQ=255)
                                                                                                           aaTCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGaa                                   <  5:145904/149‑1 (MQ=255)
                                                                                                             tCCGCGGCCCTCAACGAGGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCCCGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAaga                                 >  5:342917/1‑149 (MQ=255)
                                                                                                             tCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGaa                                   >  5:107883/1‑147 (MQ=255)
                                                                                                             tCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGAGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGAGATCCTGCTGGAGCAGATGCAGTGGCGTGaa                                   >  7:138000/1‑147 (MQ=255)
                                                                                                              ccGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAagag                                >  1:80130/1‑149 (MQ=255)
                                                                                                               cgcgGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AACTTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAagag                                <  1:271121/149‑1 (MQ=255)
                                                                                                               cgcgGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAagag                                <  5:162689/148‑1 (MQ=255)
                                                                                                                 cgGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCAt                                                                                                                 >  8:163008/1‑67 (MQ=255)
                                                                                                                 cgGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCAt                                                                                                                 <  7:163008/67‑1 (MQ=255)
                                                                                                                      cTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCgaaga                        <  2:139542/149‑1 (MQ=255)
                                                                                                                       tCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCgaaga                        >  4:201640/1‑148 (MQ=255)
                                                                                                                        cAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCt                              >  4:223715/1‑141 (MQ=255)
                                                                                                                        cAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCt                              <  3:223715/141‑1 (MQ=255)
                                                                                                                         aaCGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACT‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACt                     <  8:301636/149‑1 (MQ=255)
                                                                                                                         aaCGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAAc                      <  2:68251/148‑1 (MQ=255)
                                                                                                                           cgacgCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTg                    >  3:121144/1‑148 (MQ=255)
                                                                                                                              cgCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAgg                <  3:201640/149‑1 (MQ=255)
                                                                                                                                 ttACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGAt              >  7:28386/1‑148 (MQ=255)
                                                                                                                                 ttACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGAt              >  2:99920/1‑148 (MQ=255)
                                                                                                                                    ccAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTgcagc                                                    <  8:76618/107‑1 (MQ=255)
                                                                                                                                    ccAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTgcagc                                                    >  7:76618/1‑107 (MQ=255)
                                                                                                                                    ccAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAg          <  2:80130/149‑1 (MQ=255)
                                                                                                                                    ccAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAg          >  4:274834/1‑149 (MQ=255)
                                                                                                                                     cAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGc         >  1:46378/1‑149 (MQ=255)
                                                                                                                                     cAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGc         >  4:233883/1‑149 (MQ=255)
                                                                                                                                       gA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCg       <  8:242030/149‑1 (MQ=255)
                                                                                                                                       gA‑CGCGGCGCACTGCGCAGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAg          <  7:186352/146‑1 (MQ=255)
                                                                                                                                       gA‑CGCGGCGCACTGCGCAGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAg          >  8:186352/1‑146 (MQ=255)
                                                                                                                                       gA‑CGCGGCGCACTGCGCAGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGcc        >  8:46471/1‑148 (MQ=255)
                                                                                                                                        a‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGa      >  1:296019/1‑149 (MQ=255)
                                                                                                                                            cTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCt   >  4:131488/1‑149 (MQ=255)
                                                                                                                                             tGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCtt  >  1:304768/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                      
TCTGGATCTCGACAAGCTGGCCACTCGCTATCGCGAGCTGGCCCGCGAAGTCCATCCTGATCGCTTTGCCGACGCCTCCGAGCGTGAGCAGCGTATTGCGCTGGAAAAATCCGCGGCCCTCAACGACGCTTACCAGA‑CGCTGCGCAGTGCGCCGCGTCGCGCGCGCTACCTGCTGGCCATCAGCGGCCACG‑AA‑GTGCCCCAGGAGGTCACGGTCCACGACCCTGACTTCCTGCTGCAGCAGATGCAGTGGCGTGAAGAGCTCGAAGAACTGCAGGATGAAGCCGACCTT  >  NC_002947/980609‑980897

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