Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 707,055 C→T 8.4% intergenic (+231/+15) PP_0599 → / ← rpsT paraquat‑inducible protein B/30S ribosomal protein S20

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947707,0550CT 8.4% 164.8 / 4.2 59intergenic (+231/+15)PP_0599/rpsTparaquat‑inducible protein B/30S ribosomal protein S20
Reads supporting (aligned to +/- strand):  ref base C (26/28);  new base T (3/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.54e-01

AGCGCAGAAAGTGAGATTGCCGGGGCTGCCCCGCAGCCCTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGGATGATACCCTTGTCGGCCATACGGTCGATTACAGGTACAGCCAGAACGTAAGCG  >  NC_002947/706914‑707202
                                                                                                                                             |                                                                                                                                                   
aGCGCAGAAAGTGAGATTGCCGGGGCTGCCCCGCAGCCCTATCGCGACACAACGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGAACCTAGGGTCGCTTTTTCGACTAgc                                                                                                                                              >  8:230019/1‑149 (MQ=255)
aGCGCAGAAACTGAGATTGCCGGGGCTGCCCCGCAGCCCTAGCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCACAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAgc                                                                                                                                              <  3:6539/149‑1 (MQ=255)
       aaaGTGAGATTGCCGGGGCTGCCCCGCAGCCCTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTc                                                                                                                                                     <  5:32841/135‑1 (MQ=255)
       aaaGTGAGATTGCCGGGGCTGCCCCGCAGCCCTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTc                                                                                                                                                     >  6:32841/1‑135 (MQ=255)
           tGAGATTGCCGGGGCTGCCCCGCAGCCCTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAg                                                                                                                                   >  4:94769/1‑149 (MQ=255)
              gaTTGCCGGGGCTGCCCCGCAGCCCTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGG‑TTTTCGACTAGCGCGTCGCTTAGgcag                                                                                                                               <  5:242585/149‑1 (MQ=255)
                   ccGGGGCTGCCCCGCAGCCCTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGgcagctgc                                                                                                                           <  7:230019/149‑1 (MQ=255)
                   ccGGGGCTGCCCCGCAGCCCTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGgcagctg                                                                                                                            >  2:179687/1‑148 (MQ=255)
                    cGGGGCTGCCCCGCAGCCCTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGgcagctgca                                                                                                                          <  3:10691/149‑1 (MQ=255)
                    cGGGGCTGCCCCGCAGCCCTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGgcagctg                                                                                                                            >  8:371736/1‑147 (MQ=255)
                          tGCCCCGCAGACCTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTACAGcttc                                                                                                                     >  7:231246/1‑148 (MQ=255)
                            ccccGCAGCCCTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCttt                                                                                                                  >  2:344858/1‑149 (MQ=255)
                            ccccGCAGCCCTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCttt                                                                                                                  >  5:168565/1‑149 (MQ=255)
                            ccccGCAGCCCTATCGCGACACAAGGCCGCTCCCACAGATCGACTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTACAGcttc                                                                                                                     <  8:231246/146‑1 (MQ=255)
                              ccGCAGCCCTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCa                                                                                                                <  5:54938/149‑1 (MQ=255)
                                gCAGCCCTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAgc                                                                                                              <  7:167370/149‑1 (MQ=255)
                                    ccatATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGAACCAAGGGTCGGTTTTTCGACTAGCGCGTCGATTAGGCAGCTGCAGCTTCTTTCAGCAcca                                                                                                          >  2:355039/4‑148 (MQ=255)
                                      cTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGTAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGgca                                                                                                                                <  4:34812/125‑1 (MQ=255)
                                      cTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGgca                                                                                                                                >  3:34812/1‑125 (MQ=255)
                                         tCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGAtg                                                                                                     >  7:379635/1‑149 (MQ=255)
                                           gcgACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGAtgtg                                                                                                   <  1:179687/149‑1 (MQ=255)
                                           gcgACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGAtgtg                                                                                                   >  2:104742/1‑149 (MQ=255)
                                            cgACACAAGGCCGCTCCCACAGATCGAGTCCCGCAAGGAACTGGGCCCGCAAGAATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTAATTTGACAAGCGCGTCGCTTAGGCAGATCCAGCTTCATTCAGCGCCTTTATGTgg                                                                                                  >  5:234370/1‑149 (MQ=255)
                                                         cTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGc                                                                                     >  3:365477/1‑149 (MQ=255)
                                                          tCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGcc                                                                                    >  2:24951/1‑149 (MQ=255)
                                                            ccACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCtt                                                                                  >  6:285065/1‑149 (MQ=255)
                                                            ccACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCtt                                                                                  >  8:186662/1‑149 (MQ=255)
                                                             gtcaGTTCGAGTCCCGTATGGTTCTGGGCCCGCAAGCATATATCGCAGGCAACATTAATCCGATCCTAGGGGCGGTTTTTCGACTAGCGCGTCGGTTAGGCAGCTGCAGGTTCTTTCAGCGCCTTGATGTGGCGATTCAGACGGCCttt                                                                                 <  7:300930/147‑1 (MQ=255)
                                                             cacaGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCTACAAAAAACCTACCCTAGGGTAGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCttt                                                                                 <  1:355039/149‑1 (MQ=255)
                                                             cacaGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCttt                                                                                 <  6:168565/149‑1 (MQ=255)
                                                              acaGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTtg                                                                                <  8:73062/149‑1 (MQ=255)
                                                               gaGATCGAGTCGGGCATGGAACTTGTCCCGCAATCATAAATCGCAGGAAACAAAAAACCGACCCTTGGGTCGGTTTTTCGACTTGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTtgt                                                                               <  5:344280/148‑1 (MQ=255)
                                                                 gATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGa                                                                             <  8:379635/149‑1 (MQ=255)
                                                                  ttCGAGTGCCGCATGGTACTGGTCCCGCAAGGATATATCGCTGGCAACAAAAATACGACCCTAGGGTCTGTTTTTTGACTAGGGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTTATGTGGCCATTCAGACGGCCTTTGTGAc                                                                            <  6:234370/148‑1 (MQ=255)
                                                                     gAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGag                                                                         <  5:285065/149‑1 (MQ=255)
                                                                        tCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGagcag                                                                      <  7:371736/149‑1 (MQ=255)
                                                                          ccGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGcc                                                                    >  4:34257/1‑149 (MQ=255)
                                                                            gCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCtt                                                                  <  5:101407/149‑1 (MQ=255)
                                                                            gCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCtt                                                                  <  1:104742/149‑1 (MQ=255)
                                                                                ggAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTc                                                              >  2:377548/1‑149 (MQ=255)
                                                                                 gAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCt                                                             >  1:57683/1‑149 (MQ=255)
                                                                                  aaCTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCtt                                                            >  3:320592/1‑149 (MQ=255)
                                                                                   aCTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTtg                                                           >  2:20818/1‑149 (MQ=255)
                                                                                       ggCCGGCAAGCATAAATCGCAGGCAACAAAAATCCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGga                                                       <  2:354701/149‑1 (MQ=255)
                                                                                       ggCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTtgtg                                                         <  1:377548/147‑1 (MQ=255)
                                                                                       ggCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGga                                                       <  3:94769/149‑1 (MQ=255)
                                                                                          ccGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTATTGTTTCGCTTAGGTTGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGgatga                                                    <  1:334003/149‑1 (MQ=255)
                                                                                                       tCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGGATGATACCCTTGTCGGc                                       >  3:97119/1‑149 (MQ=255)
                                                                                                         gactgccacaaACACCCGACCCTAGGGGCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTTATGTGGCCATTCAGACGGCCTTTGTGCCGAGCAGCCTTGTTCTTGTGGATGATACCCTTGTCGGc                                       <  2:331157/140‑1 (MQ=255)
                                                                                                         gCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGGATGATACCCTTGTCGGcc                                      <  4:365477/148‑1 (MQ=255)
                                                                                                                tcgaaaaaCCGACCCTAGGGTCGGTTTTTTGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGGATGATACCCTTGTCGGCCATa                                   <  6:174724/141‑1 (MQ=255)
                                                                                                                tcgaaaaaCCGACCCTAGGGTCGGTTTTTTGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGGATGATACCCTTGTCGGCCATAcagat                              >  5:174724/4‑145 (MQ=255)
                                                                                                                  gaaaaaCCGACCCTAGGGTCGGTTTTTTGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGGATGATACCCTTGTCGGCCATACGGTCGa                            >  4:326646/2‑149 (MQ=255)
                                                                                                                  aaaaaaCCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGGATGATACCCTTGTCGGCCATACGGTCGa                            >  1:368263/1‑149 (MQ=255)
                                                                                                                   aaaaaCCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGGATGATACCCTTGTCGGCCATACGGTCGAt                           >  4:119069/1‑149 (MQ=255)
                                                                                                                     aaaCCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGGATGATACCCTTGTCGGCCATACGGTCGATTa                         >  8:320751/1‑149 (MQ=255)
                                                                                                                       aCCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGGATGATACCCTTGTCGGCCATACGGTCGATTACa                       <  2:368263/149‑1 (MQ=255)
                                                                                                                               tAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGGATGATACCCTTGTCGGCCATACGGTCGa                            >  6:82077/1‑136 (MQ=255)
                                                                                                                               tAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGGATGATACCCTTGTCGGCCATACGGTCGa                            <  5:82077/136‑1 (MQ=255)
                                                                                                                                   ggcGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGGATGATACCCTTGTCGGCCATACGGTCGATTACAGGTACAGCCAGa           <  4:320592/147‑1 (MQ=255)
                                                                                                                                            tCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGGATGATACCCTTGTCGGCCATACGGTCGATTACAGGTACAGCCAGAACGTAAGCg  >  5:203208/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
AGCGCAGAAAGTGAGATTGCCGGGGCTGCCCCGCAGCCCTATCGCGACACAAGGCCGCTCCCACAGATCGAGTCCCGCATGGAACTGGGCCCGCAAGCATAAATCGCAGGCAACAAAAAACCGACCCTAGGGTCGGTTTTTCGACTAGCGCGTCGCTTAGGCAGCTGCAGCTTCTTTCAGCGCCTTGATGTGGCCATTCAGACGGCCTTTGTGACGAGCAGCCTTGTTCTTGTGGATGATACCCTTGTCGGCCATACGGTCGATTACAGGTACAGCCAGAACGTAAGCG  >  NC_002947/706914‑707202

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