Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,541,545 C→G 6.2% H197D (CAC→GAC)  fimB → regulator for fimA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,541,5450CG 6.2% 431.3 / 11.5 127H197D (CAC→GAC) fimBregulator for fimA
Reads supporting (aligned to +/- strand):  ref base C (64/55);  new base G (4/4);  total (68/59)
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.74e-01

Only 100 of 129 total aligned reads displayed.

GGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCAT  >  NC_000913/4541398‑4541694
                                                                                                                                                   |                                                                                                                                                     
     tGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGt                                                                                                                                                >  1:1464072/1‑150 (MQ=255)
      gCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGtt                                                                                                                                               >  1:1088206/1‑150 (MQ=255)
       cTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGttt                                                                                                                                              >  1:826351/1‑150 (MQ=255)
        tttGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTGACGCTGtcgt                                                                                                                                             <  1:630858/150‑4 (MQ=255)
           ggCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTtata                                                                                                                                          >  2:132183/1‑150 (MQ=255)
             cGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATaga                                                                                                                                        <  1:414534/150‑1 (MQ=255)
                    gggAATAGATACGCGACTTATCCAGGATTTTCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATAAATCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGtt                                                                                                                                 >  2:1408124/1‑150 (MQ=255)
                     ggAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTg                                                                                                                                >  1:346695/1‑150 (MQ=255)
                     ggAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTg                                                                                                                                <  2:475000/150‑1 (MQ=255)
                     ggAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTg                                                                                                                                >  2:1275330/1‑150 (MQ=255)
                     ggAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTg                                                                                                                                >  1:803825/1‑150 (MQ=255)
                     ggAATAGATACGCGACTTATCCAGGAGTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGGCACGCTGTTTTATAGAGATAGTTg                                                                                                                                <  2:559076/150‑1 (MQ=255)
                      gAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTga                                                                                                                               <  2:347010/150‑1 (MQ=255)
                           gATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATAc                                                                                                                          >  2:182029/1‑150 (MQ=255)
                               cgcgACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAgg                                                                                                                      >  2:289742/1‑150 (MQ=255)
                               cgcgACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAgg                                                                                                                      >  2:1041348/1‑150 (MQ=255)
                                        tCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGtt                                                                                                             <  1:824139/150‑1 (MQ=255)
                                         ccAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTg                                                                                                                                <  2:39015/130‑1 (MQ=255)
                                         ccAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTg                                                                                                                                >  1:39015/1‑130 (MQ=255)
                                            ggtttATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGTATGGTGTTaaaa                                                                                                         <  1:1307388/147‑1 (MQ=255)
                                              aTTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAAc                                                                                                       <  1:1022223/150‑1 (MQ=255)
                                              aTTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAAc                                                                                                       <  1:1275330/150‑1 (MQ=255)
                                              aTTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAAc                                                                                                       >  2:1332534/1‑150 (MQ=255)
                                              aTTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAAc                                                                                                       <  1:524070/150‑1 (MQ=255)
                                               ttATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATATAGCCTTAGTACGTCAGCGTCACTCTTTTTTTTAGTGATAGTTTATATACCAGGGATGGTGTTaaaaaat                                                                                                      <  2:256468/150‑3 (MQ=255)
                                                       ggCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACattatt                                                                                              >  1:539702/1‑150 (MQ=255)
                                                          aTCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCtt                                                                                           <  1:1240606/150‑1 (MQ=255)
                                                             gCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTAAATTATTCTTACg                                                                                        <  1:346008/150‑1 (MQ=255)
                                                             gCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTATTAAAAACTACATTATTCTTACg                                                                                        >  1:1222386/1‑150 (MQ=255)
                                                                atatTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGata                                                                                     <  1:549752/150‑1 (MQ=255)
                                                                   ttCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGatatat                                                                                  >  1:1558466/1‑150 (MQ=255)
                                                                   ttCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGatatat                                                                                  <  1:1097549/150‑1 (MQ=255)
                                                                        cATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCaaaa                                                                             <  1:1332534/150‑1 (MQ=255)
                                                                        cATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCaaaa                                                                             >  2:735733/1‑150 (MQ=255)
                                                                          tACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAAt                                                                           <  2:826351/150‑1 (MQ=255)
                                                                              gTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGAtt                                                                       >  1:897286/1‑150 (MQ=255)
                                                                               tCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTa                                                                      <  1:289742/150‑1 (MQ=255)
                                                                                cTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTaa                                                                     >  2:199807/1‑150 (MQ=255)
                                                                                cTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTaa                                                                     >  2:756463/1‑150 (MQ=255)
                                                                                 tGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAAc                                                                    <  2:88011/150‑1 (MQ=255)
                                                                                  ggTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAAcc                                                                   >  1:783176/1‑150 (MQ=255)
                                                                                      taCCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCAttt                                                               >  1:1048505/1‑150 (MQ=255)
                                                                                      taCCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCAtt                                                                <  1:29453/149‑1 (MQ=255)
                                                                                        ccGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATttat                                                             >  2:217053/1‑150 (MQ=255)
                                                                                        ccGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATttat                                                             >  1:1342445/1‑150 (MQ=255)
                                                                                        ccGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGAGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATttat                                                             >  2:1452736/1‑150 (MQ=255)
                                                                                          gCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAAt                                                                           <  1:654649/134‑1 (MQ=255)
                                                                                              gCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTg                                                       <  1:505442/150‑1 (MQ=255)
                                                                                              gCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTg                                                       >  1:500095/1‑150 (MQ=255)
                                                                                              gCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTg                                                       <  1:1436175/150‑1 (MQ=255)
                                                                                              gCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTg                                                       >  1:1350914/1‑150 (MQ=255)
                                                                                                aaTGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGAt                                                     <  2:383307/150‑1 (MQ=255)
                                                                                                aaTGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGAt                                                     <  1:132183/150‑1 (MQ=255)
                                                                                                 aTGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATg                                                    <  1:756463/150‑1 (MQ=255)
                                                                                                 aTGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATg                                                    <  1:399934/150‑1 (MQ=255)
                                                                                                   gCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGtt                                                  >  1:1362344/1‑150 (MQ=255)
                                                                                                   gCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGtt                                                  >  2:1255851/1‑150 (MQ=255)
                                                                                                      gggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTAtt                                               >  2:749118/1‑150 (MQ=255)
                                                                                                      gggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTAtt                                               >  1:557829/1‑150 (MQ=255)
                                                                                                       ggCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTg                                              >  1:1149840/1‑150 (MQ=255)
                                                                                                           tttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTaat                                          >  1:703546/1‑149 (MQ=255)
                                                                                                           tttttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATATTATTGTAAg                                          <  1:1533519/150‑1 (MQ=255)
                                                                                                             tttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAgtg                                        >  1:629766/1‑150 (MQ=255)
                                                                                                              ttACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAgtgt                                       <  1:117252/150‑1 (MQ=255)
                                                                                                                aCGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGa                                     >  2:1586416/1‑150 (MQ=255)
                                                                                                                 cGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATAAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGAt                                    >  1:1399158/1‑150 (MQ=255)
                                                                                                                   gCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTc                                  >  1:233111/1‑150 (MQ=255)
                                                                                                                   gCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTc                                  <  2:1464072/150‑1 (MQ=255)
                                                                                                                     aTCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCtt                                >  2:1180438/1‑150 (MQ=255)
                                                                                                                      tCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTa                               <  2:188210/150‑1 (MQ=255)
                                                                                                                       cTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTaa                              >  1:1481464/1‑150 (MQ=255)
                                                                                                                        tGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCt                                 >  1:508908/1‑146 (MQ=255)
                                                                                                                        tGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCt                                 <  2:508908/146‑1 (MQ=255)
                                                                                                                         gggATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAAtt                            <  1:1586416/150‑1 (MQ=255)
                                                                                                                          ggATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTa                           >  1:1239764/1‑150 (MQ=255)
                                                                                                                          ggATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTa                           >  1:562317/1‑150 (MQ=255)
                                                                                                                          ggATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTAATTCTTAATTa                           <  2:897286/150‑1 (MQ=255)
                                                                                                                           gATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTAt                          <  2:151345/150‑1 (MQ=255)
                                                                                                                            aTAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATc                         >  1:538287/1‑150 (MQ=255)
                                                                                                                              agagCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCat                       <  2:470045/150‑1 (MQ=255)
                                                                                                                              agagCAAGAGGACGACAGCGTCACGCTGTTTTATATAGATACTTGAGATACCAGGGATTGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCat                       <  1:373786/150‑1 (MQ=255)
                                                                                                                               gagCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCata                      >  1:747158/1‑150 (MQ=255)
                                                                                                                               gagCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCata                      >  1:929712/1‑150 (MQ=255)
                                                                                                                                 gCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCataat                    >  1:1128635/1‑150 (MQ=255)
                                                                                                                                 gCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCataat                    >  1:153564/1‑150 (MQ=255)
                                                                                                                                  ccAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCataata                   <  1:84641/150‑1 (MQ=255)
                                                                                                                                    agagGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATaaa                 >  1:568057/1‑150 (MQ=255)
                                                                                                                                    agagGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATaaa                 >  2:46064/1‑150 (MQ=255)
                                                                                                                                     ctataaaacAGCGTGACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCa                        <  2:466144/135‑1 (MQ=255)
                                                                                                                                     ctataaaacAGCGTGACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCa                        >  1:466144/8‑142 (MQ=255)
                                                                                                                                     ctataaaacAGCGTGACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAAc                >  1:980607/8‑150 (MQ=255)
                                                                                                                                     ctataaaacAGCGTGACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAAc                >  1:948370/8‑150 (MQ=255)
                                                                                                                                           gacAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAg          <  2:803825/150‑1 (MQ=255)
                                                                                                                                            acAGCGTGACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGt         <  2:948370/150‑1 (MQ=255)
                                                                                                                                             cAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGtt        <  2:539702/150‑1 (MQ=255)
                                                                                                                                                 gTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAcc    <  2:685620/150‑1 (MQ=255)
                                                                                                                                                 gTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAcc    >  2:217659/1‑150 (MQ=255)
                                                                                                                                                 gTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAAcc    <  2:1088206/150‑1 (MQ=255)
                                                                                                                                                  tCACGCTGTTTTATAGAGATAGTTGAGATACCAGTGATTGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCa   <  2:1243223/150‑1 (MQ=255)
                                                                                                                                                   cACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCat  >  2:654881/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGTTTTGCTTTGGCGAATATGGGAATAGATACGCGACTTATCCAGGATTATCTTGGGCATCGCAATATTCGTCATACTGTCTGGTATACCGCCAGCAATGCAGGGCGTTTTTACGGCATCTGGGATAGAGCCAGAGGACGACAGCGTCACGCTGTTTTATAGAGATAGTTGAGATACCAGGGATGGTGTTAAAAACTACATTATTCTTACGATATATCAAAAATGATTAACCATTTATTAATTGATGTTATTGTAAGTGTGATTCTTAATTATCATAATAAACATTAAGTTAACCAT  >  NC_000913/4541398‑4541694

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