Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,642,373 C→G 16.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,3730CG16.0% 356.6 / 55.7 131intergenic (+306/+116)cspF/quuQQin prophage; cold shock protein/Qin prophage; putative antitermination protein Q
Reads supporting (aligned to +/- strand):  ref base C (66/43);  new base G (10/11);  total (77/54)
Fisher's exact test for biased strand distribution p-value = 3.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.85e-01

Only 100 of 136 total aligned reads displayed.

ATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCAT  >  NC_000913/1642274‑1642472
                                                                                                   |                                                                                                   
    ccAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCg                                                                                                 <  1:104585/100‑1 (MQ=255)
      aaGTTACCAGTTAAATAATAAGTTAACAGATATGAGTCCCGGGATGAGATTCAACATTACCATTGCTCTATTTAAAGCACAAAACCCGCTCATCAGCg                                                                                                 <  2:934265/98‑1 (MQ=255)
       aGTCACCAGCTAAATGATAAGTTAACAGACGTGAGTCCCGGGATGAGATTCAACATTACCATTGCGCCATTTAAAGCACAAAACCCGCTGATGAGCGGGt                                                                                              <  2:1409989/100‑1 (MQ=255)
       aGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGt                                                                                              <  2:116642/100‑1 (MQ=255)
         tCACCAGCTAAATAATCAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttt                                                                                            <  2:2129094/100‑1 (MQ=255)
          cACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGtttt                                                                                           >  1:1515937/1‑100 (MQ=255)
           aCCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTc                                                                                          <  2:543565/100‑1 (MQ=255)
           aCCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTc                                                                                          >  1:1254325/1‑100 (MQ=255)
           aCCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTc                                                                                          >  1:595027/1‑100 (MQ=255)
           aCCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTc                                                                                          <  1:996579/100‑1 (MQ=255)
            ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgt                                                                                         <  1:607641/100‑3 (MQ=255)
            ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCt                                                                                         >  2:592122/1‑100 (MQ=255)
              aGCTCAATATTGTGTTAACATGCATCAGCCCCGGGATGATATTCAACGTTACCATTGCCTCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtg                                                                                        <  2:3154316/99‑4 (MQ=19)
              aGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgc                                                                                       <  2:1582556/100‑5 (MQ=39)
              aGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTAc                                                                                       <  2:2412091/100‑1 (MQ=255)
              aGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTAc                                                                                       <  1:2334733/100‑1 (MQ=255)
              aGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAACGGGTTTTCTAc                                                                                       >  2:2582951/1‑100 (MQ=255)
               gCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACt                                                                                      >  1:3258127/1‑100 (MQ=255)
               gCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACt                                                                                      >  1:3225969/1‑100 (MQ=255)
               gCTAAATAATAAGTTAACAGACATGAGTCCCGGGAGGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgct                                                                                      <  2:3372650/100‑6 (MQ=38)
               gCTAAATAATAAGTAAACTGATATGAGTCACGGGAGAAGATTCAACATTACTATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgct                                                                                      <  2:1146306/100‑6 (MQ=25)
                cTAAATAATAAGTTAACAGATATGAGTCCCGGGATGAGAATCAACATTACAATTGCACCATTTAAAGCACGAAACTCGCTCATAAGCGGACTTTCTACtt                                                                                     >  2:2384371/1‑100 (MQ=255)
                cTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgctt                                                                                     <  1:767153/100‑7 (MQ=38)
                cTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgctt                                                                                     <  1:2641406/100‑7 (MQ=38)
                cTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACtt                                                                                     <  2:2986764/100‑1 (MQ=255)
                cTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACtt                                                                                     <  1:2816712/100‑1 (MQ=255)
                 tAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACttt                                                                                    >  2:180302/1‑100 (MQ=255)
                 tAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACttt                                                                                    >  2:3303398/1‑100 (MQ=255)
                    ataataAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAATCACAAAACCCGCTCATCAGCGGGTTTTCTACttttt                                                                                  >  2:1766496/1‑99 (MQ=255)
                     taataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                                <  2:3573760/100‑1 (MQ=255)
                     taataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                                <  1:3428140/100‑1 (MQ=255)
                     taataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                                >  1:2309868/1‑100 (MQ=255)
                     taataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                                <  1:1804769/100‑1 (MQ=255)
                       ataaTTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTa                                                                              <  2:2466049/100‑1 (MQ=255)
                       ataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTa                                                                              >  2:17066/1‑100 (MQ=255)
                       ataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTa                                                                              <  1:3576491/100‑1 (MQ=255)
                       ataaGTTAACAGACATGAGTCCCGAGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACATAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                                <  2:2372876/98‑1 (MQ=255)
                        taaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCACATTTAAAGCAAAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTaa                                                                             >  2:282215/1‑100 (MQ=255)
                          aGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACg                                                                           >  1:2519144/1‑100 (MQ=255)
                           gTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGt                                                                          >  2:2798534/1‑100 (MQ=255)
                              aaCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTCTCTACTTTTTCTTAACGTCg                                                                        >  1:888654/1‑99 (MQ=255)
                               aCAGACATGAGTCCCGTGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCggg                                                                      <  1:2960944/100‑1 (MQ=255)
                               aCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCggg                                                                      <  2:20924/100‑1 (MQ=255)
                               aCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCgg                                                                       >  2:546194/1‑99 (MQ=255)
                                   aCATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtata                                                                  <  1:1828867/100‑1 (MQ=255)
                                     aTGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACa                                                                <  2:1152637/100‑1 (MQ=255)
                                      tGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACaa                                                               <  2:2224429/100‑1 (MQ=255)
                                        aGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACaaa                                                              >  2:3416395/1‑99 (MQ=255)
                                          tCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGcc                                                           >  1:245556/1‑100 (MQ=255)
                                           cccGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGccc                                                          >  2:2088695/1‑100 (MQ=255)
                                                   gagaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGa                                                  <  1:650950/100‑1 (MQ=255)
                                                   gagaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGa                                                  <  2:2056316/100‑1 (MQ=255)
                                                     gaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaaa                                                <  2:371647/100‑1 (MQ=255)
                                                     gaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaaa                                                >  2:1736578/1‑100 (MQ=255)
                                                     gaTTCAACATTACCATTACCCAAATTAAAGCACAAAACCCGCTCATCAGCGGGTCTTCTACTTTTTTTTACAGTCGGGTATACAAACGCCAttgttg                                                   >  2:2310789/1‑97 (MQ=255)
                                                       ttCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaaaaa                                              >  1:3209313/1‑100 (MQ=255)
                                                         cAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAt                                            >  2:3612144/1‑100 (MQ=255)
                                                          aaCATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtt                                           >  2:291179/1‑100 (MQ=255)
                                                          aaCATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtt                                           >  2:746737/1‑100 (MQ=255)
                                                            cATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtttt                                         >  1:1108318/1‑100 (MQ=255)
                                                              ttACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTAt                                       <  1:3216014/100‑1 (MQ=255)
                                                                aCCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATcc                                     >  1:1921209/1‑100 (MQ=255)
                                                                 ccATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTCACGTCGGGTCTACAAAGCCAATCGTTGAAAAAATTATATc                                      >  1:2917317/1‑98 (MQ=255)
                                                                  cATTGCCCCATTTAAAGCACAAAACGCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCat                                   >  2:3606762/1‑100 (MQ=255)
                                                                  cATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCat                                   >  1:612807/1‑100 (MQ=255)
                                                                      gCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttt                               >  2:911114/1‑100 (MQ=255)
                                                                      gCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttt                               >  1:412571/1‑100 (MQ=255)
                                                                       ccccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttt                              >  2:723585/1‑100 (MQ=255)
                                                                        cccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttttt                             >  1:2098429/1‑100 (MQ=255)
                                                                        cccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttttt                             >  1:2403892/1‑100 (MQ=255)
                                                                          cATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTg                           <  2:73951/100‑1 (MQ=255)
                                                                          cATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTg                           >  1:458008/1‑100 (MQ=255)
                                                                           aTTTGAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttttt                            <  2:3568539/98‑1 (MQ=255)
                                                                              tAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaa                       >  2:161831/1‑100 (MQ=255)
                                                                              aaaagtagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaa                       >  1:480101/9‑100 (MQ=38)
                                                                                aaGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAAt                     >  2:1274411/1‑100 (MQ=255)
                                                                                 aGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATg                    >  2:67752/1‑100 (MQ=255)
                                                                                  gtagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                   >  1:279349/5‑100 (MQ=38)
                                                                                  gCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAAt                     <  1:1736586/98‑1 (MQ=255)
                                                                                   cacaAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCa                  <  1:3472379/100‑1 (MQ=255)
                                                                                    agaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaa                 >  2:3069980/3‑100 (MQ=255)
                                                                                     gaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaaa                >  1:1236846/2‑100 (MQ=255)
                                                                                     caAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATTGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaaa                >  2:3309663/1‑100 (MQ=255)
                                                                                     caAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAACAATTTTATCCATATTTTTTGAAAAATGCaaa                >  2:2036504/1‑100 (MQ=255)
                                                                                     caAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaaa                >  1:1300682/1‑100 (MQ=255)
                                                                                      aaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCTGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAc               >  2:740638/1‑100 (MQ=255)
                                                                                       aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAca              <  1:365114/100‑1 (MQ=255)
                                                                                        aaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCAGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcat             >  1:1959641/1‑100 (MQ=255)
                                                                                         aCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAGAAATTTTATCCATATTTTTTGAAAAATGCAAAcatc            >  1:2558064/1‑100 (MQ=255)
                                                                                          cccGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcatca           >  1:1774368/1‑100 (MQ=255)
                                                                                          cccGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcatca           <  1:17066/100‑1 (MQ=255)
                                                                                          cccGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcatca           >  1:206403/1‑100 (MQ=255)
                                                                                          cccGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCAGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcatca           <  2:1710770/100‑1 (MQ=255)
                                                                                            cGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATg         >  1:1037506/1‑100 (MQ=255)
                                                                                            cGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATg         >  1:901206/1‑100 (MQ=255)
                                                                                              cTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTc       <  1:2909137/100‑1 (MQ=255)
                                                                                                 atcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGcc    <  1:924367/100‑1 (MQ=255)
                                                                                                  tcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCa   <  1:3567673/100‑1 (MQ=255)
                                                                                                  tcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCa   >  2:2707762/1‑100 (MQ=255)
                                                                                                   caGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCAt  >  1:2010420/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
ATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCAT  >  NC_000913/1642274‑1642472

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: