Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,355,456 T→A 6.6% K11N (AAA→AAT) ‡ PP_1176 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,355,4560TA 6.6% 846.4 / 34.2 257K11N (AAA→AAT) ‡PP_1176hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (109/130);  new base A (6/11);  total (116/141)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 266 total aligned reads displayed.

GGAGAGCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTAGTCCC  >  NC_002947/1355308‑1355605
                                                                                                                                                    |                                                                                                                                                     
ggAGAGCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAtt                                                                                                                                                      >  1:3274426/1‑150 (MQ=255)
 gagaGCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAttt                                                                                                                                                     <  2:4818319/150‑1 (MQ=255)
 gagaGCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAttt                                                                                                                                                     <  1:2004482/150‑1 (MQ=255)
 gagaGCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAttt                                                                                                                                                     <  2:5512574/150‑1 (MQ=255)
  agagCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTg                                                                                                                                                    <  1:5396537/150‑1 (MQ=255)
     gCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCa                                                                                                                                                 <  1:3051112/150‑1 (MQ=255)
       cGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGc                                                                                                                                               <  1:3571500/150‑1 (MQ=255)
        gTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGcc                                                                                                                                              >  1:2777543/1‑150 (MQ=255)
          ttCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCtt                                                                                                                                            <  1:2068958/150‑1 (MQ=255)
          ttCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGCCtt                                                                                                                                            <  2:3599329/150‑1 (MQ=255)
           tCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGCCttt                                                                                                                                           <  2:5457407/150‑2 (MQ=255)
            cTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGCCtttt                                                                                                                                          <  1:4311295/150‑3 (MQ=255)
            cTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAAATGGCAGCCtttt                                                                                                                                          <  2:4067581/150‑3 (MQ=255)
             tGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCttgtt                                                                                                                                         <  2:660316/150‑1 (MQ=255)
             tGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCttgtt                                                                                                                                         <  1:2299367/150‑1 (MQ=255)
              gACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCttgttg                                                                                                                                        >  2:1537902/1‑150 (MQ=255)
                 gagaACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGc                                                                                                                                    >  1:1087291/1‑150 (MQ=255)
                   gaACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCt                                                                                                                                   >  2:3544561/1‑150 (MQ=255)
                        ggATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGcc                                                                                                                              >  2:509944/1‑150 (MQ=255)
                        ggAGTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGcc                                                                                                                              <  2:693740/150‑1 (MQ=255)
                              cGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtca                                                                                                                        >  2:5206801/1‑150 (MQ=255)
                               gCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtcaa                                                                                                                       >  2:4680465/1‑150 (MQ=255)
                               gCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtcaa                                                                                                                       >  2:4381503/1‑150 (MQ=255)
                                cTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtcaat                                                                                                                      <  1:1007663/150‑1 (MQ=255)
                                      gAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAACAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATcc                                                                                                                <  2:4008151/150‑1 (MQ=255)
                                         cGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAg                                                                                                             >  1:553777/1‑150 (MQ=255)
                                          gTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGc                                                                                                            <  1:2036553/150‑1 (MQ=255)
                                          gTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGc                                                                                                            >  1:1037211/1‑150 (MQ=255)
                                             cTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCa                                                                                                         <  1:3505672/150‑1 (MQ=255)
                                              ttAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGtc                                                                                                                         >  2:3428723/1‑133 (MQ=255)
                                                  gcgcCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGcc                                                                                                    <  1:1692636/150‑1 (MQ=255)
                                                    gcCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGcccc                                                                                                  >  1:4611671/1‑150 (MQ=255)
                                                     ccGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCa                                                                                                 >  1:1294403/1‑150 (MQ=255)
                                                        gAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCAccc                                                                                              >  1:2524250/1‑150 (MQ=255)
                                                            aCGTTATCGCGATGTCCAAGAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCt                                                                                          <  1:5180075/150‑1 (MQ=255)
                                                            aCGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCt                                                                                          <  1:3001070/150‑1 (MQ=255)
                                                              gTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGa                                                                                        >  1:4858758/1‑150 (MQ=255)
                                                              gTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGa                                                                                        <  2:1282124/150‑1 (MQ=255)
                                                              gTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGa                                                                                        <  1:1247151/150‑1 (MQ=255)
                                                               ttATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGAt                                                                                       >  1:1930168/1‑150 (MQ=255)
                                                                tATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATa                                                                                      <  1:4576653/150‑1 (MQ=255)
                                                                        tGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATca                                                                              <  1:4724690/150‑1 (MQ=255)
                                                                          tCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATcaca                                                                            >  1:3929436/1‑150 (MQ=255)
                                                                          tCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATcaca                                                                            >  1:2196929/1‑150 (MQ=255)
                                                                            cAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcg                                                                          >  2:4795451/1‑150 (MQ=255)
                                                                            cAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcg                                                                          >  1:928860/1‑150 (MQ=255)
                                                                             aaaaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcgc                                                                         >  1:1377419/1‑150 (MQ=255)
                                                                             aaaaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcgc                                                                         >  2:1951722/1‑150 (MQ=255)
                                                                             aaaaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCACAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcgc                                                                         <  2:613629/150‑1 (MQ=255)
                                                                             aaaaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCACAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcgc                                                                         <  2:614273/150‑1 (MQ=255)
                                                                             aaaaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCAGTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcgc                                                                         >  1:1274262/1‑150 (MQ=255)
                                                                              aaaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACAcgcg                                                                        >  2:3230561/1‑150 (MQ=255)
                                                                               aaaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGt                                                                       >  2:768590/1‑150 (MQ=255)
                                                                                aaCAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAACGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGtc                                                                      <  1:3764838/150‑1 (MQ=255)
                                                                                   aGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCa                                                                   <  2:4874937/150‑1 (MQ=255)
                                                                                     gCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCaaa                                                                 >  2:1509348/1‑150 (MQ=255)
                                                                                     gCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCaaa                                                                 >  2:1052589/1‑150 (MQ=255)
                                                                                      cAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAaag                                                                <  1:5452872/150‑1 (MQ=255)
                                                                                       aaGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAaaga                                                               <  1:5361527/150‑1 (MQ=255)
                                                                                        aGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAaagaa                                                              <  1:1473736/150‑1 (MQ=255)
                                                                                        aGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAaagaa                                                              >  2:1725561/1‑150 (MQ=255)
                                                                                        aGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAaagaa                                                              <  1:447682/150‑1 (MQ=255)
                                                                                             ggCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGttct                                                         >  2:2612873/1‑150 (MQ=255)
                                                                                              gCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGttctt                                                        <  2:4957167/150‑1 (MQ=255)
                                                                                                  ccGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTcc                                                    <  2:3813848/150‑1 (MQ=255)
                                                                                                  ccGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTcc                                                    >  1:2607469/1‑150 (MQ=255)
                                                                                                   cGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTc                                                                                   >  1:370219/1‑118 (MQ=255)
                                                                                                   cGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCct                                                   <  2:5348127/150‑1 (MQ=255)
                                                                                                    gTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCctt                                                  <  1:1564429/150‑1 (MQ=255)
                                                                                                     tAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGttctt                                                        >  1:1515622/1‑143 (MQ=255)
                                                                                                     tAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCcttc                                                 <  1:5311207/150‑1 (MQ=255)
                                                                                                      aGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCcttct                                                >  1:1200341/1‑150 (MQ=255)
                                                                                                          aaCGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCaa                                            <  1:2724859/150‑1 (MQ=255)
                                                                                                            cGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGt                                          <  2:3566381/150‑1 (MQ=255)
                                                                                                              tCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTcc                                        >  1:4166945/1‑150 (MQ=255)
                                                                                                               ccTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCa                                       <  1:685564/150‑1 (MQ=255)
                                                                                                               ccTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCa                                       <  2:4711861/150‑1 (MQ=255)
                                                                                                                cTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCAt                                      <  1:5610028/150‑1 (MQ=255)
                                                                                                                 tGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATg                                     >  2:2903493/1‑150 (MQ=255)
                                                                                                                      tgcttctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGAtctc                                <  2:2070173/150‑1 (MQ=255)
                                                                                                                       gcttctgcCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCg                               <  2:2159607/150‑1 (MQ=255)
                                                                                                                               cAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATCCTCGATCACACGCGTCTCAAAGAAGTTCGTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCa                       >  2:1381238/1‑150 (MQ=255)
                                                                                                                               cAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCa                       >  1:2430564/1‑150 (MQ=255)
                                                                                                                                 ggCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATg                     >  2:2504447/1‑150 (MQ=255)
                                                                                                                                     aTAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGaa                 >  2:3019841/1‑150 (MQ=255)
                                                                                                                                        acaaGGCTGCCAAATGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCCTGATCTCGCTCATCCATGGGACCgg              >  2:1517416/3‑150 (MQ=255)
                                                                                                                                        acaaGGCTGCCAAATGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCGTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACgg              >  1:3393764/3‑150 (MQ=255)
                                                                                                                                         caaGGCTGCCAAATGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCaa                                            >  2:4326762/2‑119 (MQ=255)
                                                                                                                                         aaaGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACggg             <  2:899102/150‑1 (MQ=255)
                                                                                                                                         aaaGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACggg             >  1:4829464/1‑150 (MQ=255)
                                                                                                                                           aGGCTGCCAAATGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGtt           >  1:1331854/1‑150 (MQ=255)
                                                                                                                                            ggCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTg          <  2:1087291/150‑1 (MQ=255)
                                                                                                                                            ggCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTg          <  2:3760006/150‑1 (MQ=255)
                                                                                                                                             gCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTgg         >  1:2858348/1‑150 (MQ=255)
                                                                                                                                               tGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACg               <  2:558248/142‑1 (MQ=255)
                                                                                                                                               tGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACg               >  1:558248/1‑142 (MQ=255)
                                                                                                                                                gCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTAg      >  1:1970987/1‑150 (MQ=255)
                                                                                                                                                   aTTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTAGTcc   >  1:138243/1‑150 (MQ=255)
                                                                                                                                                   aTTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTAGTcc   >  1:3566353/1‑150 (MQ=255)
                                                                                                                                                    tttGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTAGTccc  >  1:3154999/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGAGAGCCGTTTCTGACGAGAACGGGATTCCGCTATGTGAGCGTTCTTAAGCGCCGGAAAACGTTATCGCGATGTCCAAAAACAGGCAAGGTCGGCAGCCGTAGGAAACGTCCTGTTCTGCTTCTGCCAGGCAATAAAAAGGCTGCCATTTGGCAGCCTTGTTGATGCTCAGCCCGGTCAATCAATCCCAGCTCAACGCCCCACCCGTCTGATACTCGATCACACGCGTCTCAAAGAAGTTCTTCTCCTTCTTCAAGTCCATGATCTCGCTCATCCATGGGAACGGGTTGGTAGTCCC  >  NC_002947/1355308‑1355605

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