Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,207,908 (G)5→6 coding (408/630 nt) ycfK → e14 prophage; protein StfP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,207,9031.G77.3% 183.8 / 49.3 66coding (403/630 nt)ycfKe14 prophage; protein StfP
Reads supporting (aligned to +/- strand):  ref base . (9/6);  new base G (30/21);  total (39/27)
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.65e-01

CACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGT  >  NC_000913/1207764‑1208036
                                                                                                                                            |                                                                                                                                       
ttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGt                                                                                                                                          >  1:770432‑M1/27‑140 (MQ=255)
ttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGt                                                                                                                                          >  2:93572‑M1/27‑140 (MQ=255)
ttgttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGt                                                                                                                                          >  2:4003‑M1/27‑140 (MQ=255)
   ttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTggg                                                                                                                                      >  1:142790‑M1/24‑140 (MQ=255)
   ttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTggg                                                                                                                                      >  2:163093‑M1/24‑140 (MQ=255)
   ttatttcttcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTggg                                                                                                                                      >  1:325310‑M1/24‑140 (MQ=255)
           tcattttcataaacgTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGaa                                                                                                                             >  2:78721‑M1/16‑140 (MQ=255)
                        cgTTAGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCa                                                                                                                <  1:533611‑M1/138‑1 (MQ=255)
                         gTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCAt                                                                                                               >  2:409491‑M1/2‑140 (MQ=255)
                         gTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCAt                                                                                                               <  2:244615‑M1/139‑1 (MQ=255)
                          TTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATc                                                                                                              >  2:397216‑M1/1‑140 (MQ=255)
                           TGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCt                                                                                                             <  1:403971‑M1/140‑1 (MQ=255)
                                  GGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTg                                                                                                      <  1:385480‑M1/140‑1 (MQ=255)
                                   GAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgg                                                                                                     <  2:153946‑M1/140‑1 (MQ=255)
                                    AGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggg                                                                                                    >  2:408836‑M1/1‑140 (MQ=255)
                                     GAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgggg                                                                                                   >  2:888150‑M1/1‑140 (MQ=255)
                                     GAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgggg                                                                                                   >  1:874787‑M1/1‑140 (MQ=255)
                                       AGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgggg                                                                                                   <  1:654675‑M1/138‑1 (MQ=255)
                                       AGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTgggg                                                                                                   >  2:654675‑M1/1‑138 (MQ=255)
                                        GGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTggggggg                                                                                                >  2:762403‑M1/1‑140 (MQ=255)
                                                   cAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACa                                                                                     <  2:301385/140‑1 (MQ=255)
                                                    aaTGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAg                                                                                    <  1:9167/140‑1 (MQ=255)
                                                           caACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTAtt                                                                             >  1:540248/1‑140 (MQ=255)
                                                             aCGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTcc                                                                           <  1:413067/140‑1 (MQ=255)
                                                                 cTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGtt                                                                      <  1:784490/140‑1 (MQ=255)
                                                                  tACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGtt                                                                      <  1:711907/140‑1 (MQ=255)
                                                                     aTTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCg                                                                   <  1:52860/140‑1 (MQ=255)
                                                                      ttAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGt                                                                  <  2:3388/140‑1 (MQ=255)
                                                                           aCGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaat                                                             >  1:828478/1‑140 (MQ=255)
                                                                            cGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTtaata                                                            >  1:559493/1‑140 (MQ=255)
                                                                              cACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTaataat                                                          >  2:482979/1‑140 (MQ=255)
                                                                              cACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTaataat                                                          >  2:406870/1‑140 (MQ=255)
                                                                               aCCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATg                                                         >  1:324811/1‑140 (MQ=255)
                                                                                   gTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCtg                                                     >  1:160192/1‑140 (MQ=255)
                                                                                      tCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTgtgt                                                  >  2:348885/1‑140 (MQ=255)
                                                                                         tAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGc                                               >  2:765213/1‑140 (MQ=255)
                                                                                                 ccGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGttgttg                                       <  1:265821/140‑1 (MQ=255)
                                                                                                 ccGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGttgttg                                       <  1:30090/140‑1 (MQ=255)
                                                                                                 ccGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGttgttg                                       <  2:533611/140‑1 (MQ=255)
                                                                                                 ccGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGttgttg                                       <  2:794473/140‑1 (MQ=255)
                                                                                                    gTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCaa                                    >  1:524780/1‑140 (MQ=255)
                                                                                                      tATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAAt                                  >  1:871442/1‑140 (MQ=255)
                                                                                                         gATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTAt                               >  1:630532/1‑140 (MQ=255)
                                                                                                            tggtggTATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGAc                            >  2:35591/1‑140 (MQ=255)
                                                                                                              gtggtATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACaa                          <  2:575904/140‑1 (MQ=255)
                                                                                                              gtggtATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACaa                          <  2:296715/140‑1 (MQ=255)
                                                                                                              gtggtATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACaa                          >  2:859840/1‑140 (MQ=255)
                                                                                                              gtggtATTGAGCAATCAGTAATCATTCAGTGGGGAGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGGGTTGCCGTTGTTGCAAATTACGACaa                          >  2:816650/1‑140 (MQ=255)
                                                                                                               tggtATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAAt                         >  1:195286/1‑140 (MQ=255)
                                                                                                               tggtATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAAt                         >  2:490022/1‑140 (MQ=255)
                                                                                                               tggtATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAAt                         >  2:501775/1‑140 (MQ=255)
                                                                                                                  tATTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTc                      >  1:862462/1‑140 (MQ=255)
                                                                                                                   aTTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCa                     >  1:478626/1‑140 (MQ=255)
                                                                                                                   aTTGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCa                     >  1:552200/1‑140 (MQ=255)
                                                                                                                    ttGAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAg                    >  2:631937/1‑140 (MQ=255)
                                                                                                                    ttGAGCAATCAGTAATCATTCAGTGGGGGGCGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAg                    >  2:517740/1‑140 (MQ=255)
                                                                                                                      gAGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgcg                  <  1:424849/140‑1 (MQ=255)
                                                                                                                       aGCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAgcgc                 >  1:813571/1‑140 (MQ=255)
                                                                                                                        gCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCa                <  1:153946/140‑1 (MQ=255)
                                                                                                                        gCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCa                <  2:540248/140‑1 (MQ=255)
                                                                                                                        gCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCa                <  2:524780/140‑1 (MQ=255)
                                                                                                                        gCAATCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCa                <  1:222594/140‑1 (MQ=255)
                                                                                                                            tCAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTa            <  2:828478/140‑1 (MQ=255)
                                                                                                                             cAGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTAt           <  1:501775/140‑1 (MQ=255)
                                                                                                                              aGTAATCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATc          >  1:25290/1‑140 (MQ=255)
                                                                                                                                 aaTCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGt                       >  2:684040/1‑124 (MQ=255)
                                                                                                                                 aaTCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGt                       <  1:684040/124‑1 (MQ=255)
                                                                                                                                 aaTCATTCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAAt       >  1:704384/1‑140 (MQ=255)
                                                                                                                                      ttCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGt  >  2:460307/1‑140 (MQ=255)
                                                                                                                                      ttCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGt  >  2:255896/1‑140 (MQ=255)
                                                                                                                                      ttCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGt  >  1:94648/1‑140 (MQ=255)
                                                                                                                                      ttCAGTGGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGt  >  1:590781/1‑140 (MQ=255)
                                                                                                                                            |                                                                                                                                       
CACTGTGCAAACGGCTCTCGAAAACCTTGGTTTGGGAGAAGGAGCAAAACTCAATGCAGCAACGGCTACATTAGGACGCACCGGTTTCATAGCTATACCGGTTATGATTGGTGGTATTGAGCAATCAGTAATCATTCAGTGGGGGTGGAATGCCGCAAAAGCATCTGCCTCTGGGGGGGATGGAAATACAGTTGTATTCCCGGTTGCGTTTAATAATGCCTGTGTTGCCGTTGTTGCAAATTATGACAATGTCAGCGCACCTATCAATGCAGT  >  NC_000913/1207764‑1208036

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: