Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 747,900 G→C 9.9% intergenic (‑131/+21) abrB ← / ← ybgO regulator of aidB expression; inner membrane protein/putative fimbrial protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913747,9000GC 9.9% 310.6 / 20.7 91intergenic (‑131/+21)abrB/ybgOregulator of aidB expression; inner membrane protein/putative fimbrial protein
Reads supporting (aligned to +/- strand):  ref base G (41/41);  new base C (4/5);  total (45/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

CCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTG  >  NC_000913/747752‑748041
                                                                                                                                                     |                                                                                                                                             
ccACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑Ccggcg                                                                                                                                               <  4:227017/149‑1 (MQ=255)
         aaCTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGc                                                                                                                                      >  2:1413829/1‑149 (MQ=255)
            tGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGcatc                                                                                                                                   >  1:81494/1‑149 (MQ=255)
                 tCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑GCCGCCGGAAAGGGcat                                                                                                                                    >  2:230590/1‑143 (MQ=255)
                 tCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑GCCGCCGGAAAGGGcat                                                                                                                                    <  1:230590/143‑1 (MQ=255)
                            cTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑GCCGCCGGAAAGGGca                                                                                                                                     <  4:811127/131‑1 (MQ=255)
                            cTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑GCCGCCGGAAAGGGca                                                                                                                                     >  3:811127/1‑131 (MQ=255)
                             tCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTtata                                                                                                                  >  1:580475/1‑149 (MQ=255)
                             tCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTtata                                                                                                                  <  3:782583/149‑1 (MQ=255)
                               attattTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCaaa                                                                                                                                          >  3:1310744/1‑123 (MQ=255)
                               attattTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTtat                                                                                                                   >  1:1089568/1‑146 (MQ=255)
                               attattTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTtat                                                                                                                   <  2:1089568/146‑1 (MQ=255)
                               attattTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAACAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCaaa                                                                                                                                          <  4:1310744/123‑1 (MQ=255)
                                ttattTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGtt                                                                                                               >  1:782835/1‑149 (MQ=255)
                                  attTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTaa                                                                                                             <  1:1413829/149‑1 (MQ=255)
                                    tttCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                           <  2:580475/149‑1 (MQ=255)
                                        tctTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAg                                                                                                                                         >  4:1161464/1‑115 (MQ=255)
                                        tctTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAg                                                                                                                                         <  3:1161464/115‑1 (MQ=255)
                                        tctTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGt                                                                                                                      >  3:1397991/1‑134 (MQ=255)
                                        tctTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGt                                                                                                                      <  4:1397991/134‑1 (MQ=255)
                                        tctTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaa                                                                                                       >  3:1478586/1‑149 (MQ=255)
                                          tttGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCACAg                                                                                                                                         <  3:1256725/113‑1 (MQ=255)
                                          tttGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAg                                                                                                                                         >  4:1256725/1‑113 (MQ=255)
                                           ttGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTaaaaa                                                                                                    <  2:81494/149‑1 (MQ=255)
                                              gtgtATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                 >  4:1622346/1‑149 (MQ=255)
                                                  aTCGATACAGCATAAAGAAAGCTGGGTCAGGAAAAAATGACAGCAATCAAGGAAAAGGGGGAAATCAACAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGTGAAATCAGTTATAGTTAATCTTAAACAGcaccac                                                                                             >  3:1176214/1‑149 (MQ=255)
                                                       tACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAg                                                                                                                       >  3:204777/1‑118 (MQ=255)
                                                       tACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAg                                                                                                                       <  4:204777/118‑1 (MQ=255)
                                                           gCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGaa                                                                                    <  2:782835/149‑1 (MQ=255)
                                                            cATAAAGAAAGCTGGGCCAGAAAAAATTAGCAGCAAGCAAGGAGAAGGGGAAAATCAGCGATTCTCTGACAGAGATGCCCTTT‑CCGGCGGCAATGGGCATCATGGTAAATCAGATATAGGTAATCTTTAAAAGCACCACTGTATCGAAg                                                                                   <  4:1176214/149‑1 (MQ=255)
                                                             aTAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAgg                                                                                  >  2:749953/1‑149 (MQ=255)
                                                             aTAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAgg                                                                                  >  1:1849969/1‑149 (MQ=255)
                                                              tAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                 <  2:751804/133‑1 (MQ=255)
                                                              tAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                 >  1:751804/1‑133 (MQ=255)
                                                              tAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAggg                                                                                 <  3:1865766/149‑1 (MQ=255)
                                                                   aaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAg                                                                            <  4:1478586/149‑1 (MQ=255)
                                                                        tGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCa                                                                                                                        <  1:304428/100‑1 (MQ=255)
                                                                        tGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCa                                                                                                                        >  2:304428/1‑100 (MQ=255)
                                                                            ccAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAg                                                                                                   >  3:480509/1‑117 (MQ=255)
                                                                            ccAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAg                                                                                                   <  4:480509/117‑1 (MQ=255)
                                                                            ccAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgg                                                                   <  2:1849969/149‑1 (MQ=255)
                                                                                gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcac                                                                                                <  3:324776/116‑1 (MQ=255)
                                                                                gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcac                                                                                                >  4:324776/1‑116 (MQ=255)
                                                                                     aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCa                                                                                                                        >  3:1295548/1‑87 (MQ=255)
                                                                                     aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCa                                                                                                                        <  4:1295548/87‑1 (MQ=255)
                                                                                     aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTtata                                                                                                                  >  1:839314/1‑93 (MQ=255)
                                                                                     aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTtata                                                                                                                  <  2:839314/93‑1 (MQ=255)
                                                                                     aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGcc                                                                              <  4:1701216/129‑1 (MQ=255)
                                                                                     aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGcc                                                                              >  3:1701216/1‑129 (MQ=255)
                                                                                     aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTg                                                          <  1:749953/149‑1 (MQ=255)
                                                                                           aGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGATATGCCCTTTGCCGCCGG‑AAAGGGc                                                                                                                                      >  3:1192769/1‑67 (MQ=25)
                                                                                           aGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGATATGCCCTTT‑GCCGCCGGAAAGGGc                                                                                                                                      <  4:1192769/67‑1 (MQ=25)
                                                                                           aGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCt                                                           >  4:1366067/1‑142 (MQ=255)
                                                                                           aGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCt                                                           <  3:1366067/142‑1 (MQ=255)
                                                                                              aaTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAgg                                                                                                                                        <  2:1419948/62‑1 (MQ=255)
                                                                                              aaTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAgg                                                                                                                                        >  1:1419948/1‑62 (MQ=255)
                                                                                              aaTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                   <  2:282718/147‑1 (MQ=255)
                                                                                              aaTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                   >  1:282718/1‑147 (MQ=255)
                                                                                               aTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgggg                                                                 >  3:970/1‑132 (MQ=255)
                                                                                               aTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgggg                                                                 <  4:970/132‑1 (MQ=255)
                                                                                                tCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaa                                               >  4:419159/1‑149 (MQ=255)
                                                                                                 cAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaaa                                              <  3:1622346/149‑1 (MQ=255)
                                                                                                 cAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaaa                                              <  3:419159/149‑1 (MQ=255)
                                                                                                   aGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                 <  1:1174312/96‑1 (MQ=255)
                                                                                                   aGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                 >  2:1174312/1‑96 (MQ=255)
                                                                                                   aGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                 >  4:872349/1‑96 (MQ=255)
                                                                                                   aGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                 <  3:872349/96‑1 (MQ=255)
                                                                                                     gAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTc                                          >  2:1718664/1‑149 (MQ=255)
                                                                                                        aaGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                           >  1:440628/1‑81 (MQ=255)
                                                                                                        aaGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                           <  2:440628/81‑1 (MQ=255)
                                                                                                           gggAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCt                                                           >  4:1759877/1‑126 (MQ=255)
                                                                                                           gggAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCt                                                           <  3:1759877/126‑1 (MQ=255)
                                                                                                              aaaaTCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                   <  1:950518/131‑1 (MQ=255)
                                                                                                              aaaaTCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                   >  2:950518/1‑131 (MQ=255)
                                                                                                               aaaTCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAg                                <  1:1718664/149‑1 (MQ=255)
                                                                                                                    aGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGc                                                                                                  >  2:1216587/1‑78 (MQ=255)
                                                                                                                    aGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGc                                                                                                  <  1:1216587/78‑1 (MQ=255)
                                                                                                                              tGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTATATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGAcc                 >  1:32747/1‑149 (MQ=255)
                                                                                                                                     gaTGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGt          <  2:32747/149‑1 (MQ=255)
                                                                                                                                       tGCCCTTTGCCGCCGG‑AAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                   >  2:103437/1‑106 (MQ=255)
                                                                                                                                       tGCCCTTTGCCGCCGG‑AAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                   <  1:103437/106‑1 (MQ=255)
                                                                                                                                       tGCCCTTT‑GCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACa             >  2:1295799/1‑144 (MQ=255)
                                                                                                                                       tGCCCTTT‑GCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACa             <  1:1295799/144‑1 (MQ=255)
                                                                                                                                        gCCCTTT‑GCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCAc       >  3:963303/1‑149 (MQ=255)
                                                                                                                                        gCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCAc       >  2:1209891/1‑149 (MQ=255)
                                                                                                                                        gCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCAc       >  4:1737213/1‑149 (MQ=255)
                                                                                                                                          ccTTT‑GCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                   <  2:1535409/103‑1 (MQ=255)
                                                                                                                                          ccTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                   <  2:1616528/103‑1 (MQ=255)
                                                                                                                                          ccTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                   >  1:1616528/1‑103 (MQ=255)
                                                                                                                                           cTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGACTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTgg    >  2:739237/1‑149 (MQ=255)
                                                                                                                                             tt‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTg  <  3:1737213/149‑1 (MQ=255)
                                                                                                                                             tt‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTg  <  1:1209891/149‑1 (MQ=255)
                                                                                                                                                 cggcggCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGtttt                           >  1:929176/1‑121 (MQ=255)
                                                                                                                                                 cggcggCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGtttt                           <  2:929176/121‑1 (MQ=255)
                                                                                                                                                   gcggcAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTa                     >  2:1448839/1‑125 (MQ=255)
                                                                                                                                                   gcggcAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTa                     <  1:1448839/125‑1 (MQ=255)
                                                                                                                                                     ggcAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTa                     <  4:712750/123‑1 (MQ=255)
                                                                                                                                                     ggcAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTa                     >  3:712750/1‑123 (MQ=255)
                                                                                                                                                     ggcAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTgg    >  1:1206748/1‑140 (MQ=255)
                                                                                                                                                     ggcAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTgg    <  2:1206748/140‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                             
CCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTG  >  NC_000913/747752‑748041

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑