Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,642,373 C→G 6.7% intergenic (+306/+116) 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,3730CG 6.7% 370.4 / 14.6 119intergenic (+306/+116)cspF/quuQQin prophage; cold shock protein/Qin prophage; putative antitermination protein Q
Reads supporting (aligned to +/- strand):  ref base C (57/54);  new base G (6/2);  total (63/56)
Fisher's exact test for biased strand distribution p-value = 2.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

Only 100 of 121 total aligned reads displayed.

TAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTTCCG  >  NC_000913/1642231‑1642520
                                                                                                                                              |                                                                                                                                                   
tAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGt                                                                                                                                              <  1:929025/150‑1 (MQ=255)
tAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTAACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTAACATTGCCCCAATTAAAGCACAAAACCCGCTCATCAGCGGGt                                                                                                                                              >  2:654652/1‑150 (MQ=255)
  gggACATTGTACTTGAAGAAATCATTTTAAACATCAGGCAAATAACCAAGTCACCATCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttt                                                                                                                                            <  1:884619/150‑1 (MQ=255)
  gggACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCAACAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttt                                                                                                                                            <  1:1089602/150‑1 (MQ=255)
   ggACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGtttt                                                                                                                                           >  1:279239/1‑150 (MQ=255)
    gACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTc                                                                                                                                          <  2:27662/150‑1 (MQ=255)
      cATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGCCATGAGTCACGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTa                                                                                                                                        >  1:121171/1‑150 (MQ=255)
      cATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTa                                                                                                                                        >  2:1097719/1‑150 (MQ=255)
       aTTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTAc                                                                                                                                       <  1:658911/150‑1 (MQ=255)
         tGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACtt                                                                                                                                     <  2:797365/150‑1 (MQ=255)
            aCTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACttttt                                                                                                                                  <  1:1073219/150‑1 (MQ=255)
             cTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTc                                                                                                                                 >  2:349696/1‑150 (MQ=255)
             cTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTc                                                                                                                                 >  2:1545141/1‑150 (MQ=255)
              tGGAAGAAAACATTTTAAACATCAGTCAAATAAACAAGTCACCAGCTAAATAAAAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                                                                                <  2:7088/150‑1 (MQ=255)
              tGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                                                                                <  1:1097719/150‑1 (MQ=255)
               ggAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCtt                                                                                                                               >  1:1185510/1‑150 (MQ=255)
                  agaaAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAAc                                                                                                                            >  1:1514270/1‑150 (MQ=255)
                  agaaAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTAAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAAc                                                                                                                            >  1:940274/1‑150 (MQ=255)
                   gaaAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCGCGGGATGAGATTCAACATTACCATTGCCCCATTTAAATCACAAAACCCGCGCATCAGCGGGTTTTCTACTTTTTCTTAACg                                                                                                                           >  1:252163/1‑150 (MQ=255)
                    aaaaCATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGt                                                                                                                          >  2:774243/1‑150 (MQ=255)
                      aaCATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCAATTAAAGCACAAACCCCGCCCATCAGCGGGTTTTCTACTTTTTCTTAACGTCg                                                                                                                        <  1:256638/150‑1 (MQ=255)
                        cATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCggg                                                                                                                      <  1:33992/150‑1 (MQ=255)
                        cATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCggg                                                                                                                      <  2:727980/150‑1 (MQ=255)
                        cATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCggg                                                                                                                      <  2:578335/150‑1 (MQ=255)
                        aaTTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCggg                                                                                                                      <  2:1107376/149‑1 (MQ=255)
                        aaTTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCggg                                                                                                                      <  2:628449/149‑1 (MQ=255)
                          ttttAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGta                                                                                                                    >  1:611464/1‑150 (MQ=255)
                             tAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATAc                                                                                                                 >  2:717986/1‑150 (MQ=255)
                                   tCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGcc                                                                                                           <  1:387461/150‑1 (MQ=255)
                                    cAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGccc                                                                                                          >  1:1043157/1‑150 (MQ=255)
                                     aGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCATAACGTCGGGTATACAAAGCCCa                                                                                                         >  2:1609927/1‑150 (MQ=255)
                                     aGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTATTAACGTCGGGTATACAAAGCCCa                                                                                                         >  1:325538/1‑150 (MQ=255)
                                     aGGCAAATAACCAAGTCAACAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCAGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGACGGGTATACAAAGCCCa                                                                                                         >  1:1415356/1‑150 (MQ=255)
                                      ggCAAATAACCAAGTAACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGAACAAAACCCGCTCAGCAGCGGGTTTTCTTCTTTTTCTTAACGTCGGGTATAAATAGGCCAt                                                                                                        >  2:281664/1‑150 (MQ=255)
                                        cAAATAACCAAGTCACCAGCGAAATAATAGGTTAGGAGACATGAGTCCGGGCAGGAGTGTCAACTATACCAATGACCCATTTAAAGCACAAAACCCGCTCAACAGCGGGTTTGCTTCTTTTTCTTAAAGTCGGGATTACAAAGACCAtag                                                                                                      >  2:307106/1‑148 (MQ=255)
                                          aaTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGtt                                                                                                    <  1:491015/150‑1 (MQ=255)
                                           aTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTg                                                                                                   <  1:478296/150‑1 (MQ=255)
                                            tAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACAATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGa                                                                                                  <  2:339537/150‑1 (MQ=255)
                                                cAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaaaaa                                                                                              >  2:124057/1‑150 (MQ=255)
                                                   gTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtt                                                                                           <  1:349696/150‑1 (MQ=255)
                                                     cACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtttt                                                                                         >  2:933563/1‑150 (MQ=255)
                                                       ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTAt                                                                                       <  1:1433648/150‑1 (MQ=255)
                                                       ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTAt                                                                                       >  2:610740/1‑150 (MQ=255)
                                                         aGCTAAATAATAAGTTACCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCTCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATcc                                                                                     <  2:272560/150‑1 (MQ=255)
                                                         aGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATcc                                                                                     >  1:50256/1‑150 (MQ=255)
                                                         aGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATcc                                                                                     >  2:1659181/1‑150 (MQ=255)
                                                         aGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATcc                                                                                     >  1:1534653/1‑150 (MQ=255)
                                                            tAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCata                                                                                  >  2:778018/1‑150 (MQ=255)
                                                             aaaTAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtata                                                                                                                  <  1:1287893/117‑1 (MQ=255)
                                                             aaaTAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtata                                                                                                                  >  2:1287893/1‑117 (MQ=255)
                                                                 aataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttttt                                                                             >  2:1440530/1‑150 (MQ=255)
                                                                   tatgTTAACAGACATGAGTCCCGGTATGAGATTCAAAATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTAATTTTTCTTAACGTGGGTTTTACAAAGCCCTTGGTTTAAAAAATTTTATCCATATTTTTTg                                                                           <  2:1335462/147‑1 (MQ=255)
                                                                   taaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTg                                                                           <  1:575443/150‑1 (MQ=255)
                                                                   taaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTg                                                                           <  1:726299/150‑1 (MQ=255)
                                                                         aaCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAAt                                                                     >  2:235739/1‑150 (MQ=255)
                                                                         aaCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAAt                                                                     <  2:1715757/150‑1 (MQ=255)
                                                                           cAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCTGGTTTACAAATCCCATCGTTGAAAAAATTTTATACATATTTTTTGAAAAATGc                                                                   >  2:508438/1‑150 (MQ=255)
                                                                           cAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                                                                   <  2:534375/150‑1 (MQ=255)
                                                                             gACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaa                                                                 <  2:703680/150‑1 (MQ=255)
                                                                              aCATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAATCACAAAACCCGCTCATCAGCGTGTTTTCTACTTTTTCTTAACGTCTGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaaa                                                                <  1:332980/150‑1 (MQ=255)
                                                                               cATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAc                                                               >  1:788658/1‑150 (MQ=255)
                                                                                aTGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAca                                                              <  1:1417790/150‑1 (MQ=255)
                                                                                   aGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTAATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcatca                                                           <  1:1099356/150‑1 (MQ=255)
                                                                                    gTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcatcat                                                          >  1:817149/1‑150 (MQ=255)
                                                                                      cccGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGt                                                        >  1:515183/1‑150 (MQ=255)
                                                                                      cccGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGt                                                        >  1:480868/1‑150 (MQ=255)
                                                                                            aTGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAAATCATGTCGCCAt                                                  >  2:28794/1‑150 (MQ=255)
                                                                                               agaTTCAACATTTCCATTGCCACATTTAAAGCACATAACGCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAAATCATGTCTCGATCtt                                               <  2:1418879/150‑1 (MQ=255)
                                                                                                  ttCAACATTACCATTGCCCCATTAAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAg                                            <  2:1124391/150‑1 (MQ=255)
                                                                                                    cAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCa                                          <  1:778018/150‑1 (MQ=255)
                                                                                                    cAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCa                                          >  2:1028506/1‑150 (MQ=255)
                                                                                                         ttACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAAt                                     >  1:1438253/1‑150 (MQ=255)
                                                                                                         ttACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAAt                                     <  2:325538/150‑1 (MQ=255)
                                                                                                          tACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATc                                    <  2:1473212/150‑1 (MQ=255)
                                                                                                           aCCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAATAAATTATATAAATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCa                                   >  1:240856/1‑150 (MQ=255)
                                                                                                            ccATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCAt                                  <  2:515183/150‑1 (MQ=255)
                                                                                                              aTTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCAttt                                <  2:611464/150‑1 (MQ=255)
                                                                                                                 gCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTAtc                             <  2:1185510/150‑1 (MQ=255)
                                                                                                                 gCCCCATTTAAAGCACAAAACCCGCTCATCAGCAGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTAtc                             <  2:664478/150‑1 (MQ=255)
                                                                                                                  ccccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTAtct                            >  1:762814/1‑150 (MQ=255)
                                                                                                                   cccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTAtctc                           >  1:1594354/1‑150 (MQ=255)
                                                                                                                     cATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGt                         >  2:1030608/1‑150 (MQ=255)
                                                                                                                       tttAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCTCCATCTTCAGCAAAAATCATTTATCTCGTCa                       <  1:654652/150‑1 (MQ=255)
                                                                                                                       gaaaaaGTAGAAAACCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCa                       >  1:788653/4‑150 (MQ=255)
                                                                                                                         aaaaGTAGAAAACCCGCTGATGAGCGGGTTTTGTGCTTTTTCTTAACGTCGGGTAGACAAAGCACATCGGTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTGCCGATATTCAGCAAAAATATATTATCTCgaaacc                     >  2:326244/2‑145 (MQ=25)
                                                                                                                          aaaGCACAAAACCTGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGAAAACATCATGTCTCCATCTTCAGCAAAAATCATTTAACTCGTCACCt                    >  2:244819/1‑150 (MQ=255)
                                                                                                                           aaGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCtt                   >  1:164502/1‑150 (MQ=255)
                                                                                                                             gCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATGGTTGAATAAATTTTATCCATATTTTTTGAAAAATTCAAACATCATGTCTCCATCTTCAGCAAAAATCATTTTTCTGGTGAcatgcc                 >  2:1054957/1‑145 (MQ=255)
                                                                                                                               acaAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAACGCCCCTCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGAAAACATCATTTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTc               >  1:751492/1‑150 (MQ=255)
                                                                                                                               acaAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcatc                                                            <  2:1162725/105‑1 (MQ=255)
                                                                                                                               acaAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcatc                                                            >  1:1162725/1‑105 (MQ=255)
                                                                                                                                     cccGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTgc         <  1:684329/150‑1 (MQ=255)
                                                                                                                                      ccGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTg          <  2:1431652/148‑1 (MQ=255)
                                                                                                                                      ccGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTg          >  1:1431652/1‑148 (MQ=255)
                                                                                                                                      ccGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTg          >  1:1245084/1‑148 (MQ=255)
                                                                                                                                      ccGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTg          <  2:1245084/148‑1 (MQ=255)
                                                                                                                                      ccGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCATCAATTgcg        <  2:546900/150‑1 (MQ=255)
                                                                                                                                       cGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTgcgc       <  1:1545141/150‑1 (MQ=255)
                                                                                                                                           catcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAATCCCATCGTTGAAAAAATTTTATACATATTTTTTTAAAAATTCAAACATAATGACGCAATCTACACCAAAAATCATATTTATCGACAACTTCCTCAATTTCGATTcc   >  1:540767/1‑150 (MQ=255)
                                                                                                                                            atcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTTCcg  >  2:684692/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TAGGGACATTGTACTGGAAGAAAACATTTTAAACATCAGGCAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCATCTTCAGCAAAAATCATTTATCTCGTCACCTTCCTCAATTGCGCTTCCG  >  NC_000913/1642231‑1642520

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