Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 557,339:1 +A 8.4% intergenic (+79/‑130) rpsN → / → rpsH 30S ribosomal protein S14/30S ribosomal protein S8

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947557,3391.A 8.4% 264.5 / 18.0 71intergenic (+79/‑130)rpsN/rpsH30S ribosomal protein S14/30S ribosomal protein S8
Reads supporting (aligned to +/- strand):  ref base . (31/34);  new base A (2/4);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 6.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

AACATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAGGA  >  NC_002947/557198‑557482
                                                                                                                                              |                                                                                                                                               
aaCATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCATATAGGGGCTTg                                                                                                                                         <  1:1094137/150‑1 (MQ=255)
aaCATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGACTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCATATAGGGGCTTg                                                                                                                                         <  1:793171/150‑1 (MQ=255)
       tgcGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTgtttcgtt                                                                                                                                  >  2:1388870/1‑150 (MQ=255)
               gCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCtt                                                                                                                                          <  1:51928/134‑1 (MQ=255)
               gCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCtt                                                                                                                                          >  2:51928/1‑134 (MQ=255)
                 tGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTatat                                                                                                                                                  <  1:644609/125‑1 (MQ=255)
                 tGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTatat                                                                                                                                                  >  2:644609/1‑125 (MQ=255)
                    aaTGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTc                                                                                                                     >  2:1096246/1‑150 (MQ=255)
                    aaTGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTc                                                                                                                     <  1:201545/150‑1 (MQ=255)
                             tGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGAc                                                                                                            <  2:879098/150‑1 (MQ=255)
                                 gTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGc                                                                                                        >  2:413658/1‑150 (MQ=255)
                                    ccAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCa                                                                                                     >  1:1005663/1‑150 (MQ=255)
                                         tCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGa                                                                                                <  2:1269791/150‑1 (MQ=255)
                                             gTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTgtttcgttt                                                                                                                                 <  1:1416046/113‑1 (MQ=255)
                                             gTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTgtttcgttt                                                                                                                                 >  2:1416046/1‑113 (MQ=255)
                                              tCAAGGCCTCCTGGTAATCGATTTAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGCCCGGCTCAGGTGAGccgg                                                                                           <  1:73824/150‑3 (MQ=255)
                                                 aGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGt                                                                                        >  2:843797/1‑150 (MQ=255)
                                                  ggCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGtt                                                                                       <  2:482477/150‑1 (MQ=255)
                                                     ctcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGttttt                                                                                    >  2:370884/1‑150 (MQ=255)
                                                     ctcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGttttt                                                                                    <  2:1126597/150‑1 (MQ=255)
                                                     ctcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGttttt                                                                                    <  2:1011391/150‑1 (MQ=255)
                                                     ctcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGttttt                                                                                    <  2:1005663/150‑1 (MQ=255)
                                                      tcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGttttt                                                                                    <  1:642185/149‑1 (MQ=255)
                                                      tcctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGttttt                                                                                    >  2:642185/1‑149 (MQ=255)
                                                        ctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCAGGTGAGccgg                                                                                           <  2:259235/140‑3 (MQ=255)
                                                        ctGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCAGGTGAGccgg                                                                                           >  1:259235/1‑138 (MQ=255)
                                                         tGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGccc                                                                                >  1:974686/1‑150 (MQ=255)
                                                          ggTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCa                                                                               >  2:1120365/1‑150 (MQ=255)
                                                           gTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTgg                                                                                          >  2:551712/1‑138 (MQ=255)
                                                           gTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTgg                                                                                          <  1:551712/138‑1 (MQ=255)
                                                            tAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCaaa                                                                             <  2:998293/150‑1 (MQ=255)
                                                              aTCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGc                                                                                              >  2:316822/1‑131 (MQ=255)
                                                              aTCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGc                                                                                              <  1:316822/131‑1 (MQ=255)
                                                                 gCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGtt                                                                        <  1:413658/150‑1 (MQ=255)
                                                                          aaGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGtt                                                               <  1:1082623/150‑1 (MQ=255)
                                                                            ggCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGt                                                             <  1:843797/150‑1 (MQ=255)
                                                                                      gAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCAGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCa                                                   >  1:530834/1‑150 (MQ=255)
                                                                                       aaGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCaa                                                  >  1:213127/1‑150 (MQ=255)
                                                                                       aaGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCaa                                                  >  2:552211/1‑150 (MQ=255)
                                                                                       aaGGGGAGCAATCTTAAGACCGTTGTTGCCCTAACCAACAAACAAGCCCCTATCTGGAGTTGCAGGCAGAATGGTTTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCaa                                                  <  1:1096246/150‑1 (MQ=255)
                                                                                         ggggAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAgg                                                >  1:1172508/1‑150 (MQ=255)
                                                                                           ggAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGc                                                                                                                             >  1:138875/1‑71 (MQ=255)
                                                                                           ggAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGc                                                                                                                             <  2:138875/71‑1 (MQ=255)
                                                                                           ggAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTc                                                                       <  2:61812/125‑1 (MQ=255)
                                                                                           ggAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTc                                                                       >  1:61812/1‑125 (MQ=255)
                                                                                            gAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCt                                                                                                  >  2:451040/1‑97 (MQ=255)
                                                                                            gAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCt                                                                                                  <  1:451040/97‑1 (MQ=255)
                                                                                            gAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTgg                                                                                          <  2:13070/105‑1 (MQ=255)
                                                                                            gAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTgg                                                                                          >  1:13070/1‑105 (MQ=255)
                                                                                             aGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCtaa                                            <  1:496726/150‑1 (MQ=255)
                                                                                             aGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCtaa                                            >  2:23083/1‑150 (MQ=255)
                                                                                                aaTCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCtaataa                                         <  1:370884/150‑1 (MQ=255)
                                                                                                    ttCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCAGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCt                                     <  2:530834/150‑1 (MQ=255)
                                                                                                     tCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCtt                                    <  1:1388870/150‑1 (MQ=255)
                                                                                                      ccGACCGCCGGTGACCTTGCCAACAAACAAGCCCCAAACGTGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTg                                   <  1:75257/150‑1 (MQ=255)
                                                                                                        gACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTa                                 >  2:604604/1‑150 (MQ=255)
                                                                                                              cgGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGCCTAGTGTCCCGGTCGCACCTGAGCCTGGGTTTTTGCCCACACGTTTTGTGCGGTCGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGa                           <  1:324211/150‑1 (MQ=255)
                                                                                                                   aCCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATc                      <  1:1199580/150‑1 (MQ=255)
                                                                                                                   aCCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATc                      <  1:23083/150‑1 (MQ=255)
                                                                                                                     cTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTa                    >  1:1221556/1‑150 (MQ=255)
                                                                                                                         ccAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGccc                >  1:15110/1‑150 (MQ=255)
                                                                                                                           aacaaacaAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCAt              >  2:681946/1‑150 (MQ=255)
                                                                                                                           aacaaacaAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCAt              >  2:810095/1‑150 (MQ=255)
                                                                                                                            cgaaacaaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCTCAAGGCTAATAATTCt                                     <  2:350079/124‑1 (MQ=255)
                                                                                                                            cgaaacaaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCt                                     >  1:350079/3‑126 (MQ=255)
                                                                                                                                  aaGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCt                                     >  1:144240/1‑120 (MQ=255)
                                                                                                                                  aaGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCt                                     <  2:144240/120‑1 (MQ=255)
                                                                                                                                  aaGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGa            >  1:1132669/1‑145 (MQ=255)
                                                                                                                                  aaGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGa            <  2:1132669/145‑1 (MQ=255)
                                                                                                                                  aaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCtt                                    >  2:1047401/1‑121 (MQ=255)
                                                                                                                                  aaGCCCCATATAGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCtt                                    <  1:1047401/121‑1 (MQ=255)
                                                                                                                                      cccTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAgg   >  2:1072375/1‑150 (MQ=255)
                                                                                                                                       ccTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAGGa  <  2:506155/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                               
AACATGCTGCGTCAAGCTGCAATGCGCGGTGACGTACCAGGTCTGGTCAAGGCCTCCTGGTAATCGCTGCAGGCAAGGCGCGGGTGGAAGGGGAGCAATCTTCCGACCGCCGGTGACCTTGCCAACAAACAAGCCCCTATATGGGGCTTGTTTCGTTTCTGCCTTGTGTCTAGAATGACCGGCTCACCTGAGCCTGGGTTTTTGCCCAAACGGTTCGGGTGGTTGTCATAGCCGCAAGGCTAATAATTCTTGTATCAGGAGCATCTAGCCCATGAGTATGCAGGA  >  NC_002947/557198‑557482

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