Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,768,018 T→C 6.7% Q223R (CAG→CGG)  ydiJ ← putative FAD‑linked oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,768,0180TC 6.7% 267.8 / 7.1 89Q223R (CAG→CGG) ydiJputative FAD‑linked oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (38/45);  new base C (2/4);  total (40/49)
Fisher's exact test for biased strand distribution p-value = 6.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.77e-01

CGGTAATAAAGGCCAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGA  >  NC_000913/1767869‑1768166
                                                                                                                                                     |                                                                                                                                                    
cGGTAATAAAGGCCAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCt                                                                                                                                                      >  1:749864/1‑150 (MQ=255)
   taataaAGGCCAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGAc                                                                                                                                                   <  2:123428/150‑1 (MQ=255)
    aataaAGGCCAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACg                                                                                                                                                  <  2:1168551/150‑1 (MQ=255)
     ataaAGGCCAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACgg                                                                                                                                                 <  2:667414/150‑1 (MQ=255)
           gCCAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGt                                                                                                                                           <  1:296185/150‑1 (MQ=255)
            ccAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGtt                                                                                                                                          >  1:568968/1‑150 (MQ=255)
            ccAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGtt                                                                                                                                          >  1:662466/1‑150 (MQ=255)
            ccAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCCGTCAGCAACgcc                                                                                                                                          <  1:1078997/150‑3 (MQ=255)
             cAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTg                                                                                                                                         <  2:609119/150‑1 (MQ=255)
             cAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCCGTCAGCAACgcca                                                                                                                                         <  1:1116279/150‑4 (MQ=255)
               gCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCCGTCAGCAACgccagt                                                                                                                                       <  1:1423032/150‑6 (MQ=255)
                 gTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATaa                                                                                                                                     >  1:1116127/1‑150 (MQ=255)
                 gTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCCGTCAGCAACgccagtta                                                                                                                                     <  2:454302/150‑8 (MQ=255)
                   cccTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAAc                                                                                                                                   >  2:1255171/1‑150 (MQ=255)
                    ccTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGCCAACGTTGATAAAcc                                                                                                                                  >  1:768025/1‑150 (MQ=255)
                        cTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCgtgt                                                                                                                              >  2:168806/1‑150 (MQ=255)
                          gAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTta                                                                                                                            >  2:597660/1‑150 (MQ=255)
                          gAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTta                                                                                                                            >  2:936467/1‑150 (MQ=255)
                          gAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTta                                                                                                                            <  2:993907/150‑1 (MQ=255)
                            aCCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTtata                                                                                                                          <  1:506182/150‑1 (MQ=255)
                             cccGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATGAACCGTGTTATaa                                                                                                                         >  2:1009214/1‑150 (MQ=255)
                             cccGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATaa                                                                                                                         <  2:790249/150‑1 (MQ=255)
                               cGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAAt                                                                                                                       >  2:154293/1‑150 (MQ=255)
                                  cAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCg                                                                                                                    <  2:1321615/150‑1 (MQ=255)
                                   aGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGc                                                                                                                   >  2:651436/1‑150 (MQ=255)
                                   aGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGc                                                                                                                   <  1:939270/150‑1 (MQ=255)
                                    gAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATaaa                                                                                                                        >  1:315528/1‑144 (MQ=255)
                                    gAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATaaa                                                                                                                        <  2:315528/144‑1 (MQ=255)
                                         cgcgTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGAtt                                                                                                             >  2:157020/1‑150 (MQ=255)
                                             tCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTgg                                                                                                         <  1:100023/150‑1 (MQ=255)
                                             tCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTgg                                                                                                         <  1:660030/150‑1 (MQ=255)
                                                  tCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTg                                                                                                    <  2:1237502/150‑1 (MQ=255)
                                                       cTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGAttt                                                                                               <  1:1209068/150‑1 (MQ=255)
                                                       cTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGAttt                                                                                               >  1:51893/1‑150 (MQ=255)
                                                        tCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTa                                                                                              <  2:568968/150‑1 (MQ=255)
                                                                ctcatcGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACCTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGt                                                                                      >  2:588767/1‑150 (MQ=255)
                                                                     cGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCg                                                                                 <  1:157020/150‑1 (MQ=255)
                                                                       ttAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGcc                                                                               <  2:749864/150‑1 (MQ=255)
                                                                       ttAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGcc                                                                               <  1:168806/150‑1 (MQ=255)
                                                                        tAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCa                                                                              >  1:150975/1‑150 (MQ=255)
                                                                         aaaGACATTACGCAGATCGTAACAGGTAAGAAAGCGGTTAATTTTGGGGAAGTTTTCGATAATTAACTGGCGTTGCTGAATGCAACGTTGATAAACCGTGTTATAAATTCGCCAGATTGTGGTATTGGATTTACCCAGAGTTTCCGCCaa                                                                             >  1:1229265/1‑150 (MQ=255)
                                                                          aaGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAAt                                                                            >  1:164207/1‑150 (MQ=255)
                                                                            gACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTc                                                                          <  1:597660/150‑1 (MQ=255)
                                                                            gACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTc                                                                          <  2:400787/150‑1 (MQ=255)
                                                                              cATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACAGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGa                                                                        >  1:1212922/1‑150 (MQ=255)
                                                                               attaCGCAGATCGTAACCGGTAATAAAGCGTTTAAGTTTGGGGAAGTTGTCTATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACAGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGAc                                                                       <  2:51893/147‑1 (MQ=255)
                                                                               aTGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGAc                                                                       <  1:39417/150‑1 (MQ=255)
                                                                                tGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACg                                                                      <  2:150975/150‑1 (MQ=255)
                                                                                tGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACg                                                                      <  2:1163864/150‑1 (MQ=255)
                                                                                    gCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTa                                                                  >  1:137819/1‑150 (MQ=255)
                                                                                     cAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTaa                                                                 >  2:1320975/1‑150 (MQ=255)
                                                                                        aTCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAgg                                                              >  1:289239/1‑150 (MQ=255)
                                                                                           gTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTg                                                           >  1:972600/1‑150 (MQ=255)
                                                                                           gTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTg                                                           <  1:1255171/150‑1 (MQ=255)
                                                                                           gTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTg                                                           >  1:1303725/1‑150 (MQ=255)
                                                                                             aaCCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCg                                                         <  2:426112/150‑1 (MQ=255)
                                                                                               ccGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTa                                                       >  2:803811/1‑150 (MQ=255)
                                                                                                cGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTAt                                                      <  2:164207/150‑1 (MQ=255)
                                                                                                 ggTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTAtg                                                     <  1:803811/150‑2 (MQ=255)
                                                                                                   tAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCga                                                   <  2:289239/150‑1 (MQ=255)
                                                                                                      gAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGaa                                                >  1:1225998/1‑150 (MQ=255)
                                                                                                      gAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGaa                                                >  2:605978/1‑150 (MQ=255)
                                                                                                      gAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGaa                                                >  1:81335/1‑150 (MQ=255)
                                                                                                       aaaGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAat                                               <  2:267164/150‑1 (MQ=255)
                                                                                                           cGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCg                                           <  2:395175/150‑1 (MQ=255)
                                                                                                            ggTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGc                                          >  1:1386056/1‑150 (MQ=255)
                                                                                                             gTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGcc                                         >  2:324859/1‑150 (MQ=255)
                                                                                                                  gTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTTACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCca                                    >  2:914247/1‑150 (MQ=255)
                                                                                                                    ttGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCcaac                                  <  1:605978/150‑1 (MQ=255)
                                                                                                                     tGGGGAAGTTGTCGATAATTAACTGGCGTTGCTTACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCcaaca                                 <  1:914247/150‑1 (MQ=255)
                                                                                                                     tGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCcaaca                                 <  2:1225998/150‑1 (MQ=255)
                                                                                                                        ggAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAacaaca                              >  1:507723/1‑150 (MQ=255)
                                                                                                                         gAAGTTGTCGATATTTAACTGGCGTTGCTGACGGCAAGGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAacac                             <  2:768025/150‑1 (MQ=255)
                                                                                                                             ttGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACcgcg                         <  1:1320975/150‑1 (MQ=255)
                                                                                                                              tGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACcgcgc                        <  1:588767/150‑1 (MQ=255)
                                                                                                                               gTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACcgcgcg                       >  2:257358/1‑150 (MQ=255)
                                                                                                                                       ttAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCaa               >  1:825131/1‑150 (MQ=255)
                                                                                                                                       ttAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCaa               <  1:154293/150‑1 (MQ=255)
                                                                                                                                        tAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAg              <  1:324859/150‑1 (MQ=255)
                                                                                                                                         aaCTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGt             <  1:651436/150‑1 (MQ=255)
                                                                                                                                             ggCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGt         >  1:494819/1‑150 (MQ=255)
                                                                                                                                             ggCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGt         <  2:1116127/150‑1 (MQ=255)
                                                                                                                                             aacGTTGCCGTCAGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGt         >  2:154705/3‑150 (MQ=255)
                                                                                                                                             aacGTTGCCGTCAGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCACCAACAACACCGCGCGTACGCCAAGTACGt         >  1:912454/3‑150 (MQ=255)
                                                                                                                                                 ttGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATc     <  2:137819/150‑1 (MQ=255)
                                                                                                                                                 ttGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATc     <  2:1107000/150‑1 (MQ=255)
                                                                                                                                                 ttGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATc     <  2:81335/150‑1 (MQ=255)
                                                                                                                                                  tGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCt    >  1:941602/1‑150 (MQ=255)
                                                                                                                                                   gCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTg   >  1:274582/1‑150 (MQ=255)
                                                                                                                                                    cTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGa  >  1:1111642/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGTAATAAAGGCCAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGA  >  NC_000913/1767869‑1768166

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