Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,720,934 T→C A85A (GCT→GCC)  insJ → IS150 transposase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,720,9340TC98.2% 406.7 / ‑1.9 112A85A (GCT→GCC) insJIS150 transposase A
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (61/49);  total (62/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.97e-02

Only 100 of 115 total aligned reads displayed.

AACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGATCCAGAAAAAGCGGCATCA  >  NC_000913/3720787‑3721081
                                                                                                                                                   |                                                                                                                                                   
aaCCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgccgct                                                                                                                                                  <  3:930386/151‑1 (MQ=255)
  cccAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgccgctgc                                                                                                                                                <  2:564340/151‑1 (MQ=255)
   ccAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGgccgctgct                                                                                                                                               <  1:359507/151‑1 (MQ=255)
     aGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTc                                                                                                                                              >  3:599537/1‑150 (MQ=255)
          aGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACttt                                                                                                                                         >  2:418256/1‑150 (MQ=255)
          aGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACttt                                                                                                                                         >  1:71014/1‑150 (MQ=255)
           ggACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTAt                                                                                                                                       <  2:884966/151‑1 (MQ=255)
            gACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTAtg                                                                                                                                      >  2:1100130/1‑151 (MQ=255)
            gACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTAtg                                                                                                                                      >  3:386697/1‑151 (MQ=255)
                 gggTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTAtgcttgc                                                                                                                                 >  2:489919/1‑151 (MQ=255)
                 gggTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTAtgcttgc                                                                                                                                 >  3:285611/1‑151 (MQ=255)
                 gggTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTAtgcttg                                                                                                                                  >  3:1143055/1‑150 (MQ=255)
                       cccTCTATGAAAAACATGGAGAAATAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTTAAAGCTTTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAg                                                                                                                           <  2:185020/151‑1 (MQ=255)
                       cccTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGt                                                                                                                             >  3:760196/1‑149 (MQ=255)
                           ctATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTg                                                                                                                         <  3:760405/149‑1 (MQ=255)
                             aTGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGtt                                                                                                                      <  1:181585/150‑1 (MQ=255)
                               gAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTg                                                                                                                   >  2:245384/1‑151 (MQ=255)
                                aaaaaCATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGt                                                                                                                  <  2:554741/151‑1 (MQ=255)
                                  aaaCATGGAGAAAAAGGTTCAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGt                                                                                                                  >  2:974497/1‑149 (MQ=255)
                                   aaCATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGc                                                                                                               >  3:319311/1‑151 (MQ=255)
                                   aaCATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGc                                                                                                               <  3:688560/151‑1 (MQ=255)
                                     cATGGAGAAAACGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATACAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGcc                                                                                                              >  3:1036768/1‑150 (MQ=255)
                                     cATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCa                                                                                                             <  3:588582/151‑1 (MQ=255)
                                      aTGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGcc                                                                                                              >  2:775256/1‑149 (MQ=255)
                                        ggAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGt                                                                                                          <  3:991525/151‑1 (MQ=255)
                                         gagaAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTg                                                                                                         >  2:166361/1‑151 (MQ=255)
                                            aaaaaGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCt                                                                                                      >  3:247845/1‑151 (MQ=255)
                                              aaaGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTg                                                                                                     >  2:1245777/1‑150 (MQ=255)
                                                aGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGa                                                                                                    >  1:421368/1‑149 (MQ=255)
                                                    ttAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTc                                                                                               <  2:556615/150‑1 (MQ=255)
                                                    ttAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCt                                                                                              >  2:249034/1‑151 (MQ=255)
                                                     tAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCt                                                                                              <  3:1100179/150‑1 (MQ=255)
                                                      aaTTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCt                                                                                              <  4:334291/149‑1 (MQ=255)
                                                       aTTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATg                                                                                           >  3:1111851/1‑151 (MQ=255)
                                                       aTTCCCAAACCTAAAGGCGTTAGCGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATg                                                                                           >  3:829173/1‑151 (MQ=255)
                                                         tCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGa                                                                                          <  4:329280/150‑1 (MQ=255)
                                                         tCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGa                                                                                          <  4:6067/150‑1 (MQ=255)
                                                          cccAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGa                                                                                          >  3:370483/1‑149 (MQ=255)
                                                          cccAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAag                                                                                        <  2:951154/151‑1 (MQ=255)
                                                            cAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAag                                                                                        >  3:879454/1‑149 (MQ=255)
                                                             aaaCCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAgc                                                                                     >  3:725723/1‑151 (MQ=255)
                                                             aaaCCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAgc                                                                                     >  2:416564/1‑151 (MQ=255)
                                                              aaCCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAgcg                                                                                    <  3:64920/151‑1 (MQ=255)
                                                               aCCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAgcgc                                                                                   >  3:430265/1‑151 (MQ=255)
                                                                 cTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCgg                                                                                 >  3:127564/1‑151 (MQ=255)
                                                                   aaaGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCgg                                                                                 >  2:212098/1‑149 (MQ=255)
                                                                   aaaGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCgg                                                                                 <  3:418287/149‑1 (MQ=255)
                                                                    aaGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCgg                                                                                 >  3:218905/1‑148 (MQ=255)
                                                                     aGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGgaga                                                                              >  2:471485/1‑150 (MQ=255)
                                                                     aGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGgaga                                                                              >  4:141340/1‑150 (MQ=255)
                                                                            aGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGtt                                                                      >  3:336310/1‑151 (MQ=255)
                                                                                  gATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTAcgc                                                                 <  2:1371094/150‑1 (MQ=255)
                                                                                   aTCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTAcgcgc                                                               >  4:396272/1‑151 (MQ=255)
                                                                                   aTCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTAcgcgc                                                               >  3:811180/1‑151 (MQ=255)
                                                                                   aTCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTAcgcg                                                                >  4:16710/1‑150 (MQ=255)
                                                                                      cAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCt                                                            <  2:920071/151‑1 (MQ=255)
                                                                                      cAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCt                                                            >  3:842640/1‑151 (MQ=255)
                                                                                       agagTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTc                                                           <  1:120851/151‑1 (MQ=255)
                                                                                       agagTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTc                                                           >  3:1381302/1‑151 (MQ=255)
                                                                                             gCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGAtt                                                     <  3:600874/151‑1 (MQ=255)
                                                                                              cGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTg                                                    >  2:1377349/1‑151 (MQ=255)
                                                                                                 aTTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGc                                                  <  4:100557/150‑1 (MQ=255)
                                                                                                 aTTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCa                                                 <  2:1143000/151‑1 (MQ=255)
                                                                                                 aTTAAGGTCGTGAAAGCTGGGATCGACCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGAAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTg                                                    >  4:466582/1‑148 (MQ=255)
                                                                                                    aaGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCa                                              <  2:285585/151‑1 (MQ=255)
                                                                                                       gTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCaaaa                                           <  3:475678/151‑1 (MQ=255)
                                                                                                         cGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCaaa                                            >  3:1280624/1‑148 (MQ=255)
                                                                                                          gTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGaa                                        <  3:245409/151‑1 (MQ=255)
                                                                                                          gTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGaa                                        <  4:471818/151‑1 (MQ=255)
                                                                                                          gTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGaa                                        <  2:1147974/151‑1 (MQ=255)
                                                                                                            gAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAAc                                      >  1:274434/1‑151 (MQ=255)
                                                                                                            gAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAAc                                      >  2:855748/1‑151 (MQ=255)
                                                                                                            gAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAAc                                      >  3:656813/1‑151 (MQ=255)
                                                                                                             aaaGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAAc                                      >  2:1210420/1‑150 (MQ=255)
                                                                                                              aaGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACAt                                    <  2:829138/151‑1 (MQ=255)
                                                                                                              aaGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACAt                                    >  3:291958/1‑151 (MQ=255)
                                                                                                                gCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTg                                  >  1:424238/1‑151 (MQ=255)
                                                                                                                gCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTg                                  <  3:1286798/151‑1 (MQ=255)
                                                                                                                  tgtgATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGc                                 >  2:460985/1‑150 (MQ=255)
                                                                                                                   gtgATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCa                                >  3:1057218/1‑150 (MQ=255)
                                                                                                                   gtgATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCa                                <  2:655156/150‑1 (MQ=255)
                                                                                                                     gATCGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAat                              <  3:544429/150‑1 (MQ=255)
                                                                                                                      aacGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAatat                            <  1:72480/149‑1 (MQ=255)
                                                                                                                        cGAGCAGCACATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCa                          <  2:1371417/151‑1 (MQ=255)
                                                                                                                             agcaCATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACTAAAAGAAACATTGCAATATCAGTTg                      >  2:937868/1‑150 (MQ=255)
                                                                                                                             agcaCATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGa                     <  2:1111798/151‑1 (MQ=255)
                                                                                                                             agcaCATGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGa                     <  2:386661/151‑1 (MQ=255)
                                                                                                                               cacaTGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGATc                   >  2:234078/1‑151 (MQ=255)
                                                                                                                               cacaTGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGATc                   >  2:982534/1‑151 (MQ=255)
                                                                                                                                acaTGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGATcc                  <  3:1249159/151‑1 (MQ=255)
                                                                                                                                  aTGTCCCTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGATcc                  >  3:358961/1‑149 (MQ=255)
                                                                                                                                    gTCCCTTAATCAGGCTGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGATCCAGaa              <  2:599511/151‑1 (MQ=255)
                                                                                                                                      cccTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGATCCAGaaa             <  3:968283/150‑1 (MQ=255)
                                                                                                                                       ccTTAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGATCCAGaaaaa           <  2:430234/151‑1 (MQ=255)
                                                                                                                                          tAATCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGATCCAGAAAAAGc         <  3:775292/150‑1 (MQ=255)
                                                                                                                                           aaTCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGATCCAGAAAAAGCgg       >  2:1288353/1‑151 (MQ=255)
                                                                                                                                           aaTCAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGATCCAGAAAAAGCgg       >  2:562173/1‑151 (MQ=255)
                                                                                                                                              cAGGCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGATCCAGAAAAAGCGGCAt    >  1:378967/1‑151 (MQ=255)
                                                                                                                                                ggCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGATCCAGAAAAAGCGGCATCa  <  3:117361/151‑1 (MQ=255)
                                                                                                                                                ggCCGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGATCCAGAAAAAGCGGCATCa  <  2:127548/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AACCCAGGTCAGGACATGGGTTGCCCTCTATGAAAAACATGGAGAAAAAGGTTTAATTCCCAAACCTAAAGGCGTTAGTGCTGATCCAGAGTTGCGTATTAAGGTCGTGAAAGCTGTGATCGAGCAGCACATGTCCCTTAATCAGGCTGCTGCTCACTTTATGCTTGCTGGTAGTGGTTCTGTAGCCAGGTGGCTGAAGGTCTATGAAGAGCGCGGAGAAGCTGGTTTACGCGCGCTCAAGATTGGCACCAAAAGAAACATTGCAATATCAGTTGATCCAGAAAAAGCGGCATCA  >  NC_000913/3720787‑3721081

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑