Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,402,563 T→A 7.2% intergenic (+44/‑108) PP_3010 → / → PP_3011 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,402,5630TA 7.2% 262.7 / 9.6 83intergenic (+44/‑108)PP_3010/PP_3011hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/38);  new base A (4/2);  total (43/40)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.44e-01

GGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAGCAAGAAGCGTA  >  NC_002947/3402421‑3402712
                                                                                                                                              |                                                                                                                                                      
ggCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTtgtgt                                                                                                                                                 <  1:957022/150‑1 (MQ=255)
ggCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACAAAAAGCCCGCTCTCGACGGCGGGCTATTggtgt                                                                                                                                                 >  2:116311/1‑150 (MQ=38)
    cGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACAAAAAGCCCGCTCTCGACGGCGGGCTATTggtgtgggg                                                                                                                                             <  1:116311/150‑4 (MQ=37)
         agcgagcgCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGc                                                                                                                                        <  1:1794131/150‑1 (MQ=255)
               cgcCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGCCAGGTTGCAGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGccc                                                                                                                                  >  1:857384/1‑150 (MQ=255)
                  cACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGccc                                                                                                                                  >  2:2193982/1‑147 (MQ=255)
                  cACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGccc                                                                                                                                  <  1:2193982/147‑1 (MQ=255)
                      tCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACg                                                                                                                           <  2:2311813/150‑1 (MQ=255)
                       cACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCTCTCGACGGCGGGCTATTggtgtgggggt                                                                                                                                           <  2:1682880/133‑6 (MQ=37)
                       cACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCTCTCGACGGCGGGCTATTggtgtgggggt                                                                                                                                           >  1:1682880/1‑128 (MQ=37)
                               aTCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCaaaaaa                                                                                                                 >  2:1166826/1‑150 (MQ=255)
                               aTCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCaaaaaa                                                                                                                 >  2:1166779/1‑150 (MQ=255)
                                tCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCaaaaaaa                                                                                                                >  2:2162037/1‑150 (MQ=255)
                                tCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCaaaaaaa                                                                                                                >  1:1298591/1‑150 (MQ=255)
                                tCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCaaaaaaa                                                                                                                <  1:264259/150‑1 (MQ=255)
                                tCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATATGTCAACAACCGTAGTCATTCAACCCCCACATCACTAGGCAGCCGTCGCGAGCGGGCTTTTTGCGTGCCCGTCGTCTTGCAGGCCACGTTctacagca                                                                                                                 >  2:30535/1‑143 (MQ=255)
                                  aaCGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAgg                                                                                                              >  2:2123629/1‑150 (MQ=255)
                                           ccTGGTACGACCATAGCAGGTTGCTGGATGCAGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGACGTCGAGAGCGGGCTTATTGTGTGCCCATCGACTGCCCGGAAACATTCAAAAAAAGGGCGGGTGGa                                                                                                     <  2:28414/150‑1 (MQ=255)
                                            cTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAg                                                                                                    >  2:2262432/1‑150 (MQ=255)
                                               gTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCg                                                                                                 <  2:1888254/150‑1 (MQ=255)
                                                tACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCgg                                                                                                >  2:1528319/1‑150 (MQ=255)
                                                  cGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTATTCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTa                                                                                              <  1:147434/150‑1 (MQ=255)
                                                  cGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGTTa                                                                                              <  2:1298591/150‑1 (MQ=255)
                                                   gACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTAt                                                                                             >  2:432816/1‑150 (MQ=255)
                                                      cATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGAAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGc                                                                                          <  2:2127186/150‑1 (MQ=255)
                                                        tGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCt                                                                                        <  2:970819/150‑1 (MQ=255)
                                                        tGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCt                                                                                        >  2:963102/1‑150 (MQ=255)
                                                             ggTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGagca                                                                                   <  1:1166826/150‑1 (MQ=255)
                                                             ggTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGagca                                                                                   <  1:1166779/150‑1 (MQ=255)
                                                              gTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCg                                                                                                 <  2:1942439/135‑1 (MQ=255)
                                                              gTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCg                                                                                                 >  1:1942439/1‑135 (MQ=255)
                                                                 gCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCAATCGCCTGCCCGGCAACCTTCGCAAAACGGGCGGAAGCGTCCGGGATGGCCTGAGCAGCGt                                                                               >  2:22019/1‑150 (MQ=255)
                                                                  cTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGtt                                                                              <  2:662233/150‑1 (MQ=255)
                                                                  cTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCACAACCAGAATAACCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGtt                                                                              <  2:1270449/150‑1 (MQ=255)
                                                                    ggATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGc                                                                            >  2:647920/1‑150 (MQ=255)
                                                                     gATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCt                                                                           >  2:973289/1‑150 (MQ=255)
                                                                       tgctgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCATACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCttt                                                                         <  1:1893086/150‑1 (MQ=255)
                                                                       tgctgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTTCCCATCGCCTGCCCGGCAACGTTCAAAAAAAAGGTGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCttt                                                                         <  1:1769421/150‑1 (MQ=255)
                                                                        gctgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTg                                                                        <  1:963102/150‑1 (MQ=255)
                                                                          tgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCa                                                                      >  1:952974/1‑150 (MQ=255)
                                                                          tgcACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATAGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCa                                                                      >  1:603815/1‑150 (MQ=255)
                                                                              cTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATAGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTg                                                                  <  2:603815/150‑1 (MQ=255)
                                                                                 aTGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCg                                                               <  1:2123629/150‑1 (MQ=255)
                                                                                      ttCAACAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGcccc                                                          >  1:2306218/1‑150 (MQ=255)
                                                                                      ttCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGcccc                                                          >  1:2381418/1‑150 (MQ=255)
                                                                                      ttCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGcccc                                                          >  2:575857/1‑150 (MQ=255)
                                                                                         aaaaGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTg                                                                                                                                     >  2:361626/1‑73 (MQ=255)
                                                                                         aaaaGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTg                                                                                                                                     <  1:361626/73‑1 (MQ=255)
                                                                                         aaaaGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGc                                                                                                   <  1:2211840/106‑1 (MQ=255)
                                                                                         aaaaGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGc                                                                                                   >  2:2211840/1‑106 (MQ=255)
                                                                                         aaaaGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGTCGGCcccgcc                                                       >  1:2260420/1‑150 (MQ=255)
                                                                                          aaaGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCcccgccc                                                      <  1:432816/150‑1 (MQ=255)
                                                                                                      tcCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAAGGGCGGGTGGAGCCGGTATg                                                                                            >  2:1141493/1‑101 (MQ=255)
                                                                                                      tcCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGa                                                                                      <  1:2357147/106‑1 (MQ=255)
                                                                                                      tcCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGa                                                                                      >  2:2357147/1‑106 (MQ=255)
                                                                                                      gccAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAAGGGCGGGTGGAGCCGGTATg                                                                                            <  1:1141493/100‑1 (MQ=255)
                                                                                                         aaCCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCGTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGTCTTAGCAGCGTTGCTTTGCATTTTGCGGCCCCGCCCGTAATATTCGTTGTa                                       >  1:881619/1‑150 (MQ=255)
                                                                                                         aaCCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCa                                                                      >  1:2276919/1‑119 (MQ=255)
                                                                                                         aaCCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCa                                                                      <  2:2276919/119‑1 (MQ=255)
                                                                                                          atgggcacacaAAAAGCCCGCTCTCGACGGCGGGCTATTGGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTAGGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAt                                      >  2:8293/6‑150 (MQ=37)
                                                                                                            ccccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAttt                                    >  1:1563626/1‑150 (MQ=255)
                                                                                                            ccccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAttt                                    >  1:961219/1‑150 (MQ=255)
                                                                                                            ccccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAttt                                    <  1:2162037/150‑1 (MQ=255)
                                                                                                            ccccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAttt                                    <  1:647920/150‑1 (MQ=255)
                                                                                                            ccccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAttt                                    >  1:1402788/1‑150 (MQ=255)
                                                                                                              ccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCttt                                                                         >  2:585007/1‑111 (MQ=255)
                                                                                                              ccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCttt                                                                         <  1:585007/111‑1 (MQ=255)
                                                                                                              ccACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTGGCTTTGCAGTTTCCGGCCCCGCCCGTAATATTCGTTGTATTtgg                                  >  2:1476053/1‑150 (MQ=255)
                                                                                                               cacaCCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGTGGGGGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtgga                                 >  2:1288559/1‑150 (MQ=255)
                                                                                                                acacCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTAAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTTCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggat                                >  2:505518/1‑150 (MQ=255)
                                                                                                                  acaaaaaGCCCGCTCTCGACGGCGGGCTATTGGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTAt                                      >  2:1799798/7‑142 (MQ=38)
                                                                                                                  acCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTtggatgg                              >  1:836882/1‑150 (MQ=255)
                                                                                                                        aGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTc                        <  2:1555784/150‑1 (MQ=255)
                                                                                                                          cccGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGc                      <  1:2262432/150‑1 (MQ=255)
                                                                                                                          cccGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGc                      >  2:1652165/1‑150 (MQ=255)
                                                                                                                            cgccgTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGtt                                          <  2:663980/128‑1 (MQ=255)
                                                                                                                            cgccgTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGtt                                          >  1:663980/1‑128 (MQ=255)
                                                                                                                              ccgTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGt                  <  2:2306218/150‑1 (MQ=255)
                                                                                                                              ccgTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTGGGATGGAGAGTCGCAGGt                  <  1:973289/150‑1 (MQ=255)
                                                                                                                                  cGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGaa              <  1:1652165/150‑1 (MQ=255)
                                                                                                                                   gagaGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAg             >  1:100705/1‑150 (MQ=255)
                                                                                                                                      agCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAGCaa          <  1:1520609/150‑1 (MQ=255)
                                                                                                                                         gggCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAGCaagaa       <  1:200983/150‑1 (MQ=255)
                                                                                                                                             tttttGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGgag                            >  1:1726440/1‑125 (MQ=255)
                                                                                                                                             tttttGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGgag                            <  2:1726440/125‑1 (MQ=255)
                                                                                                                                             tttttGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAGCAAGAAGCGt   >  2:1318193/1‑150 (MQ=255)
                                                                                                                                              ttttGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAGCAAGAAGCGTa  <  2:683416/150‑1 (MQ=255)
                                                                                                                                              ttttGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAGCAAGAAGCGTa  <  2:1181993/150‑1 (MQ=255)
                                                                                                                                              ttttGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAGCAAGAAGCGTa  <  2:836882/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                      
GGCACGGGAAGCGAGCGCCACTTCACGGCCTATCAACGTGATGCCTGGTACGACCATGGCAGGTTGCTGGATGCTGCACTGATGTGTTCAAAAGCGTAATCATCCAACCCCCACACCAATAGCCCGCCGTCGAGAGCGGGCTTTTTGTGTGCCCATCGCCTGCCCGGCAACGTTCAAAAAAAGGGCGGGTGGAGCCGGTATGGCCTGAGCAGCGTTGCTTTGCATTTGCCGGCCCCGCCCGTAATATTCGTTGTATTTGGATGGAGAGTCGCAGGTGGAAGCAAGAAGCGTA  >  NC_002947/3402421‑3402712

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