Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,768,020 A→T 6.4% R222R (CGT→CGA) ‡ ydiJ ← putative FAD‑linked oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,768,0200AT 6.4% 276.1 / 4.5 91R222R (CGT→CGA) ‡ydiJputative FAD‑linked oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (40/45);  new base T (2/4);  total (42/49)
Fisher's exact test for biased strand distribution p-value = 6.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

TAATAAAGGCCAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACG  >  NC_000913/1767872‑1768168
                                                                                                                                                    |                                                                                                                                                    
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                                                                                                                          <  2:315528/144‑1 (MQ=255)
                                 gAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATaaa                                                                                                                          >  1:315528/1‑144 (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                                                                                 <  1:168806/150‑1 (MQ=255)
                                                                    ttAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGcc                                                                                 <  2:749864/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:1163864/150‑1 (MQ=255)
                                                                             tGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACg                                                                        <  2:150975/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                                                  >  2:605978/1‑150 (MQ=255)
                                                                                                   gAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGaa                                                  >  1:1225998/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:154293/150‑1 (MQ=255)
                                                                                                                                    ttAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCaa                 >  1:825131/1‑150 (MQ=255)
                                                                                                                                     tAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAg                <  1:324859/150‑1 (MQ=255)
                                                                                                                                      aaCTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGt               <  1:651436/150‑1 (MQ=255)
                                                                                                                                          ggCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGt           <  2:1116127/150‑1 (MQ=255)
                                                                                                                                          ggCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGt           >  1:494819/1‑150 (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)
                                                                                                                                                   gACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACg  >  1:815591/1‑150 (MQ=255)
                                                                                                                                                   gACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCAAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACg  >  2:469609/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAATAAAGGCCAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACG  >  NC_000913/1767872‑1768168

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: