Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,642,373 C→G 17.8% 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.8% 343.2 / 65.8 136intergenic (+306/+116)cspF/quuQQin prophage; cold shock protein/Qin prophage; putative antitermination protein Q
Reads supporting (aligned to +/- strand):  ref base C (68/43);  new base G (12/12);  total (80/56)
Fisher's exact test for biased strand distribution p-value = 3.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.61e-01

Only 100 of 148 total aligned reads displayed.

ATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCA  >  NC_000913/1642274‑1642471
                                                                                                   |                                                                                                  
agaaCCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCAGTGAAAGCACAAAACCCGCtcatc                                                                                                    <  1:3806085/98‑1 (MQ=255)
aTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatc                                                                                                    <  1:15399/100‑1 (MQ=255)
 tAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatca                                                                                                   >  2:469653/1‑100 (MQ=255)
 tAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatca                                                                                                   <  2:796106/100‑1 (MQ=255)
   aCCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGc                                                                                                 >  1:2332156/1‑100 (MQ=255)
   aCCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGc                                                                                                 >  2:138259/1‑100 (MQ=255)
    ccAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCg                                                                                                <  1:1213054/100‑1 (MQ=255)
      aaGTCACCAGCTATATATTTAGTTAACAGACATGAGTCCCGGAATGAGATTCAACATTACCCTTGCCCCCTTTACAGCACACAACCCGCTCATAAGCggg                                                                                              <  1:687830/100‑1 (MQ=255)
      aaGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCggg                                                                                              <  2:2593543/100‑1 (MQ=255)
      aaGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCggg                                                                                              <  1:4106591/100‑1 (MQ=255)
       aGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGt                                                                                             >  2:3534900/1‑100 (MQ=255)
          cACCAGGTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttt                                                                                           <  1:315589/99‑1 (MQ=255)
           aCCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGttct                                                                                          >  1:1618177/1‑97 (MQ=255)
           aCCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTc                                                                                         >  1:1617518/1‑100 (MQ=255)
           aCCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTc                                                                                         >  2:1596774/1‑100 (MQ=255)
            ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgt                                                                                        <  1:1218475/100‑3 (MQ=255)
            ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttg                                                                                         <  1:4363334/99‑2 (MQ=255)
            ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCt                                                                                        >  2:3310022/1‑100 (MQ=255)
             cAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtg                                                                                       <  2:1075446/100‑4 (MQ=39)
             cAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtg                                                                                       <  2:4182467/100‑4 (MQ=255)
             cAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTa                                                                                       >  2:2495907/1‑100 (MQ=255)
             cAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTa                                                                                       >  2:907999/1‑100 (MQ=255)
               gcgaaaTAAGAAGTTAACAGACATGAGGCCCGGGATGAGAGTCAACATTACCATTGGCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtg                                                                                       <  1:2050883/95‑4 (MQ=37)
               gCTAAAAAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgc                                                                                      <  2:2007253/99‑5 (MQ=38)
                cTAAATAATACGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCACTGCCCCCTTTCAAGCACACAACCCGCTCATCAGCGGGCTTTCTACt                                                                                     >  1:2125792/1‑99 (MQ=255)
                   aataatAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACttttt                                                                                 >  1:3529618/1‑100 (MQ=255)
                    ataataAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTc                                                                                >  2:3622785/1‑100 (MQ=255)
                          aGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACg                                                                          >  2:2534660/1‑100 (MQ=255)
                          aGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACg                                                                          >  2:566982/1‑100 (MQ=255)
                          aGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACg                                                                          <  2:3529637/100‑1 (MQ=255)
                          aGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACg                                                                          <  2:977761/100‑1 (MQ=255)
                            ttAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCtttttttt                                                                        >  1:390321/1‑94 (MQ=255)
                            ttAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTc                                                                        >  2:70432/1‑100 (MQ=255)
                               aCAGACATGAGTTCCGGGATGAGAAACAATATTACCATTGCCACATTTAAAGCACAAAACCCACTCATCACCTGCATATCTAATTTTACTTAAc                                                                           >  1:3014945/1‑94 (MQ=255)
                               aCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCggg                                                                     <  2:3008200/100‑1 (MQ=255)
                                cAGACATGAGTCCCGGGATGAGATTCAACATGACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGt                                                                    <  2:2808141/100‑1 (MQ=255)
                                  gACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtat                                                                  >  1:453519/1‑100 (MQ=255)
                                     aTGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACa                                                               <  2:3113639/100‑1 (MQ=255)
                                     aTGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACa                                                               <  1:3281912/100‑1 (MQ=255)
                                      tGAGTGCCTAGATGAGATTCAACATTACCATTTCATCATTTAAAGCACAAAACCCGCTCATCAGCAGGTTTTCTACTTTTTCTTAGCCTCGGGTATACa                                                               <  1:4167162/99‑1 (MQ=255)
                                         gTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGc                                                           >  2:4053935/1‑100 (MQ=255)
                                         gTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGc                                                           >  1:4333693/1‑100 (MQ=255)
                                         gTCCCGGGAGGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGc                                                           <  1:3246256/100‑1 (MQ=255)
                                           cccGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGccc                                                         >  2:437337/1‑100 (MQ=255)
                                            ccGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCTTGTATACAAAACCCa                                                        >  1:808098/1‑100 (MQ=255)
                                               ggATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCg                                                     >  2:4150678/1‑100 (MQ=255)
                                                  tGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTg                                                  <  1:3175593/100‑1 (MQ=255)
                                                  tGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTg                                                  <  1:1020776/100‑1 (MQ=255)
                                                   gagaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGa                                                 <  1:136744/100‑1 (MQ=255)
                                                     gaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaaa                                               >  2:2862236/1‑100 (MQ=255)
                                                     gaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaaa                                               >  1:857413/1‑100 (MQ=255)
                                                      aTTCAACATTACCATTGCCCCATGTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaaaa                                              <  2:188633/100‑1 (MQ=255)
                                                       ttCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAACCCATTGTTGaaaaa                                             >  1:4179084/1‑100 (MQ=255)
                                                         cAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAt                                           >  1:3706507/1‑100 (MQ=255)
                                                         cAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAt                                           >  2:1098469/1‑100 (MQ=255)
                                                          aaCATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCTGGTATACAAAGCCCATCGTTGAAAAAAt                                           <  2:3321580/99‑1 (MQ=255)
                                                           aCATTACCATTGGCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAttt                                         >  2:2901384/1‑100 (MQ=255)
                                                           aCATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAttt                                         <  2:1743441/100‑1 (MQ=255)
                                                            cATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtttt                                        >  2:2701716/1‑100 (MQ=255)
                                                             aTTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTa                                       <  2:1037475/100‑1 (MQ=255)
                                                              ttACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACATCGGGTATACAAAGCCCATCGTTGAAAAAATTTTAt                                      <  2:3108672/100‑1 (MQ=255)
                                                                aCCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATcc                                    <  2:2116699/100‑1 (MQ=255)
                                                                  cattcacccATTTAAAGCACAAAACCCGCGCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGATGAAAAAATTTTATCCat                                  >  1:3545186/7‑100 (MQ=255)
                                                                  cATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCat                                  >  2:3726063/1‑100 (MQ=255)
                                                                   aTTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCata                                 >  2:3051581/1‑100 (MQ=255)
                                                                       ctccATTTAAAGCACAAAACCCGTTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGTGTATACTTAGCCCATCGTTGAAAAAATTTTATCCATAtttt                             <  1:2729808/98‑1 (MQ=255)
                                                                       ccccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTCTCCATAtttt                             >  1:4242558/1‑100 (MQ=255)
                                                                       ccccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttt                             <  1:3504641/100‑1 (MQ=255)
                                                                       ccccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttt                             >  2:2748551/1‑100 (MQ=255)
                                                                        cccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttttt                            >  1:3592196/1‑100 (MQ=255)
                                                                         ccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttttt                           >  2:843276/1‑100 (MQ=255)
                                                                         ccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttttt                           <  1:494243/100‑1 (MQ=255)
                                                                         ccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttttt                           <  2:3650318/100‑1 (MQ=255)
                                                                          cATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttttt                           <  2:1589198/99‑1 (MQ=255)
                                                                          cATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttttt                           <  1:2091490/99‑1 (MQ=255)
                                                                          cATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTg                          >  2:2919979/1‑100 (MQ=255)
                                                                          cATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTg                          <  2:2307105/100‑1 (MQ=255)
                                                                           aTTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGa                         >  1:3879359/1‑100 (MQ=255)
                                                                              tAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaa                      >  1:2736913/1‑100 (MQ=255)
                                                                                aaGAACAAAACCCGCTCATCAGCGAGTTCTCTACTTTTTCTTAACGTCGTGTATACATAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAAt                    <  1:1526023/100‑1 (MQ=255)
                                                                                 agtagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATg                   >  2:1010395/6‑100 (MQ=38)
                                                                                  gtagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGACAAAATTTTATCCATATTTTTTGAAAAATGc                  >  2:938087/5‑100 (MQ=38)
                                                                                  gtagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAAt                    >  1:1918353/5‑98 (MQ=255)
                                                                                  gCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                  >  2:1757864/1‑100 (MQ=255)
                                                                                   tagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAATTTTTGGTATACAAAGCCCATCGTTGAAAAAATTTTATTCATATTTTTTGAAAAATGc                  >  1:2911548/4‑99 (MQ=37)
                                                                                   cacaAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCa                 >  1:4130433/1‑100 (MQ=255)
                                                                                    agaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaa                >  1:2766079/3‑100 (MQ=255)
                                                                                    agaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaa                >  1:2451636/3‑100 (MQ=255)
                                                                                     gaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaaa               >  1:330209/2‑100 (MQ=255)
                                                                                     caAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaaa               >  1:606090/1‑100 (MQ=255)
                                                                                       acaCCCGCTGATGAGCGGGTTTTGTGCTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTAATACATATTTTTTGAAAAATCCAAAc              >  1:2293050/3‑99 (MQ=37)
                                                                                       aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAca             >  2:4046893/1‑100 (MQ=255)
                                                                                       aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAca             <  2:447028/100‑1 (MQ=255)
                                                                                            cGCTCATCAGCGGGTTTTCTTCTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATAATg        >  1:2024343/1‑100 (MQ=255)
                                                                                               tcatcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCg     <  1:522888/100‑1 (MQ=255)
                                                                                                catcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCATACATCATGTCGc    <  1:140966/100‑1 (MQ=255)
                                                                                                catcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGc    >  2:1254867/1‑100 (MQ=255)
                                                                                                 atcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGcc   <  1:3495657/100‑1 (MQ=255)
                                                                                                 atcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGcc   <  1:985436/100‑1 (MQ=255)
                                                                                                  tcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCa  <  2:3675917/100‑1 (MQ=255)
                                                                                                   |                                                                                                  
ATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCA  >  NC_000913/1642274‑1642471

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: