Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,544,911 T→A 6.3% V203E (GTG→GAG)  fimC → periplasmic chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,544,9110TA 6.3% 219.1 / 4.4 81V203E (GTG→GAG) fimCperiplasmic chaperone
Reads supporting (aligned to +/- strand):  ref base T (35/40);  new base A (3/2);  total (38/43)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.21e-01

GCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAG  >  NC_000913/4544763‑4545060
                                                                                                                                                    |                                                                                                                                                     
gCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTg                                                                                                                                                      <  2:115059/150‑1 (MQ=255)
gCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTg                                                                                                                                                      <  2:572216/150‑1 (MQ=255)
   aaaTTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCt                                                                                                                                                   >  1:1690361/1‑150 (MQ=255)
     aTTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTcc                                                                                                                                                 <  2:810318/150‑1 (MQ=255)
        aGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAAt                                                                                                                                              <  1:1670820/150‑1 (MQ=255)
        aGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAAt                                                                                                                                              <  1:1327951/150‑1 (MQ=255)
           gTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGAGGCACCAATgca                                                                                                                                           >  2:1226410/1‑148 (MQ=255)
           gTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGAGAATTCTCTGACGCTGATTAAACCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGAATTGGTGCCTCCAATggg                                                                                                                                           >  1:1031819/1‑150 (MQ=255)
             tGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCg                                                                                                                                         <  2:872526/150‑1 (MQ=255)
             tGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCg                                                                                                                                         <  1:1707076/150‑1 (MQ=255)
              gCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATAACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGAGGCACCAATgcattt                                                                                                                                        <  1:282146/150‑6 (MQ=255)
              gCCACCCGATCAGGCCGCACAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGGTGAATGACGGAACCCGGGTTCTTGAAAATGCATTGGAGGCACCAATgcattt                                                                                                                                        <  2:1522102/150‑6 (MQ=255)
               ccACCCGATCAGGCCTCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGaa                                                                                                                                       >  2:1427310/1‑150 (MQ=255)
               ccACCCGATCAGGCCGTAGAAAAATTAAGATTTCTTCGTAGCGCGAATTCTATGACGCTGATTAACCCGACACCCTATTACCAGACTGTAACATAGTTGAAAGCCGGAACCATGGTTCTATAAAATGCATTCGTGACTCAAATGGGAGaa                                                                                                                                       >  2:1532639/1‑150 (MQ=255)
               ccACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGaa                                                                                                                                       >  2:1210555/1‑150 (MQ=255)
                cACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGaaa                                                                                                                                      >  1:1613942/1‑150 (MQ=255)
                 aCCCGATCAGGCAGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTACCCCGACACCCTATTACCTGACCGTAACAGAGTTGAATGCCGGACCCCGGGTGCTTGAAAATGCATTGGTGCCTCCTATGGGCGAAAg                                                                                                                                     >  1:1080950/1‑150 (MQ=255)
                      aTCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACgg                                                                                                                                <  2:1321098/150‑1 (MQ=255)
                         aGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTTGTGCCTCCAATGGGCGAAAGCACGGTTa                                                                                                                             >  2:1412047/1‑150 (MQ=255)
                               cAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGc                                                                                                                       <  2:1547182/150‑1 (MQ=255)
                                aGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGcc                                                                                                                      <  2:362760/150‑1 (MQ=255)
                                   aaaaTTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTc                                                                                                                   <  2:1488023/150‑1 (MQ=255)
                                    aaaTTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTCCCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCt                                                                                                                  <  2:1341771/150‑1 (MQ=255)
                                    aaaTTAAGATTTCGGCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCt                                                                                                                  <  1:638843/150‑1 (MQ=255)
                                        tAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATg                                                                                                              <  2:1613942/150‑1 (MQ=255)
                                           taTTTCGGAGTAGCGCGAATTCTCTGAAGCTTATTAACCCGACACCCTTTTAACTGACGGTAACAGAGTTGAATGCCGGAAGCCGGGGTGTTGAAAATGCATTGGGGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAg                                                                                                           <  2:335892/149‑1 (MQ=255)
                                             tttCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATCGCCTTCTGATGCAGGa                                                                                                         <  1:1532299/150‑1 (MQ=255)
                                              ttCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGaa                                                                                                        >  2:365543/1‑150 (MQ=255)
                                               tcgtcgTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGTCGGTAACAGAGTTGAATTCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCTGGAAg                                                                                                       <  1:530655/150‑1 (MQ=255)
                                                cgtcgtAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGc                                                                                                      <  2:1214435/150‑1 (MQ=255)
                                                     tAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATGGCCTTCTGATGCAGGAAGCAatat                                                                                                 >  1:1408274/1‑150 (MQ=255)
                                                      aGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATAtt                                                                                                <  2:1640916/150‑1 (MQ=255)
                                                          cgAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATAttactt                                                                                            <  1:1435758/150‑1 (MQ=255)
                                                          cgAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATAttactt                                                                                            <  2:880297/150‑1 (MQ=255)
                                                          cgAATTCTATGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATAttactt                                                                                            <  1:1111965/150‑1 (MQ=255)
                                                            aaTTCTCTGACGCTGATTACCCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATAttacttac                                                                                          <  1:1664069/150‑1 (MQ=255)
                                                              ttCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCg                                                                                        <  1:7820/150‑1 (MQ=255)
                                                                 tctGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAAc                                                                                     >  2:1478782/1‑150 (MQ=255)
                                                                 tctGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAAc                                                                                     <  2:1132955/150‑1 (MQ=255)
                                                                 tctGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAAc                                                                                     <  2:1127037/150‑1 (MQ=255)
                                                                  ctGACGCTGTTTAACCCGACACCCTATTACGTGACGGTAACAGAGTTGAATGCCGGAACCCGGGGTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACa                                                                                    <  2:1006742/150‑1 (MQ=255)
                                                                      cGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCCGGAAGCAATATTACTTACCGACCAATaa                                                                                >  1:854299/1‑150 (MQ=255)
                                                                      cGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATaa                                                                                >  1:243172/1‑150 (MQ=255)
                                                                      cGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATaa                                                                                >  1:1647093/1‑150 (MQ=255)
                                                                             tAACCCTACACCCTATTCCCTGACGGTAACAGAGTTGAATGCCGTAACCCGGGTGCGTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTa                                                                         <  2:1080950/150‑1 (MQ=255)
                                                                              aaCCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTAt                                                                        <  1:1155533/150‑1 (MQ=255)
                                                                               aCCGGAGACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAAGCCGGGTTCTTGAAAATGCTTTGGTGCCTCCATTTGGCGAAAGCACGGTTAAATTGCATTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATg                                                                       <  1:777782/150‑1 (MQ=255)
                                                                               aCCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATg                                                                       >  2:295104/1‑150 (MQ=255)
                                                                                cccGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATgg                                                                      <  1:476013/150‑1 (MQ=255)
                                                                                       cccTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACtt                                                               >  1:1442749/1‑150 (MQ=255)
                                                                                       cccTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACtt                                                               >  1:1539309/1‑150 (MQ=255)
                                                                                          tATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTAcc                                                            <  1:876029/150‑1 (MQ=255)
                                                                                           aTTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTAccc                                                           <  1:465010/150‑1 (MQ=255)
                                                                                           aTTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTAccc                                                           <  1:365543/150‑1 (MQ=255)
                                                                                            ttACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTAcccc                                                          >  1:450158/1‑150 (MQ=255)
                                                                                              aCCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCAAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCa                                                        >  1:117017/1‑150 (MQ=255)
                                                                                               ccTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCaa                                                       >  1:1658459/1‑150 (MQ=255)
                                                                                                     ggTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAAGGAc                                                 >  1:536190/1‑150 (MQ=255)
                                                                                                       tAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACgg                                               >  1:805004/1‑150 (MQ=255)
                                                                                                       tAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACgg                                               >  1:1428282/1‑150 (MQ=255)
                                                                                                        aaCAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACggg                                              <  1:1055774/150‑1 (MQ=255)
                                                                                                         aCAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACTATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGc                                             <  1:150512/150‑1 (MQ=255)
                                                                                                          cAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCg                                            >  1:302719/1‑150 (MQ=255)
                                                                                                                tGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATgg                                      >  1:1456586/1‑150 (MQ=255)
                                                                                                                tGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATgg                                      >  1:1712013/1‑150 (MQ=255)
                                                                                                                tGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATgg                                      >  1:522622/1‑150 (MQ=255)
                                                                                                                   aTGCCGGAACCCGGGTTCTTGAAAATGCATTGGAGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAAt                                   >  1:897670/1‑150 (MQ=255)
                                                                                                                     gCCGGATCCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGACAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCACCAAAATGACGGGCGTAATGGAATaa                                 >  1:1001947/1‑150 (MQ=255)
                                                                                                                            cccGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAggg                          >  1:1400652/1‑150 (MQ=255)
                                                                                                                            cccGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAggg                          >  2:993417/1‑150 (MQ=255)
                                                                                                                             ccGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAgggg                         <  2:1664759/150‑1 (MQ=255)
                                                                                                                             ccGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAgggg                         >  2:1021389/1‑150 (MQ=255)
                                                                                                                                      tGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTc                >  1:1651692/1‑150 (MQ=255)
                                                                                                                                         aactGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGCGCGTAATGGAATAACGCAGGGGGAATTTTTCGcc             <  1:815085/147‑1 (MQ=255)
                                                                                                                                         aaaTGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGcc             >  1:1406636/1‑150 (MQ=255)
                                                                                                                                          cgcccATTGGAGGCACCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCt            >  1:875599/5‑150 (MQ=255)
                                                                                                                                           aTGCATTGGTGCATCCATTGGGCGAAAGCAAGGTTTAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATTATTATGGCGAACTTACCCCCAAAATGACGGGCGTAATGTAATAACGCAGGGGGAATTTTTCGCCTg           <  1:393670/150‑1 (MQ=255)
                                                                                                                                             gCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGaa         >  2:827147/1‑150 (MQ=255)
                                                                                                                                              cATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAAt        <  1:1210555/150‑1 (MQ=255)
                                                                                                                                               aTTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATa       <  1:1427310/150‑1 (MQ=255)
                                                                                                                                                  ggTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATaaaa    >  2:1579874/1‑150 (MQ=255)
                                                                                                                                                    tGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAg  >  2:1356565/1‑150 (MQ=255)
                                                                                                                                                    tGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAg  >  2:704848/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCTAAATTAGCGTTGCCACCCGATCAGGCCGCAGAAAAATTAAGATTTCGTCGTAGCGCGAATTCTCTGACGCTGATTAACCCGACACCCTATTACCTGACGGTAACAGAGTTGAATGCCGGAACCCGGGTTCTTGAAAATGCATTGGTGCCTCCAATGGGCGAAAGCACGGTTAAATTGCCTTCTGATGCAGGAAGCAATATTACTTACCGAACAATAAATGATTATGGCGCACTTACCCCCAAAATGACGGGCGTAATGGAATAACGCAGGGGGAATTTTTCGCCTGAATAAAAAG  >  NC_000913/4544763‑4545060

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