Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,768,022 G→A 6.7% R222C (CGT→TGT) ‡ ydiJ ← putative FAD‑linked oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,768,0220GA 6.7% 275.6 / 6.5 90R222C (CGT→TGT) ‡ydiJputative FAD‑linked oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (40/43);  new base A (2/4);  total (43/47)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

AATAAAGGCCAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGTT  >  NC_000913/1767873‑1768170
                                                                                                                                                     |                                                                                                                                                    
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                                                                                                                       <  1:939270/150‑1 (MQ=255)
                               aGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGc                                                                                                                       >  2:651436/1‑150 (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                                                                              <  2:400787/150‑1 (MQ=255)
                                                                        gACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTc                                                                              <  1:597660/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:1255171/150‑1 (MQ=255)
                                                                                       gTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTg                                                               >  1:972600/1‑150 (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:81335/1‑150 (MQ=255)
                                                                                                  gAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGaa                                                    >  1:1225998/1‑150 (MQ=255)
                                                                                                  gAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGaa                                                    >  2:605978/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             >  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)
                                                                                                                                                  gACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACg    >  1:815591/1‑150 (MQ=255)
                                                                                                                                                  gACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCAAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACg    >  2:469609/1‑150 (MQ=255)
                                                                                                                                                    cGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGtt  <  1:31382/150‑1 (MQ=255)
                                                                                                                                                    cGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGtt  >  1:694838/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AATAAAGGCCAGCGTCCCTTCTGAACCCGTCAGAATGCGCGTCAGGTCGAACTCGGTCATCTCATCGTTAAAGACATGACGCAGATCGTAACCGGTAAGAAAGCGGTTAAGTTTGGGGAAGTTGTCGATAATTAACTGGCGTTGCTGACGGCAACGTTGATAAACCGTGTTATAAATTCGCCCGATTGTGGTATTGGATTTACCCAGCGTTTCCGCCAATTCGACGGGTAAAGGTTGCGTATCGAGAATATCGCCCCCCAACAACACCGCGCGTACGCCAAGTACGTGATCTGACGTT  >  NC_000913/1767873‑1768170

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: