Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 747,897 G→C 9.9% intergenic (‑128/+24) 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,8970GC 9.9% 312.5 / 22.0 91intergenic (‑128/+24)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.41e-01

CCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTG  >  NC_000913/747751‑748041
                                                                                                                                                   |                                                                                                                                                
cccACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑Ccggc                                                                                                                                                >  4:1865766/1‑149 (MQ=255)
 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                                                                                                                                         <  3:1161464/115‑1 (MQ=255)
                                         tctTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAg                                                                                                                                         >  4:1161464/1‑115 (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                                                                                                 >  1:751804/1‑133 (MQ=255)
                                                               tAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                 <  2:751804/133‑1 (MQ=255)
                                                               tAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAggg                                                                                 <  3:1865766/149‑1 (MQ=255)
                                                                    aaagCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAg                                                                            <  4:1478586/149‑1 (MQ=255)
                                                                         tGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCa                                                                                                                        >  2:304428/1‑100 (MQ=255)
                                                                         tGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCa                                                                                                                        <  1:304428/100‑1 (MQ=255)
                                                                             ccAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAg                                                                                                   <  4:480509/117‑1 (MQ=255)
                                                                             ccAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAg                                                                                                   >  3:480509/1‑117 (MQ=255)
                                                                             ccAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgg                                                                   <  2:1849969/149‑1 (MQ=255)
                                                                                 gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcac                                                                                                >  4:324776/1‑116 (MQ=255)
                                                                                 gAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGcac                                                                                                <  3:324776/116‑1 (MQ=255)
                                                                                      aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCa                                                                                                                        >  3:1295548/1‑87 (MQ=255)
                                                                                      aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCa                                                                                                                        <  4:1295548/87‑1 (MQ=255)
                                                                                      aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTtata                                                                                                                  <  2:839314/93‑1 (MQ=255)
                                                                                      aaTAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTtata                                                                                                                  >  1:839314/1‑93 (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                                                                                                                                        >  1:1419948/1‑62 (MQ=255)
                                                                                               aaTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAgg                                                                                                                                        <  2:1419948/62‑1 (MQ=255)
                                                                                               aaTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                   >  1:282718/1‑147 (MQ=255)
                                                                                               aaTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                   <  2:282718/147‑1 (MQ=255)
                                                                                                aTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgggg                                                                 <  4:970/132‑1 (MQ=255)
                                                                                                aTCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATgggg                                                                 >  3:970/1‑132 (MQ=255)
                                                                                                 tCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaa                                               >  4:419159/1‑149 (MQ=255)
                                                                                                  cAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaaa                                              <  3:419159/149‑1 (MQ=255)
                                                                                                  cAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGaaa                                              <  3:1622346/149‑1 (MQ=255)
                                                                                                    aGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                 >  4:872349/1‑96 (MQ=255)
                                                                                                    aGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                 >  2:1174312/1‑96 (MQ=255)
                                                                                                    aGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                 <  3:872349/96‑1 (MQ=255)
                                                                                                    aGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGca                                                                                                 <  1:1174312/96‑1 (MQ=255)
                                                                                                      gAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTc                                          >  2:1718664/1‑149 (MQ=255)
                                                                                                         aaGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                           <  2:440628/81‑1 (MQ=255)
                                                                                                         aaGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATc                                                                                                           >  1:440628/1‑81 (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                                                                                                  <  1:1216587/78‑1 (MQ=255)
                                                                                                                     aGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGc                                                                                                  >  2:1216587/1‑78 (MQ=255)
                                                                                                                               tGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTATATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGAcc                 >  1:32747/1‑149 (MQ=255)
                                                                                                                                      gaTGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGt          <  2:32747/149‑1 (MQ=255)
                                                                                                                                        tGCCCTTTGCCGCCGG‑AAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                   <  1:103437/106‑1 (MQ=255)
                                                                                                                                        tGCCCTTTGCCGCCGG‑AAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                   >  2:103437/1‑106 (MQ=255)
                                                                                                                                        tGCCCTTT‑GCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACa             <  1:1295799/144‑1 (MQ=255)
                                                                                                                                        tGCCCTTT‑GCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACa             >  2:1295799/1‑144 (MQ=255)
                                                                                                                                         gCCCTTT‑GCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCAc       >  3:963303/1‑149 (MQ=255)
                                                                                                                                         gCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCAc       >  4:1737213/1‑149 (MQ=255)
                                                                                                                                         gCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCAc       >  2:1209891/1‑149 (MQ=255)
                                                                                                                                           ccTTT‑GCCGCCGGAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                   <  2:1535409/103‑1 (MQ=255)
                                                                                                                                           ccTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                   >  1:1616528/1‑103 (MQ=255)
                                                                                                                                           ccTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGa                                                   <  2:1616528/103‑1 (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)
                                                                                                                                                   |                                                                                                                                                
CCCACTGCAAAACTGGCATCCTTCGCTCCCTCATTATTTTCTCTTTGGTGTATCGATACAGCATAAAGAAAGCTGGGCCAGGAAAAAATAACAGCAATCAAGGAAAAGGGGAAAATCAGCAATTTTCTGAAAGAGATGCCCTTT‑CCGGCGGCAAAGGGCATCATGGTAAATCAGTTATAGTTAATCTTAAAAAGCACCACTGTATCGAAGGGGCCAGCCTCGATGGGGGTGCCTGGAATGGACGAAAGTTCAGCCGTATAGGTTTTTTCGTAGACCAACATGTCACTGGTG  >  NC_000913/747751‑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: ‑