Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 289,942 C→A 7.2% intergenic (+26/‑140) oprE → / → PP_0235 outer‑membrane porin E/peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947289,9420CA 7.2% 278.0 / 8.9 92intergenic (+26/‑140)oprE/PP_0235outer‑membrane porin E/peroxidase
Reads supporting (aligned to +/- strand):  ref base C (41/44);  new base A (4/3);  total (45/47)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.10e-01

AAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATG  >  NC_002947/289795‑290084
                                                                                                                                                   |                                                                                                                                              
aaGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCt                                                                                                                                               >  5:348013/1‑149 (MQ=255)
 aGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTg                                                                                                                                              <  2:303831/149‑1 (MQ=255)
   gCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGAt                                                                                                                                            >  1:42442/1‑149 (MQ=255)
     aCTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTGCACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGc                                                                                                                                          >  2:224288/1‑149 (MQ=255)
       tctcAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGCCCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGc                                                                                                                                          >  6:203401/1‑147 (MQ=255)
         tcAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAgg                                                                                                                                      >  7:337839/1‑149 (MQ=255)
         tcAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAgg                                                                                                                                      >  2:139609/1‑149 (MQ=255)
           aaaGGCGTGGTCTTCAGCCTGCGTCAAGGTTTCTACCGTGTTAGCGGCGAAAGCGCCCTGTACAAAGACCATACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAg                                                                                                                                       <  4:381210/146‑1 (MQ=255)
           aaaGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCt                                                                                                                                    >  2:85555/1‑149 (MQ=255)
             aGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCATCAGGCGAGGCttc                                                                                                                                  <  5:181048/149‑2 (MQ=255)
              ggCCTGGGCTTCAGCCTGCGGCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCATCAGGCGAGGCttct                                                                                                                                 <  5:189327/149‑3 (MQ=255)
                ccTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttttt                                                                                                                               <  6:260136/149‑1 (MQ=255)
                    ggCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGc                                                                                                                           <  1:214971/149‑1 (MQ=255)
                    ggCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGc                                                                                                                           <  7:81062/149‑1 (MQ=255)
                               gCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTta                                                                                                                >  8:353080/1‑149 (MQ=255)
                                cGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAATAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTta                                                                                                                <  1:81961/148‑1 (MQ=255)
                                       gtgtCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttttt                                                                                                                               >  8:153700/1‑126 (MQ=255)
                                       gtgtCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttttt                                                                                                                               <  7:153700/126‑1 (MQ=255)
                                       gtgtCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTAtt                                                                                                        <  7:164471/149‑1 (MQ=255)
                                            tACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGt                                                                                                   >  6:323819/1‑149 (MQ=255)
                                             aCCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTATCTGAAAAAGAATCCTCGACTGATCCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGttt                                                                                                  >  2:149692/1‑149 (MQ=255)
                                              ccGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTgg                                                                                                    <  1:5245/146‑1 (MQ=255)
                                               cGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTtat                                                                                                <  5:323856/149‑1 (MQ=255)
                                                      gcggcgAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTtat                                                                                                               <  3:50946/127‑1 (MQ=255)
                                                      gcggcgAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTtat                                                                                                               >  4:50946/1‑127 (MQ=255)
                                                            aaaGCGCCCTGGACAAGGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGtta                                                                                   >  1:318894/1‑149 (MQ=255)
                                                             aaGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGGAGATTTATATGTATTTTGGTTATAAATAAATTGttat                                                                                  <  1:70668/149‑1 (MQ=255)
                                                                cgcCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGttattta                                                                               >  4:238942/1‑149 (MQ=255)
                                                                   ccTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTc                                                                            >  5:212281/1‑149 (MQ=255)
                                                                   ccTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTc                                                                            >  2:137673/1‑149 (MQ=255)
                                                                       gACAAAGACCAGACTCGCTTCATCGTGAACTGCACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTCGTTATAAATAAATTGTTATTTATTCttt                                                                         <  1:224288/148‑1 (MQ=255)
                                                                       gACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCTAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCttt                                                                         <  5:203401/148‑1 (MQ=255)
                                                                       gACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCttt                                                                         <  2:320144/148‑1 (MQ=255)
                                                                        aCAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCttttt                                                                       <  6:148163/149‑1 (MQ=255)
                                                                         cAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATaaataaat                                                                                        <  1:363037/131‑1 (MQ=255)
                                                                         cAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATaaataaat                                                                                        >  2:363037/1‑131 (MQ=255)
                                                                          aaaGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGt                                                                     >  3:61358/1‑149 (MQ=255)
                                                                          aaaGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGt                                                                     >  2:125205/1‑149 (MQ=255)
                                                                          aaaGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGt                                                                     >  4:287357/1‑149 (MQ=255)
                                                                             gACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGtttt                                                                  <  5:323819/149‑1 (MQ=255)
                                                                                 aGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAAt                                                              <  3:238942/149‑1 (MQ=255)
                                                                                 aGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAAt                                                              <  8:112445/149‑1 (MQ=255)
                                                                                  gACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATc                                                             <  1:85555/149‑1 (MQ=255)
                                                                                    cTCGCTTCATCGTGAACAACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTTTTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTGTTTAATcca                                                           >  4:22892/1‑147 (MQ=255)
                                                                                     tCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGAt                                                          >  5:186632/1‑149 (MQ=255)
                                                                                       gcatCATCGTGACCTACACTTACCACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTg                                                        <  8:253552/146‑1 (MQ=255)
                                                                                         ttCATCGTGAACTACACTTAAAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGAtt                                                         <  1:4437/146‑1 (MQ=255)
                                                                                           cATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAg                                                     >  2:223182/1‑148 (MQ=255)
                                                                                           cATCGTGAACTACACTTACAACATCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCa                                                      <  8:335441/147‑1 (MQ=255)
                                                                                            aTCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAgg                                                    <  8:337839/148‑1 (MQ=255)
                                                                                               gTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttt                                                                                                                                 <  3:372086/68‑1 (MQ=255)
                                                                                               gTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGcac                                                 <  3:287357/148‑1 (MQ=255)
                                                                                               gTGAACTAAACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCtttt                                                                                                                                 >  4:372086/1‑68 (MQ=255)
                                                                                                tGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACAt                                               >  1:29051/1‑149 (MQ=255)
                                                                                                  aaCTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACAtt                                              <  2:29051/148‑1 (MQ=255)
                                                                                                    cTACACTTACCACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGaa                                           <  2:333873/149‑1 (MQ=255)
                                                                                                       cacTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGACTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATTTATTTTGGTTATAAATAAATTGTTATTTATTATTTTTTTTTTTAATAGATTGAAGGAACATTGaaac                                         >  7:221339/1‑146 (MQ=255)
                                                                                                       caattACAACTTACTGTAAGCCGTATCTGAAAAAGAAGCATCGCCTGATGAGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGTTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAAcc                                         <  8:221339/145‑1 (MQ=255)
                                                                                                          ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                     >  6:359009/1‑149 (MQ=255)
                                                                                                          ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                     >  7:226024/1‑149 (MQ=255)
                                                                                                          ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                     <  1:137673/149‑1 (MQ=255)
                                                                                                          ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                     >  7:211476/1‑149 (MQ=255)
                                                                                                          ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                     >  8:3212/1‑149 (MQ=255)
                                                                                                          ttACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATa                                     >  5:187967/1‑149 (MQ=255)
                                                                                                            aCAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACg                                   >  5:145312/1‑149 (MQ=255)
                                                                                                              aaCTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACgg                                  >  7:105520/1‑148 (MQ=255)
                                                                                                               aCTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTa                                <  4:61358/149‑1 (MQ=255)
                                                                                                                 ttCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTa                                >  1:312137/1‑147 (MQ=255)
                                                                                                                       tAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcc                         >  2:295484/1‑148 (MQ=255)
                                                                                                                             gTATCTGAAAAAGAAGCCTCGACTGATGCGATGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagccca                  <  5:351521/149‑1 (MQ=255)
                                                                                                                             gTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagccca                  <  8:105520/149‑1 (MQ=255)
                                                                                                                             gTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagccca                  <  2:318894/149‑1 (MQ=255)
                                                                                                                              tAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagcccag                 <  2:42442/149‑1 (MQ=255)
                                                                                                                                gCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGa               <  6:187967/149‑1 (MQ=255)
                                                                                                                                gCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAAAAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGa               <  1:125205/149‑1 (MQ=255)
                                                                                                                                aggcacaaaaaaGCCTCGCATCAGGCGAGGCTTCTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagcccag                 >  4:98545/7‑147 (MQ=38)
                                                                                                                                  tGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGACTGCAGCTTTATATGTATTTTGGTTATAAATTAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAGCCCAGCCGAGGAGc             >  6:329742/1‑149 (MQ=255)
                                                                                                                                   gAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGcc            >  7:223977/1‑149 (MQ=255)
                                                                                                                                      aaaaaaGCCTCGCATCAGGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCACAGGAGCCg           >  5:352010/1‑147 (MQ=255)
                                                                                                                                      aaaGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagccca                  <  2:272171/140‑1 (MQ=255)
                                                                                                                                      aaaGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagccca                  >  1:272171/1‑140 (MQ=255)
                                                                                                                                       aagaagCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAg        <  1:223182/149‑1 (MQ=255)
                                                                                                                                        agaagCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagcc                    <  2:92682/136‑1 (MQ=255)
                                                                                                                                        agaagCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAAcccagcc                    >  1:92682/1‑136 (MQ=255)
                                                                                                                                        agaagCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGc       >  7:281398/1‑149 (MQ=255)
                                                                                                                                        agaagCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGACGCAGc       >  8:7711/1‑149 (MQ=255)
                                                                                                                                        aaaagCCTCGCATCAGGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGc       >  4:188049/3‑149 (MQ=255)
                                                                                                                                           agCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGccc     <  7:353080/148‑1 (MQ=255)
                                                                                                                                            gCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCAt   <  6:186632/149‑1 (MQ=255)
                                                                                                                                            gCCTCGCATCAGGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGa               <  2:256333/137‑1 (MQ=255)
                                                                                                                                            gCCTCGCATCAGGCGAGGCTTTTTTGTGCCTGCAAATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGa               >  1:256333/1‑137 (MQ=255)
                                                                                                                                              cTCGCCTGATGCGATGCTTTTTTGTGCATGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATg  <  1:149692/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
AAGGCACTCTCAAAGGCCTGGGCTTCAGCCTGCGTCAAGGTGTCTACCGTGGTAGCGGCGAAAGCGCCCTGGACAAAGACCAGACTCGCTTCATCGTGAACTACACTTACAACTTCCTGTAAGCCGTAGCTGAAAAAGAAGCCTCGCCTGATGCGAGGCTTTTTTGTGCCTGCAGATTTATATGTATTTTGGTTATAAATAAATTGTTATTTATTCTTTTTGTTTTTAATCGATTGCAGGCACATTGAACCCATACGGTACAGAACCCAGCCCAGGAGCCGCAGCCCATG  >  NC_002947/289795‑290084

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: