Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 866,388 A→T 6.0% intergenic (‑24/‑180) moeA ← / → iaaA molybdopterin molybdenumtransferase; molybdopterin biosynthesis protein/Isoaspartyl peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913866,3880AT 6.0% 246.4 / 2.9 91intergenic (‑24/‑180)moeA/iaaAmolybdopterin molybdenumtransferase; molybdopterin biosynthesis protein/Isoaspartyl peptidase
Reads supporting (aligned to +/- strand):  ref base A (38/47);  new base T (4/2);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 4.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAA  >  NC_000913/866244‑866536
                                                                                                                                                |                                                                                                                                                    
ccAGAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGa                                                                                                                                                  <  1:221162/149‑1 (MQ=255)
 cAGAATACGAGCATATCACTCTACCAGTGCCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCTATCCGGTGCTAAATTCCATGAAAACACTCCTGCGGAGGCAAAAACGaa                                                                                                                                                 <  2:106176/149‑1 (MQ=255)
 cAGAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTAAGTGGGTTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCTGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGaa                                                                                                                                                 <  2:101092/149‑1 (MQ=255)
  aGAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGGGGGGTGACGCGTGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAAt                                                                                                                                                >  1:1054658/1‑149 (MQ=255)
   gAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAAtt                                                                                                                                               <  1:418977/149‑1 (MQ=255)
    aaTACGACCAAAACTCTGTACCATTGGCAGCGTTTCGTGGGCGGTCAGTGGTGTGAGGCGAGAAAGATTCTCATTCAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAATTTCCATGAAAACACTCCTGCGGAGGAAAAATCGAAttt                                                                                                                                              <  2:793651/149‑1 (MQ=255)
    aaTACGACCAAAACACTGTACCAGGGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAAttt                                                                                                                                              <  2:1507216/149‑1 (MQ=255)
     aTACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGGCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTg                                                                                                                                             <  1:745126/149‑1 (MQ=255)
      aaCGACCAAAACACAGAACCAGTGGCAGCCTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACAACAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGc                                                                                                                                            <  1:918216/148‑1 (MQ=255)
         gACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCTGGCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGATTTTGCCtc                                                                                                                                         <  2:356325/149‑2 (MQ=255)
         gACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCta                                                                                                                                         <  1:107801/149‑1 (MQ=255)
         gACCAAAACACTGTACCAGTGGCAGCGTTTCCAGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGAATTTGCCta                                                                                                                                         >  1:1210835/1‑149 (MQ=255)
         caCCAACACACGGTACCAGTGGCAGCGGTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAATCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGAATTTGCCtc                                                                                                                                         <  1:696997/148‑2 (MQ=255)
          aCCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCtat                                                                                                                                        >  2:699915/1‑149 (MQ=255)
           ccAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCtatt                                                                                                                                       <  1:53656/149‑1 (MQ=255)
           ccAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCtatt                                                                                                                                       <  1:197448/149‑1 (MQ=255)
            cAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGCGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATATGCCtatta                                                                                                                                      >  1:1392602/1‑149 (MQ=255)
                acacTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCTAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTc                                                                                                                                  <  1:864100/149‑1 (MQ=255)
                acacTGTACCAGTGGCAGCGGTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTc                                                                                                                                  <  2:639538/149‑1 (MQ=255)
                       aCCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGATAAAGCAGCGCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACTCTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGaaaaa                                                                                                                           <  2:76374/149‑1 (MQ=255)
                        ccAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAACCATCTCATAAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAAc                                                                                                                          <  2:311734/149‑1 (MQ=255)
                         cAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACg                                                                                                                         >  2:609697/1‑149 (MQ=255)
                          aGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGc                                                                                                                        <  1:538535/149‑1 (MQ=255)
                          aGTGGCAGCGTTTCCTGGGCGGGCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCCGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGc                                                                                                                        <  1:1518668/149‑1 (MQ=255)
                            tGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCca                                                                                                                      <  1:607861/149‑1 (MQ=255)
                            tGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCca                                                                                                                      >  1:1202396/1‑149 (MQ=255)
                             ggCAGCTTTTCCTGGGCGTTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCcac                                                                                                                     <  2:891184/149‑1 (MQ=255)
                              gCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGCGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGccccc                                                                                                                    >  2:1154379/1‑146 (MQ=255)
                              gCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCcaca                                                                                                                    >  1:202068/1‑149 (MQ=255)
                                aGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTTAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGa                                                                                                                  <  2:1145438/149‑1 (MQ=255)
                                  cGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACt                                                                                                                <  2:1060465/149‑1 (MQ=255)
                                      tCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTAt                                                                                                            <  1:486199/149‑1 (MQ=255)
                                      tCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTAt                                                                                                            >  1:414543/1‑149 (MQ=255)
                                        cTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGc                                                                                                          >  2:977457/1‑149 (MQ=255)
                                         tGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGcc                                                                                                         >  2:997050/1‑149 (MQ=255)
                                         tGGGCGGGCAGTGGGGGGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTTAATTCCATGAAAACACTCCTGCGGAGTCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGcc                                                                                                         <  2:116959/149‑1 (MQ=255)
                                          gggCGGTCAGTTGGTTGACGCGTGAAACCATCTCATTAAGCGCGGTGTCGAGAGACATCAATCAGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGcga                                                                                                        <  2:633967/149‑3 (MQ=255)
                                           ggCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCAc                                                                                                       >  2:266273/1‑149 (MQ=255)
                                                tCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCgg                                                                                                  <  1:1113251/149‑1 (MQ=255)
                                                   gTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCg                                                                                               >  1:882810/1‑149 (MQ=255)
                                                      gggTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCTATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCACAATCGAATTTGCCTTTTATGTCAGCAAACCGCCACAGACTGTATTCCACCTCGGGCTTAg                                                                                            >  1:1059752/1‑149 (MQ=255)
                                                             gcgAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCGTATTATGTAAGAAAAACGCCACAGACTGTATGCCACCTCCGGCGTAgcagtggt                                                                                     >  2:1564717/1‑143 (MQ=255)
                                                             gcgAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGGAAATTCCATGAAAACCCTCCTGCGGAGCCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCCCAGACTGTATGCCCCCTCGGGCGGTGCGCTggg                                                                                     >  2:155856/1‑149 (MQ=255)
                                                             gcgAGAAAGCAACTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCAAGAAAACACTCCTGCAGAGGCAAAATCGAATTTGCCTATTATGACAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGAAGCGCTgcg                                                                                     >  1:413530/1‑147 (MQ=255)
                                                                 gaAAGCATCTCATTAAGCGTGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCt                                                                                 <  2:1324578/149‑1 (MQ=255)
                                                                   aaGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGc                                                                               <  2:583844/149‑1 (MQ=255)
                                                                   aaGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGc                                                                               <  2:369118/149‑1 (MQ=255)
                                                                    aGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGcc                                                                              >  1:1283276/1‑149 (MQ=255)
                                                                      cATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCtt                                                                            <  1:1464661/149‑1 (MQ=255)
                                                                      cATCGCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAACATCGAATTTGCCTATTATGTCAGAAAACCGCCACAGACTGTATGCAACCTCGGGCGTAGCGCTGGCTCCTGCCtt                                                                            >  1:1411430/1‑149 (MQ=255)
                                                                         ctcATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGAATGCCACCTCGGGCGTTGCGCTGGGTCCTGCCTTTAc                                                                         >  2:1321540/1‑149 (MQ=255)
                                                                         ctcATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTAc                                                                         >  1:639888/1‑149 (MQ=255)
                                                                         cgcATTAAGCGCGGTGTCGAGCGACATCAATCCGGAGGTTAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAACCGCCGCAGACTGTATGCCAGCTCGGGCGTAGCGCTGGGTCCTGCCTTTaa                                                                         >  2:1346536/3‑148 (MQ=255)
                                                                            aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATg                                                                      <  2:50123/149‑1 (MQ=255)
                                                                            aTTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATg                                                                      >  1:50123/1‑149 (MQ=255)
                                                                             ttAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGc                                                                     <  1:44829/149‑1 (MQ=255)
                                                                              tAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGcc                                                                    <  2:1202396/149‑1 (MQ=255)
                                                                                 gcgcGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCatt                                                                 >  1:370130/1‑148 (MQ=255)
                                                                                   gcgGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGCGGCAAAATCGAAATTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCAGGGTCCTGCCTTTACATGACAtctc                                                               >  1:1200200/1‑149 (MQ=255)
                                                                                   gcgGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATc                                                               <  2:389602/149‑1 (MQ=255)
                                                                                       tgtCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATc                                                           >  1:787863/1‑149 (MQ=255)
                                                                                           gAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAATTCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTc                                                       >  2:80414/1‑149 (MQ=255)
                                                                                            aGCGACATCAATCCGGTGGTAAATTCCATGATAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCCTCTTTCt                                                      <  2:347231/149‑1 (MQ=255)
                                                                                            aGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCt                                                      <  2:581313/149‑1 (MQ=255)
                                                                                             gCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCta                                                     >  1:798321/1‑149 (MQ=255)
                                                                                              cGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCtat                                                    <  2:1054658/149‑1 (MQ=255)
                                                                                                aCATCAATCCGGTGGTAAATTCCATGAAACCACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAACAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGTTCCTGCCTTTACATGCCATATCCATCTTTCTatat                                                  >  2:1376903/1‑149 (MQ=255)
                                                                                                aCATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTatat                                                  >  2:790783/1‑149 (MQ=255)
                                                                                                aCATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTACTGCCTTTACATGCCATATCCATCTTTCTatat                                                  >  1:588673/1‑149 (MQ=255)
                                                                                                 cATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATAtt                                                 <  1:604224/149‑1 (MQ=255)
                                                                                                 aaTCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATAtt                                                 <  2:1445484/148‑1 (MQ=255)
                                                                                                     aaTCCGGTGGTAAATTCCATGAAAACACTCCGGCGGAGGCAAAATCGAATTTGCCTATTTTGTCAGATACACGCCACAGACTGTATGCCCCCTCGGGCGTACCGCTGGGTCCTGCCGTTACATGCCATAGCCATCTTTCGAAATTCaaa                                             >  2:896476/1‑149 (MQ=255)
                                                                                                          ggtggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTg                                        >  2:1036800/1‑149 (MQ=255)
                                                                                                          ggtggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTg                                        >  2:492255/1‑149 (MQ=255)
                                                                                                            gggtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGaa                                      <  2:1481814/148‑1 (MQ=255)
                                                                                                              gtAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATg                                    <  2:1356926/149‑1 (MQ=255)
                                                                                                                    tCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAAt                              <  2:414543/149‑1 (MQ=255)
                                                                                                                     ccATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCCATATTCAAAAATTGAATGAGTAAtt                             >  1:247516/1‑149 (MQ=255)
                                                                                                                      cATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTc                            <  2:410291/149‑1 (MQ=255)
                                                                                                                      cATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTc                            >  1:7638/1‑149 (MQ=255)
                                                                                                                       aTGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCa                           >  1:612108/1‑149 (MQ=255)
                                                                                                                            aaCACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATaaaa                      >  1:558097/1‑149 (MQ=255)
                                                                                                                                cTCCTTCGGAGTCAAAATCGAATTTGCATTTTATGTCAGAAAAACGCCACAGACTGTATGACACTTCGGGCGTATCGCTGTGTCCTTCCTTTACATGCCATATCCATCTTTCTATATTAAAAAATTGAATGAGTAATTCATAAAAATTc                  <  2:296346/149‑1 (MQ=255)
                                                                                                                                cTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTc                  <  1:699915/149‑1 (MQ=255)
                                                                                                                                 tCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGATTGAGTAATTCATAAAAATTCt                 >  1:1108606/1‑149 (MQ=255)
                                                                                                                                 tCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCt                 <  1:874984/149‑1 (MQ=255)
                                                                                                                                  ccTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTg                <  1:609697/149‑1 (MQ=255)
                                                                                                                                  acataataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTg                >  2:1265481/8‑149 (MQ=255)
                                                                                                                                    tGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACTTGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGat              >  2:350337/1‑149 (MQ=255)
                                                                                                                                       ataGGCAAATTCGATTTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATTTCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATAtt           >  2:74464/3‑149 (MQ=255)
                                                                                                                                         aGGCAAAATCGACTTTGCCTATTATGTCAGAAAAACGCCACAGTCTGTATGCCACCTCGGGCGTAGCGCTGGGGCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGACTGAGTAATTCATAAAAATTCTGATATTta         <  2:1342934/149‑1 (MQ=255)
                                                                                                                                          ttaaaaaTAGAATTTGCCTATTATGTCAGAAAAAAGACAATGACTTTATTCCAACTATTTCGTATCGCTTGGTACTGACTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTtat        <  1:97651/146‑1 (MQ=255)
                                                                                                                                                aTCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCCCCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCaaa  >  2:1496581/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCAGAATACGACCAAAACACTGTACCAGTGGCAGCGTTTCCTGGGCGGTCAGTGGGGTGACGCGAGAAAGCATCTCATTAAGCGCGGTGTCGAGCGACATCAATCCGGTGGTAAATTCCATGAAAACACTCCTGCGGAGGCAAAATCGAATTTGCCTATTATGTCAGAAAAACGCCACAGACTGTATGCCACCTCGGGCGTAGCGCTGGGTCCTGCCTTTACATGCCATATCCATCTTTCTATATTCAAAAATTGAATGAGTAATTCATAAAAATTCTGATATTTATAGCAAA  >  NC_000913/866244‑866536

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