Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 618,627 C→A 5.1% P125H (CCT→CAT)  fepE → regulator of length of O‑antigen component of lipopolysaccharide chains

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913618,6270CA 5.1% 481.8 / 9.8 134P125H (CCT→CAT) fepEregulator of length of O‑antigen component of lipopolysaccharide chains
Reads supporting (aligned to +/- strand):  ref base C (62/65);  new base A (4/3);  total (66/68)
Fisher's exact test for biased strand distribution p-value = 7.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

GTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCT  >  NC_000913/618482‑618774
                                                                                                                                                 |                                                                                                                                                   
gTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTa                                                                                                                                                  <  4:600737/149‑1 (MQ=255)
 ttCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTAt                                                                                                                                                 <  2:1352407/149‑1 (MQ=255)
     gTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCt                                                                                                                                                    <  2:132013/142‑1 (MQ=255)
     gTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCt                                                                                                                                                    >  1:132013/1‑142 (MQ=255)
     gTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGa                                                                                                                                             <  4:820728/149‑1 (MQ=255)
        gCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTAtgt                                                                                                                                               >  3:197666/1‑144 (MQ=255)
        gCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTAtgt                                                                                                                                               <  4:197666/144‑1 (MQ=255)
         cAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGa                                                                                                                                         <  2:476708/149‑1 (MQ=255)
                  gagaAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACATAAGGTGATgaac                                                                                                                                        >  2:578658/1‑138 (MQ=255)
                   agaAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAGATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTaaaa                                                                                                                               <  2:255414/149‑1 (MQ=255)
                      aaTCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAagag                                                                                                                            >  1:771960/1‑149 (MQ=255)
                       aTCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAgg                                                                                                                           >  2:1179496/1‑149 (MQ=255)
                              aCTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAAt                                                                                                                    <  1:402829/149‑1 (MQ=255)
                              aCTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAAt                                                                                                                    >  2:607169/1‑149 (MQ=255)
                              aCTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAAt                                                                                                                    >  3:1131296/1‑149 (MQ=255)
                               cTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATc                                                                                                                   <  1:169701/149‑1 (MQ=255)
                                        gtgtGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTg                                                                                                          >  2:990644/1‑149 (MQ=255)
                                            gCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGAtt                                                                                                      <  1:1368726/149‑1 (MQ=255)
                                                 atctggatATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGa                                                                                                                                         >  2:749063/1‑109 (MQ=255)
                                                 atctggatATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTaaaa                                                                                                                               >  2:816657/1‑119 (MQ=255)
                                                 atctggatATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACt                                                                                                           <  1:600020/139‑1 (MQ=255)
                                                 atctggatATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCAt                                                                                                 >  1:1292356/1‑149 (MQ=255)
                                                     ggatATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAagag                                                                                                                            >  1:808933/1‑118 (MQ=255)
                                                          tCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAAtt                                                                                        >  1:997850/1‑149 (MQ=255)
                                                            aaaaTTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGa                                                                                                                                         >  4:709660/1‑98 (MQ=255)
                                                            aaaaTTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATcgac                                                                                                                <  4:502709/123‑1 (MQ=255)
                                                            aaaaTTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATcgac                                                                                                                >  3:502709/1‑123 (MQ=255)
                                                            aaaaTTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACATAAGGTGATgaa                                                                                                                                         >  2:104876/1‑96 (MQ=255)
                                                            aaaaTTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACATAAGGTGATgaa                                                                                                                                         >  2:1160245/1‑96 (MQ=255)
                                                            aaaaTTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGCTCATCACATAAGGTGATgaa                                                                                                                                         <  1:1160245/98‑3 (MQ=255)
                                                               aTTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCCATTAAAAGAgg                                                                                                                           <  3:837916/109‑1 (MQ=255)
                                                               aTTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGc                                                                                   <  1:246813/149‑1 (MQ=255)
                                                                   aTCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTg                                                                               >  1:230259/1‑149 (MQ=255)
                                                                    tCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGa                                                                              <  3:77019/149‑1 (MQ=255)
                                                                      gTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGt                                                                                      <  1:446314/139‑1 (MQ=255)
                                                                      gTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGt                                                                                      >  2:446314/1‑139 (MQ=255)
                                                                            aaGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATc                                                                                                <  2:1209266/123‑1 (MQ=255)
                                                                            aaGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATc                                                                                                >  1:1209266/1‑123 (MQ=255)
                                                                                aTTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTAt                                                                                                                                                 <  1:1015530/70‑1 (MQ=255)
                                                                                  ttAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGaaa                                                                <  1:990644/149‑1 (MQ=255)
                                                                                      ccTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGaa                                                                 >  1:1203029/1‑144 (MQ=255)
                                                                                       cTGTTTATCAAGAAGTTTCAGTCGGTTATCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAgg                                                                                                                           <  2:599537/85‑1 (MQ=255)
                                                                                       cTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAgg                                                                                                                           >  1:599537/1‑85 (MQ=255)
                                                                                       cTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGt                                                           >  3:1202955/1‑149 (MQ=255)
                                                                                        tGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGAc                                                                                                                                        <  3:984216/71‑1 (MQ=255)
                                                                                        tGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCAtt                                                                                <  4:730520/127‑1 (MQ=255)
                                                                                        tGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCAtt                                                                                >  3:730520/1‑127 (MQ=255)
                                                                                        tGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGtt                                                          >  3:348656/1‑149 (MQ=255)
                                                                                           ttATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTtgat                                                       <  2:230259/149‑1 (MQ=255)
                                                                                           ttATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTtgat                                                       <  4:1202955/149‑1 (MQ=255)
                                                                                                aagaagTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACaa                                                  >  2:795998/1‑149 (MQ=255)
                                                                                                aagaagTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACaa                                                  >  4:937284/1‑149 (MQ=255)
                                                                                                 agaagTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAAt                                                 >  3:1034974/1‑149 (MQ=255)
                                                                                                   aagTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGc                                               <  4:1034974/149‑1 (MQ=255)
                                                                                                     gTTTAAGTCCGTTAGCTTCCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCCTTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCa                                             <  4:1131296/149‑1 (MQ=255)
                                                                                                      tttCAGTCGGTTAGCTTGCTGGAAGAGTGCCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACa                                                   <  1:675170/142‑1 (MQ=255)
                                                                                                         cAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTaa                                         <  4:348656/149‑1 (MQ=255)
                                                                                                               gTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTc                                                                                     >  4:328291/1‑99 (MQ=255)
                                                                                                               gTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTc                                                                                     <  3:328291/99‑1 (MQ=255)
                                                                                                                ttAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGaaaaaa                                  <  4:79551/149‑1 (MQ=255)
                                                                                                                ttAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGaaaaaa                                  <  1:795998/149‑1 (MQ=255)
                                                                                                                  aGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGAAAAAAGa                                <  1:486312/149‑1 (MQ=255)
                                                                                                                  aGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGAAAAAAGa                                <  3:937284/149‑1 (MQ=255)
                                                                                                                    cTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGc                                                                            >  4:481065/1‑103 (MQ=255)
                                                                                                                    cTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGc                                                                            <  3:481065/103‑1 (MQ=255)
                                                                                                                    cTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGtt                                                          <  4:961994/121‑1 (MQ=255)
                                                                                                                     ttGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTtga                                                        <  1:1095892/122‑1 (MQ=255)
                                                                                                                      tGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTtgatg                                                      <  3:1453074/123‑1 (MQ=255)
                                                                                                                      tGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTtgatg                                                      >  4:1453074/1‑123 (MQ=255)
                                                                                                                       gCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGAAAAAAGATGAAc                           <  2:771960/149‑1 (MQ=255)
                                                                                                                         tGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGaaaaaa                                  <  2:1025644/140‑1 (MQ=255)
                                                                                                                         tGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGaaaaaa                                  >  1:1025644/1‑140 (MQ=255)
                                                                                                                          ggAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGaaaaaa                                  <  1:1227521/139‑1 (MQ=255)
                                                                                                                          ggAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGaaaaaa                                  >  2:1227521/1‑139 (MQ=255)
                                                                                                                           gAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCa                                             >  1:1016512/1‑127 (MQ=255)
                                                                                                                           gAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCa                                             <  2:1016512/127‑1 (MQ=255)
                                                                                                                            aaGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGa                                                                              >  1:611170/1‑93 (MQ=255)
                                                                                                                            aaGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGa                                                                              <  2:611170/93‑1 (MQ=255)
                                                                                                                              gagTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAAt                                                 >  2:343109/1‑120 (MQ=255)
                                                                                                                              gagTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGaaaaaa                                  <  1:113199/135‑1 (MQ=255)
                                                                                                                              gagTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGaaaaaa                                  >  2:113199/1‑135 (MQ=255)
                                                                                                                                  aCCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGAAAAAAGATGAACCGTCACTGtat                <  3:706164/149‑1 (MQ=255)
                                                                                                                                   ccTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGtt                                                          <  3:818377/106‑1 (MQ=255)
                                                                                                                                   ccTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGtt                                                          >  4:818377/1‑106 (MQ=255)
                                                                                                                                   ccTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGtt                                                          >  2:93783/1‑106 (MQ=255)
                                                                                                                                   ccTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGtt                                                          <  1:93783/106‑1 (MQ=255)
                                                                                                                                     tGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTtgat                                                       >  3:1046557/1‑107 (MQ=255)
                                                                                                                                        gtccatcaCATAAGGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGAAAAAAGATGa                             >  3:1462965/4‑130 (MQ=255)
                                                                                                                                        gtccatcaCATAAGGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGAAAAAAGATGa                             <  4:1462965/127‑1 (MQ=255)
                                                                                                                                            atcaCCTTATGTGATGGACCAATTAAAAGAGGTGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGAAAAAAGATGAACCGt                        <  3:723660/131‑1 (MQ=255)
                                                                                                                                            atcaCCTTATGTGATGGACCAATTAAAAGAGGTGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGAAAAAAGATGAACCGt                        >  4:723660/1‑131 (MQ=255)
                                                                                                                                            atcaCCTTATGTGATGGACCAATTAAAAGAGGCGAAAATcga                                                                                                                 >  4:573563/1‑42 (MQ=255)
                                                                                                                                            atcaCCTTATGTGATGGACCAATTAAAAGAGGCGAAAATcga                                                                                                                 <  3:573563/42‑1 (MQ=255)
                                                                                                                                             tcaCCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAAt                                                                    <  2:1232494/86‑1 (MQ=255)
                                                                                                                                             tcaCCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGt                                                           <  2:659109/95‑1 (MQ=255)
                                                                                                                                             tcaCCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGt                                                           >  1:659109/1‑95 (MQ=255)
                                                                                                                                             tcaCCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAAt                                                 >  4:1202231/1‑105 (MQ=255)
                                                                                                                                             tcaCCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAAt                                                 <  3:1202231/105‑1 (MQ=255)
                                                                                                                                             tcaCCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGAAAAAAGATGAACCGTCACTGTATACCTCCTGGAc     >  2:494092/1‑149 (MQ=255)
                                                                                                                                             tcaCCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTCGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGAAAAAAGATGAACCGTCACTGTATACCTCCTGGAc     >  2:425144/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GTTCAGTGGCAAGAGTTGGAGAAATCATTCACTAAGCTTCGTGTGCTGGATCTGGATATCAAAATTGATCGTACAGAAGCATTTAACCTGTTTATCAAGAAGTTTCAGTCGGTTAGCTTGCTGGAAGAGTACCTGCGTTCATCACCTTATGTGATGGACCAATTAAAAGAGGCGAAAATCGACGAACTGGATTTGCATCGCGCAATTGTCGCATTGAGCGAAAAAATGAAAGCGGTTGATGACAATGCCAGTAAGAAAAAAGATGAACCGTCACTGTATACCTCCTGGACGCT  >  NC_000913/618482‑618774

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: