Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 142,889 C→T 6.4% intergenic (+11/‑164) PP_0134 → / → PP_5428 hypothetical protein/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947142,8890CT 6.4% 227.5 / 3.0 78intergenic (+11/‑164)PP_0134/PP_5428hypothetical protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (33/39);  new base T (3/2);  total (36/42)
Fisher's exact test for biased strand distribution p-value = 6.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

ATGGTGCGGCTCTCGGGCGAGGTGGTCAGCGCGGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTACCTGTTCAGGTTTGAACAAAAAATAC  >  NC_002947/142747‑143036
                                                                                                                                                |                                                                                                                                                   
aTGGTGCGGCTCTCGGGCGAGGTGGTCAGCGCGGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGa                                                                                                                                               <  2:271170/149‑1 (MQ=255)
       ggCTCTCGGGCGAGGTGGTCAGCGCGGAGCGAGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACcaca                                                                                                                                          >  1:244931/1‑148 (MQ=255)
       ggCTCTCGGGCGAGGTGGTCAGCGCGGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACcaca                                                                                                                                          >  3:163047/1‑147 (MQ=255)
        gCTCTCGGGCGAGGTGGTCAGCGCGGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAg                                                                                                                                       >  8:146135/1‑149 (MQ=255)
           ctcGGGCGAGGTGGTCAGTGCGGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCATTATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGc                                                                                                                                     >  1:116972/1‑149 (MQ=255)
            tcGGGCGAGGTGGTCAGCGCGGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGc                                                                                                                                     <  2:244931/147‑1 (MQ=255)
                gCGAGGTGGTCAGCGCGGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATc                                                                                                                               >  6:31456/1‑149 (MQ=255)
                   aGGTGGTCAGCGCGGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTg                                                                                                                            <  4:380547/149‑1 (MQ=255)
                     gtggtCAGCGCGGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTg                                                                                                                            >  6:395392/1‑147 (MQ=255)
                      tggtCAGCGCGGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCGAATTGCACGATcag                                                                                                                                           >  7:201678/1‑130 (MQ=255)
                      tggtCAGCGCGGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCGAATTGCACGATcag                                                                                                                                           <  8:201678/131‑2 (MQ=255)
                        gtCAGCGCGGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTg                                                                                                                            <  7:43698/144‑1 (MQ=255)
                        gtCAGCGCGGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAgg                                                                                                                        <  5:40469/148‑1 (MQ=255)
                        gtCAGCGCGGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGCCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTg                                                                                                                            >  8:43698/1‑144 (MQ=255)
                         tCAGTGCGGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGAt                                                                                                                      <  2:116972/149‑1 (MQ=255)
                               cgGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCAc                                                                                                                                                 <  4:179033/116‑1 (MQ=255)
                               cgGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCAc                                                                                                                                                 >  3:179033/1‑116 (MQ=255)
                                    c‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTcaca                                                                                                           <  7:146135/149‑1 (MQ=255)
                                    c‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTcaca                                                                                                           <  5:31456/149‑1 (MQ=255)
                                            ccATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGtttt                                                                                                    >  1:104749/1‑149 (MQ=255)
                                               tCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTta                                                                                                 <  1:249146/149‑1 (MQ=255)
                                                  ccGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTtaata                                                                                              <  1:205356/149‑1 (MQ=255)
                                                   cGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTtaataa                                                                                             <  8:374226/149‑1 (MQ=255)
                                                      ggAGATCGCCCGGCA‑AATCATCGGCGTACGCGGCGTGGAAGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCCTCCAGGATACGTCCTCAAAGAGTTTTTTAATAATCa                                                                                          >  2:363528/1‑149 (MQ=255)
                                                       gagaTCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCac                                                                                         <  4:371977/149‑1 (MQ=255)
                                                          aTCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGt                                                                                                                                   >  3:305888/1‑104 (MQ=255)
                                                          aTCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGt                                                                                                                                   <  4:305888/104‑1 (MQ=255)
                                                              cccGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTGTTTAATAATCACTAACTAg                                                                                  >  6:242677/1‑149 (MQ=255)
                                                                 ggCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGtt                                                                                >  6:367893/1‑148 (MQ=255)
                                                                 ggCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTg                                                                               <  8:60914/149‑1 (MQ=255)
                                                                   cA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTGATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCAGTAACTAGTTGAt                                                                             <  4:117545/149‑1 (MQ=255)
                                                                    a‑aaTCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTcaca                                                                                                           >  7:285880/1‑118 (MQ=255)
                                                                    a‑aaTCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTcaca                                                                                                           <  8:285880/118‑1 (MQ=255)
                                                                    a‑aaTCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTc                                                                                                              <  1:28796/115‑1 (MQ=255)
                                                                    a‑aaTCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTc                                                                                                              >  2:28796/1‑115 (MQ=255)
                                                                    a‑aaTCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGAtt                                                                            >  6:174690/1‑149 (MQ=255)
                                                                    a‑aaTCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGAtt                                                                            <  1:113765/149‑1 (MQ=255)
                                                                         catcGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTtat                                                                         >  1:156690/1‑148 (MQ=255)
                                                                          atcGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCAGACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATaa                                                                       <  2:278295/149‑1 (MQ=255)
                                                                              gCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGAt                                                                   <  5:367893/149‑1 (MQ=255)
                                                                                   cGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGTATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATtttat                                                              <  1:245367/149‑1 (MQ=255)
                                                                                    ggggCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATtttatt                                                             <  5:174690/149‑1 (MQ=255)
                                                                                      ggCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATtttattt                                                            >  6:144062/1‑148 (MQ=255)
                                                                                         gTGGATGCTGACCTGCTGAAAGTCGCGACCATGGTGGAAGGCTGATGGTGCAATTGGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCAGAgtgt                                                                                                       <  5:330582/102‑1 (MQ=255)
                                                                                         gTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACACGCGTCATCGTGCAGGATACGTCCTCACAgtgt                                                                                                       >  6:330582/1‑102 (MQ=255)
                                                                                          tGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATc                                                       <  5:395392/149‑1 (MQ=255)
                                                                                              tGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCaaaa                                                   <  4:197290/149‑1 (MQ=255)
                                                                                               gCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGttt                                                                                                     <  5:178025/98‑1 (MQ=255)
                                                                                               gCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGttt                                                                                                     >  6:178025/1‑98 (MQ=255)
                                                                                                   aactgcttAAAGTCGCGACCATGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGTTTCGTCCTCACAGTGTTTTTTAATAATGACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTg                                              <  1:363528/147‑1 (MQ=255)
                                                                                                   aCCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTg                                              >  5:401867/1‑149 (MQ=255)
                                                                                                    ccTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCaaa                                                    <  1:309139/142‑1 (MQ=255)
                                                                                                    ccTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCaaa                                                    >  2:309139/1‑142 (MQ=255)
                                                                                                     ctgctgAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTgg                                             >  5:14628/1‑148 (MQ=255)
                                                                                                     ctgctgAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGc                                            <  5:144062/149‑1 (MQ=255)
                                                                                                       gctgAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCAt                                          >  3:99584/1‑149 (MQ=255)
                                                                                                          gAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATga                                        >  3:243669/1‑147 (MQ=255)
                                                                                                          gAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGc                                        <  4:243669/148‑1 (MQ=255)
                                                                                                              gTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGa                                                                              >  4:246463/1‑106 (MQ=255)
                                                                                                              gTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCt                                   <  6:64949/149‑1 (MQ=255)
                                                                                                                   gACCAAGGTTGAAGGCTGATCGTGCAATTCGCACGAACACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGAt                                                                             <  7:142947/102‑1 (MQ=255)
                                                                                                                   gACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATAGTGCAGGATACGTCCTCACAGTCTTTTTTAATAAACACTAACTAGTAGAt                                                                             >  8:142947/1‑102 (MQ=255)
                                                                                                                     ccAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTa                            >  5:45986/1‑149 (MQ=255)
                                                                                                                      cAAGGTGGAAGGCTGATCGTGCGATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAAt                                                        >  7:261051/1‑120 (MQ=255)
                                                                                                                      cAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGAtttt                                                                >  7:116844/1‑112 (MQ=255)
                                                                                                                      cAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGAtttt                                                                <  8:116844/112‑1 (MQ=255)
                                                                                                                      cAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAAt                                                        <  8:261051/120‑1 (MQ=255)
                                                                                                                        aGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCaaa                                                    >  3:371552/1‑122 (MQ=255)
                                                                                                                        aGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCAACGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCaaa                                                    <  4:371552/122‑1 (MQ=255)
                                                                                                                          gTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTACCTGt                       <  7:8481/149‑1 (MQ=255)
                                                                                                                                    gATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGGTACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTACCTGTTCAGGTTTg              >  3:183304/1‑148 (MQ=255)
                                                                                                                                     gtCGTGCGAATTGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTACCTGTTCAGGTTTGa             >  8:38679/2‑148 (MQ=255)
                                                                                                                                     aTCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTACCTGTTCAGGTTTGa             <  2:104749/148‑1 (MQ=255)
                                                                                                                                      tCGTGCGAATTGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTg                                  >  8:123297/1‑126 (MQ=255)
                                                                                                                                      tCGTGCGAATTGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTg                                  <  7:123297/126‑1 (MQ=255)
                                                                                                                                           cAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTgg                                             <  1:340022/110‑1 (MQ=255)
                                                                                                                                           cAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAAATAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTgg                                             >  2:340022/1‑110 (MQ=255)
                                                                                                                                            aaTTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTACCTGTTCAGGTTTGAACaaaaaa     <  6:401867/149‑1 (MQ=255)
                                                                                                                                                cGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTACCTGTTCAGGTTTGAACAAAAAATAc  <  6:14628/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                   
ATGGTGCGGCTCTCGGGCGAGGTGGTCAGCGCGGAGC‑AGAAAACCATCGCCGTGGAGATCGCCCGGCA‑AATCATCGGCGTACGGGGCGTGGATGCTGACCTGCTGAAAGTCGCGACCAAGGTGGAAGGCTGATCGTGCAATTCGCACGACCACACAGGCGTCATCGTGCAGGATACGTCCTCACAGTGTTTTTTAATAATCACTAACTAGTTGATTTATAAGGATTTTATTTAAATCAAAAGGCTGGCATGCAGGCTGCAATTACCTGTTCAGGTTTGAACAAAAAATAC  >  NC_002947/142747‑143036

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