Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,642,373 C→G 17.0% 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,3730CG17.0% 319.4 / 57.7 118intergenic (+306/+116)cspF/quuQQin prophage; cold shock protein/Qin prophage; putative antitermination protein Q
Reads supporting (aligned to +/- strand):  ref base C (52/46);  new base G (9/11);  total (61/57)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.56e-01

Only 100 of 125 total aligned reads displayed.

ATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGC  >  NC_000913/1642274‑1642469
                                                                                                   |                                                                                                
aTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatc                                                                                                  >  2:1094179/1‑100 (MQ=255)
 tAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatca                                                                                                 >  2:1472435/1‑100 (MQ=255)
  aaCCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCACtcatca                                                                                                 >  2:2873440/1‑99 (MQ=255)
   aCCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGc                                                                                               >  2:943427/1‑100 (MQ=255)
       aGTCACCAGCTAAATGATAAGTTAACATACATGAGTCCCGGGATGAGATTCAACATTACCATGGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGt                                                                                           <  2:1843996/100‑1 (MQ=255)
        gTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGtt                                                                                          <  2:2051389/100‑1 (MQ=255)
         tCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGTAGAAAACCCGCTGATGAGCGGGttt                                                                                         <  1:2560940/100‑1 (MQ=38)
          cACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGtttt                                                                                        <  2:1770606/100‑1 (MQ=255)
          cACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGtttt                                                                                        >  2:2172483/1‑100 (MQ=255)
            ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgt                                                                                      <  2:1659655/100‑3 (MQ=255)
            ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCt                                                                                      >  1:2591355/1‑100 (MQ=255)
             cAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTa                                                                                     >  2:798663/1‑100 (MQ=255)
              aGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgc                                                                                    <  2:344420/100‑5 (MQ=38)
              aGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgc                                                                                    <  2:2733262/100‑5 (MQ=39)
              aGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTAc                                                                                    >  1:2267910/1‑100 (MQ=255)
              aGCTAAATAATAAGTTAACAGACATGACCCCCGGTATGAGATTCCACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgc                                                                                    <  1:1898976/100‑5 (MQ=37)
               gCTAAATAATAAGTTAATAGACATTACATCCCGGATGAGATTCAACTTTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgc                                                                                    <  1:305065/99‑5 (MQ=25)
                cTAAATAATAAGTTAACAGACATGAGTCCCGTGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGTTTTGTGCtt                                                                                  <  1:1069567/100‑1 (MQ=38)
                 tAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACttt                                                                                 >  2:1577968/1‑100 (MQ=255)
                  aaaGAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCAGTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACtttt                                                                                <  1:1517646/100‑1 (MQ=255)
                   aataatAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACttttt                                                                               <  2:1484680/100‑1 (MQ=255)
                    ataataAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTc                                                                              >  2:1823837/1‑100 (MQ=255)
                    ataataAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTc                                                                              >  2:1725810/1‑100 (MQ=255)
                     taataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                             <  2:2099533/100‑1 (MQ=255)
                     taataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                             <  1:3139493/100‑1 (MQ=255)
                       ataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTa                                                                           >  2:695979/1‑100 (MQ=255)
                       ataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTa                                                                           <  2:2556554/100‑1 (MQ=255)
                       ataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGAGGGTTTTCGACTTTGTCt                                                                             <  1:286875/98‑1 (MQ=255)
                        taaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTaa                                                                          <  2:1217454/100‑1 (MQ=255)
                        taaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTaa                                                                          <  2:2591375/100‑1 (MQ=255)
                         aaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAAc                                                                         <  1:1390483/100‑1 (MQ=255)
                          aGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTa                                                                           >  1:338044/1‑97 (MQ=255)
                          aGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACg                                                                        >  1:2709797/1‑100 (MQ=255)
                            ttAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACACGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTc                                                                      >  1:3189095/1‑100 (MQ=255)
                             tAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCg                                                                     <  1:3339262/100‑1 (MQ=255)
                              aaCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCg                                                                     >  1:2205752/1‑99 (MQ=255)
                                cAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGTTTTTCTACTTTTTCTTAACGTCGGGt                                                                  <  2:2069775/100‑1 (MQ=255)
                                  gACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtat                                                                >  1:2665823/1‑100 (MQ=255)
                                  gACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtat                                                                <  2:1265188/100‑1 (MQ=255)
                                  gACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtat                                                                >  2:3104256/1‑100 (MQ=255)
                                   aCATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtata                                                               <  1:2496506/100‑1 (MQ=255)
                                    cATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATAc                                                              >  2:391956/1‑100 (MQ=255)
                                      tGAGTCCCGGGATGAGATTCAACATTACCATTGCCTCATTTAAAGCACAGAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGAtataa                                                              >  1:3346691/1‑97 (MQ=255)
                                      tGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATCCaa                                                            >  1:1954978/1‑100 (MQ=255)
                                      tGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACaa                                                            >  2:1445590/1‑100 (MQ=255)
                                         gTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACaaacc                                                         <  2:2278327/100‑3 (MQ=255)
                                         gTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTAATTTTTCTTAACGTCGGGTATACAAAGc                                                         >  2:3373803/1‑100 (MQ=255)
                                          tCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGcc                                                        <  1:3104230/100‑1 (MQ=255)
                                           cccGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGccc                                                       >  1:361827/1‑100 (MQ=255)
                                            ccGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCa                                                      <  1:3159517/100‑1 (MQ=255)
                                               ggATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCg                                                   <  2:2803213/100‑1 (MQ=255)
                                                  tGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTg                                                >  1:2919962/1‑100 (MQ=255)
                                                   gagaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTg                                                <  2:2770710/99‑1 (MQ=255)
                                                   gagaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGa                                               <  1:823275/100‑1 (MQ=255)
                                                   gagaTGCAACATTACCAGTGCCCCATTTAGAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGa                                               <  1:1662011/100‑1 (MQ=255)
                                                      aTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaaaa                                            >  1:1792229/1‑100 (MQ=255)
                                                          aaCATTACCATTGCCCCATTTAAAGCACACAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtt                                        <  1:2407482/100‑1 (MQ=255)
                                                          aaCATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtt                                        >  1:55158/1‑100 (MQ=255)
                                                          aaCATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtt                                        >  1:1617194/1‑100 (MQ=255)
                                                           aCATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAttt                                       >  1:821492/1‑100 (MQ=255)
                                                            cATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtttt                                      >  1:1098011/1‑100 (MQ=255)
                                                            cATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtttt                                      <  2:1850697/100‑1 (MQ=255)
                                                                      gCCCCATTTAAATCACAAAACCCGCTCATCAGCGGGTTTTCTACTTGTTCTTAACGTGTGTTATTCAAAGCCCATCTTTTAAAAAATTTTATCCatat                              >  1:2373864/1‑98 (MQ=255)
                                                                      gCCCCATTTAAAGCACACAACCCGCTCATCAGCGGGTTTTCTACTTTGTCTTAACGTCTGGTATACAAAGCCCATCGTTGAAAAAATTTTAt                                    <  1:3218838/92‑1 (MQ=255)
                                                                      gCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttt                            >  2:666367/1‑100 (MQ=255)
                                                                      gCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttt                            <  2:1959363/100‑1 (MQ=255)
                                                                       ccccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttt                           >  2:821808/1‑100 (MQ=255)
                                                                       ccccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttt                           >  1:445849/1‑100 (MQ=255)
                                                                        cccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttttt                          <  2:126293/100‑1 (MQ=255)
                                                                         ccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttttt                         >  1:1002316/1‑100 (MQ=255)
                                                                          cATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttttttn                        <  2:2858618/100‑2 (MQ=255)
                                                                          cATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTg                        <  2:1058445/100‑1 (MQ=255)
                                                                          cATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTg                        <  2:2098175/100‑1 (MQ=255)
                                                                             ttAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaa                     >  1:475679/1‑100 (MQ=255)
                                                                              aaaaGTAGAAAACCCGCTGATGAGCGGGTTTTGTGCTTTTTCTTAACGTCGGGTCTACAAAGCCCATCGTTGAAAGAAATTAATCCATATTTTGTGaaa                     >  1:3036870/2‑99 (MQ=18)
                                                                               aaaGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaaa                   <  2:2120773/100‑1 (MQ=255)
                                                                                aaGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAAt                  <  2:1002738/100‑1 (MQ=255)
                                                                                 agtagaaaaCCCGCTGATGAGCGGGTTTTGTGCATTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAGAAAATTTTATCCATATTTTTTGAAAAAt                  >  1:2038915/6‑99 (MQ=25)
                                                                                 agtagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaaa                   >  2:2772221/6‑98 (MQ=38)
                                                                                 aGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATg                 <  1:1577961/100‑1 (MQ=255)
                                                                                  gtagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                >  2:2920045/5‑100 (MQ=39)
                                                                                  gCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                >  2:1294526/1‑100 (MQ=255)
                                                                                  gCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                <  1:2448437/100‑1 (MQ=255)
                                                                                   tagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCa               >  2:1859448/4‑100 (MQ=39)
                                                                                   cacaAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCa               >  1:225258/1‑100 (MQ=255)
                                                                                    acaAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                <  1:1066263/98‑1 (MQ=255)
                                                                                    acaAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaa              <  2:2182014/100‑1 (MQ=255)
                                                                                    acaAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTATTTGAAAAATGCaa              <  1:513556/100‑1 (MQ=255)
                                                                                       aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAca           >  1:2003303/1‑100 (MQ=255)
                                                                                       aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAca           >  1:1516094/1‑100 (MQ=255)
                                                                                       aaaCCCGCTCATCAGCGGGTTTTCTACTATTTCTTAACGTCGGGTATACAAAGCGCATCGTTGAAAAAATAAAATCAATATTTTTTGAAAAATGCAAAca           >  1:656151/1‑100 (MQ=255)
                                                                                        aaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcat          >  2:1268993/1‑100 (MQ=255)
                                                                                        aaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGCAAAAAGTTTATCCAAATTTTTTGAAAAATGCAAAcat          >  1:2287630/1‑100 (MQ=255)
                                                                                          cccGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcatca        >  2:1598427/1‑100 (MQ=255)
                                                                                           ccGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcatcat       >  1:1695744/1‑100 (MQ=255)
                                                                                            cGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATg      >  1:422010/1‑100 (MQ=255)
                                                                                            cGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATg      >  2:1899336/1‑100 (MQ=255)
                                                                                             gCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGt     <  1:3186544/100‑1 (MQ=255)
                                                                                             gCTCATCAGCGGGATTTCTACTTTTTCTTAACGTCGGGTATACAAAACCCAGCGTTGAAAAAATTTTATCCATATTTTTTGTAAAATGCAGACCTCATGt     >  1:2521828/1‑100 (MQ=255)
                                                                                                catcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTccc  >  1:2841336/1‑98 (MQ=255)
                                                                                                   |                                                                                                
ATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGC  >  NC_000913/1642274‑1642469

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