Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,642,370 C→G 6.8% intergenic (+303/+119) cspF → / ← quuQ Qin prophage; cold shock protein/Qin prophage; putative antitermination protein Q

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,642,3700CG 6.8% 212.7 / 5.2 80intergenic (+303/+119)cspF/quuQQin prophage; cold shock protein/Qin prophage; putative antitermination protein Q
Reads supporting (aligned to +/- strand):  ref base C (32/40);  new base G (4/2);  total (36/44)
Fisher's exact test for biased strand distribution p-value = 4.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGTGTTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATT  >  NC_000913/1642224‑1642511
                                                                                                                                                  |                                                                                                                                             
ttGTGTTTAGGGACATTGTACTGGAAGAAAAGATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCGGGATGAGATTCAACATTACCATTGCCCCTTTTAAAGCACAAAACCCGAgcatc                                                                                                                                            >  1:298939/1‑149 (MQ=255)
    gtTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAAACAAGTCACCAGATAAATAATAAGTTAACAGAATTGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGc                                                                                                                                         <  1:715306/149‑1 (MQ=255)
     tttAGTGACATTGAACAGTAAGAAAACATTTTAAACATCAGGAAAATACCCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCTTTTACAGCACAAAACCCGCTCATCAGCg                                                                                                                                        <  1:416954/149‑1 (MQ=255)
      ttAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCgg                                                                                                                                       <  2:165479/149‑1 (MQ=255)
        aGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGt                                                                                                                                     <  1:213345/149‑1 (MQ=255)
          ggACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACTGACATGAGTCCCGGGATGAGATTCAACATTACGATTGCCCCATTTAACGCACACAACCCGCTCCTCAGCGGGttt                                                                                                                                   >  2:572585/1‑149 (MQ=255)
          ggACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACACGCTCATCAGCAGGatt                                                                                                                                   >  2:250175/1‑149 (MQ=255)
           gACATTGTACTTGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGtttt                                                                                                                                  <  2:512218/149‑1 (MQ=255)
            aCATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGttttg                                                                                                                                 >  2:755459/1‑148 (MQ=255)
                 gTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACtt                                                                                                                            >  1:568255/1‑149 (MQ=255)
                  tACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACttt                                                                                                                           <  2:945109/149‑1 (MQ=255)
                   aCTGGAAGAAAACATATTAAACAACACGCAAATAACCAACTCACCAGCCAAATAATAACTTAACAGACATGAGTCCCGGGATGCGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACtttt                                                                                                                          <  1:646645/149‑1 (MQ=255)
                    cTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACttttt                                                                                                                         <  2:353559/149‑1 (MQ=255)
                    attgaagaaAACATTTTAAAAATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACttttt                                                                                                                         <  2:224351/146‑1 (MQ=255)
                     tGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTc                                                                                                                        <  2:25833/149‑1 (MQ=255)
                      ggAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATATAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                                                                       <  2:580610/149‑1 (MQ=255)
                       gaagaaAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCCTCAGCGGGTGTTCTACTTGTTCtt                                                                                                                      >  1:251960/1‑149 (MQ=255)
                                aTTTTAACCATCAGGCATATAACCAAGTCAACAACAAAATAATAAGATAACATAAATGAGTACCGGCATGACATTCAACATTACAATTCCCGCATTTAAAGCACAAAACCCGATCATCAGCGGGTATTCCAGTTTTTCTTAACCTCggg                                                                                                             >  2:605885/1‑149 (MQ=255)
                                aTTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCggg                                                                                                             <  2:308543/149‑1 (MQ=255)
                                aTTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTATTCCACTTTTTCTTAACGTCggg                                                                                                             >  2:605884/1‑149 (MQ=255)
                                  tttAAACATCAGGCAAATACCCAAGTCACCAGCTAAATAATAAGTTAACAGCCATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGAGGGTTTTCTACTTTTTCTTAACGTCGGGta                                                                                                           >  2:299851/1‑149 (MQ=255)
                                   ttAAACATCAGCCAAATAACCACGGCACCAGCTAAATAATAAGTTAACAGACCTGAGTCCCGGGATGAGAATCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGTGTTTTCTACTTTTTCTTAACGTCGGGtat                                                                                                          <  2:180955/149‑1 (MQ=255)
                                     aaaCATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATAc                                                                                                        >  2:423553/1‑149 (MQ=255)
                                         aTCAGGCAAATAACCAAGTCAGCAGCTATAGAATATGTTAACAGACAGGAGTCCCGGGAGGCGAGTCCACATTACCATTGCCCCATTTTAAGCACAAAACCCGCTGATCAGCGGGGGGTCTACGTTTTCTTAACGGGGGGTATCCaacg                                                                                                    >  2:561168/1‑147 (MQ=255)
                                           cAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCTTCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGcc                                                                                                  <  2:325584/149‑1 (MQ=255)
                                            aGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGacc                                                                                                 >  2:439815/1‑149 (MQ=255)
                                                 aaTAACCAAGTCACCATCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGt                                                                                            <  1:621201/149‑1 (MQ=255)
                                                 aaTAACCAAGTCACCAGATAAATAAAAAGTTAACAGACACGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAAGGTCGGGTATACAAAGCCCATCGt                                                                                            <  2:266663/149‑1 (MQ=255)
                                                  aTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGtt                                                                                           <  1:596677/149‑1 (MQ=255)
                                                  aTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGtt                                                                                           <  2:769514/149‑1 (MQ=255)
                                                   tAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTg                                                                                          <  1:999264/149‑1 (MQ=255)
                                                              ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTa                                                                               <  2:635420/149‑1 (MQ=255)
                                                                    aaaTAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCata                                                                         >  2:169523/1‑149 (MQ=255)
                                                                        aataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttt                                                                     >  1:67727/1‑149 (MQ=255)
                                                                        aataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGATTTCTAGTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCAtttttt                                                                     >  1:717374/1‑149 (MQ=255)
                                                                        aataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCACCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTATCGTCGGGTATACAAAGCACATCGTTGAAAAAATTTTATACATAtttt                                                                     >  1:613569/1‑149 (MQ=255)
                                                                        aataaGTTAACAGACATGAGTCCCGCGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACACGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGTAAAAATTTTATCCATAtttt                                                                     >  2:501744/1‑149 (MQ=255)
                                                                         ataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttttt                                                                    <  1:979919/149‑1 (MQ=255)
                                                                          taaGTTAACAGACCTGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGAGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttttt                                                                   <  2:361804/149‑1 (MQ=255)
                                                                          taaGTTAACAGACATGAGTCCCGGGATGAGATTCTACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttttt                                                                   <  1:77637/149‑1 (MQ=255)
                                                                          taaGTTAACAGACATGAGGCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttttt                                                                   <  2:824364/149‑1 (MQ=255)
                                                                           aaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTAATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTg                                                                  <  1:761462/149‑1 (MQ=255)
                                                                             gTTAACAGACATGAGTCCCGGGATGAGCTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATTTTTTTTGaa                                                                >  2:224247/1‑149 (MQ=255)
                                                                              ttAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaa                                                               >  1:403425/1‑149 (MQ=255)
                                                                                  cAGACATGAGTCCATGTTTGAGATTCAAAATTAACCTTGCCCCATTTTAAGCACAAAACCCGCGCATCAGAGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCGATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATg                                                           <  2:244047/149‑1 (MQ=255)
                                                                                   ataCATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                                                          <  2:814184/147‑1 (MQ=255)
                                                                                    gACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCa                                                         <  1:71444/149‑1 (MQ=255)
                                                                                      cATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaaa                                                       <  1:321020/149‑1 (MQ=255)
                                                                                       aTGAGTCCAGGGTTGAGATTCAACATTACCATTGCCCCATTTAAAGCACAATCCCCTCTCATCAGCGGGTTTTCTACTTTTTCTTCACGGCGGGTATACAAAGCCCATCGTTGAAAAAATATTATCCATATTTTTTGAAAAATGCAAAc                                                      <  2:239060/149‑1 (MQ=255)
                                                                                         gAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTAGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcaa                                                    >  2:398871/1‑148 (MQ=255)
                                                                                            tCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcatcat                                                 >  2:805889/1‑149 (MQ=255)
                                                                                                  gATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGcc                                           <  1:805889/149‑1 (MQ=255)
                                                                                                  gATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATAGTTGAAAAAATTTTTTCCTATTTTTTTGAAAAATGAAAACAGAATTTCGca                                           >  2:367805/1‑148 (MQ=255)
                                                                                                    tGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATTAAAACATCATGTCGCCAt                                         >  1:326260/1‑149 (MQ=255)
                                                                                                      agaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAATATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCAtgt                                       >  2:99112/1‑147 (MQ=255)
                                                                                                         ttCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAACGCCCATCGATGAAAAAATTTTAAACATATTTTTTGAAAAATACAAACATCATATCTCCAtcatca                                    >  2:855598/1‑149 (MQ=255)
                                                                                                         ttCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCa                                    <  1:436398/149‑1 (MQ=255)
                                                                                                          tCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAg                                   <  1:800095/149‑1 (MQ=255)
                                                                                                          tCAAAATTACCATTGCACCATTTAAAGCACAAAACCCGATCATCAGCGGTTTTTCTACTTTTTCTTAACGTCTGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAg                                   <  1:668031/149‑1 (MQ=255)
                                                                                                            aaCATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGTTATACAAAGCCCATCGTTGAAAAAATTTTATCCATTTTTTTTGAAAAATTCAAACATCATGTCGCCATCTTCAGCa                                 >  1:347369/1‑149 (MQ=255)
                                                                                                               aTTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCaaaa                              >  1:316686/1‑149 (MQ=255)
                                                                                                                    cATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCAtt                        >  2:450984/1‑149 (MQ=255)
                                                                                                                     aTTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCAtt                        <  2:539333/149‑1 (MQ=255)
                                                                                                                      ttGCCCCATTTAAAGCACAAAACCCGCTGATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCAttt                       <  1:896152/149‑1 (MQ=255)
                                                                                                                         ccccATTTAAAGCACAAAACCCTCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTAtc                    <  2:568255/149‑1 (MQ=255)
                                                                                                                         ccccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGTGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTAtc                    <  1:420149/149‑1 (MQ=255)
                                                                                                                          cccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTAtct                   >  1:452674/1‑149 (MQ=255)
                                                                                                                              gaaaaaGTAGAAAACCCGCTGATGAGCGGGTTTTGTGCTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATTTTTTTTGAAAAATGCAAACATAATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTc               >  1:543862/4‑149 (MQ=38)
                                                                                                                               aaaaaGTAGAAAACCCGCTGATGAGCGGGTTTTGTGCTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCa              >  2:270294/3‑149 (MQ=38)
                                                                                                                                  aaGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCt           >  2:224096/1‑149 (MQ=255)
                                                                                                                                    gCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCGATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTGACCTTc         <  2:319781/149‑1 (MQ=255)
                                                                                                                                    gCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCTTCTTCAGCATAAAGCATTTATCTCGTCACCTTc         >  2:303852/1‑149 (MQ=255)
                                                                                                                                    gCACAAAACACTCTCATCATCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGGCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTc         <  1:683804/149‑1 (MQ=255)
                                                                                                                                     tagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTcc        >  1:1025985/4‑149 (MQ=255)
                                                                                                                                      acaAACCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATAGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCt       <  2:251960/149‑1 (MQ=255)
                                                                                                                                       caAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTc      <  2:408724/149‑1 (MQ=255)
                                                                                                                                        aaaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCa     >  2:511076/1‑149 (MQ=255)
                                                                                                                                        aaaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCa     >  1:137997/1‑149 (MQ=255)
                                                                                                                                          aaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAAt   >  1:1021602/1‑149 (MQ=255)
                                                                                                                                          aaCCCGCTAATCAGCGGGTTTTCTACTTTTTCTTACCGTCGGGTATACAAAGCCCATCGTTGAAAAATTTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAAt   <  2:712253/149‑1 (MQ=255)
                                                                                                                                           aCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAAtt  <  1:423553/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
TTGTGTTTAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATT  >  NC_000913/1642224‑1642511

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