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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,642,3700CG17.0% 320.5 / 57.1 118intergenic (+303/+119)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 131 total aligned reads displayed.

CAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGC  >  NC_000913/1642271‑1642469
                                                                                                   |                                                                                                   
cAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtc                                                                                                     >  2:2257459/1‑100 (MQ=255)
cAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtc                                                                                                     >  2:2901302/1‑100 (MQ=255)
 aaaTAACCATGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtca                                                                                                    <  2:622351/100‑1 (MQ=255)
 aaaTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtca                                                                                                    >  1:3115786/1‑100 (MQ=255)
 aaaTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtca                                                                                                    <  2:1614345/100‑1 (MQ=255)
 aaaTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtca                                                                                                    <  1:530798/100‑1 (MQ=255)
   aTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatc                                                                                                  <  1:27077/100‑1 (MQ=255)
   aTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatc                                                                                                  >  2:1094179/1‑100 (MQ=255)
   aTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatc                                                                                                  <  2:666822/100‑1 (MQ=255)
   aTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatc                                                                                                  >  2:2855424/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:1790607/100‑1 (MQ=255)
             cACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGtttt                                                                                        <  2:1770606/100‑1 (MQ=255)
               ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgt                                                                                      <  2:1597141/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)
                 aGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTAc                                                                                    >  1:2267910/1‑100 (MQ=255)
                 aGCTAAATAATAAGTTAACAGACATGACCCCCGGTATGAGATTCCACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgc                                                                                    <  1:1898976/100‑5 (MQ=37)
                  gCTAAATAATAAGTTAATAGACATTACATCCCGGATGAGATTCAACTTTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgc                                                                                    <  1:305065/99‑5 (MQ=25)
                  gCTAAATAATAAGTAAACAGACATGAGTCCCGGAATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgc                                                                                    <  1:1898124/99‑5 (MQ=38)
                   cTAAATAATAAGTTAACAGACATGAGTCCCGTGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGTTTTGTGCtt                                                                                  <  1:1069567/100‑1 (MQ=38)
                    tAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACttt                                                                                 >  1:1617118/1‑100 (MQ=255)
                     aaaGAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCAGTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACtttt                                                                                <  1:1517646/100‑1 (MQ=255)
                      aataatAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACttttt                                                                               <  2:1484680/100‑1 (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: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)
                                  aCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCggg                                                                   <  1:37003/100‑1 (MQ=255)
                                   cAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGTTTTTCTACTTTTTCTTAACGTCGGGt                                                                  <  2:2069775/100‑1 (MQ=255)
                                    aGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGta                                                                 >  1:875104/1‑100 (MQ=255)
                                     gACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtat                                                                <  2:1265188/100‑1 (MQ=255)
                                      aCATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtata                                                               <  1:2496506/100‑1 (MQ=255)
                                         tGAGTCCCGGGATGAGATTCAACATTACCATTGCCTCATTTAAAGCACAGAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGAtataa                                                              >  1:3346691/1‑97 (MQ=255)
                                         tGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACaa                                                            >  2:1445590/1‑100 (MQ=255)
                                            gTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTAATTTTTCTTAACGTCGGGTATACAAAGc                                                         >  2:3373803/1‑100 (MQ=255)
                                               ccGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCa                                                      <  1:3159115/100‑1 (MQ=255)
                                                  ggATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCg                                                   <  2:2803213/100‑1 (MQ=255)
                                                      gagaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTg                                                <  2:2770710/99‑1 (MQ=255)
                                                      gagaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGa                                               <  1:823275/100‑1 (MQ=255)
                                                      gagaTGCAACATTACCAGTGCCCCATTTAGAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGa                                               <  1:1662011/100‑1 (MQ=255)
                                                       agaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaa                                              >  2:3317415/1‑100 (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)
                                                                 ttACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTAt                                    <  2:2630528/100‑1 (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)
                                                                             cATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTg                        <  1:229736/100‑1 (MQ=255)
                                                                             cATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTg                        <  2:2098175/100‑1 (MQ=255)
                                                                             cATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTg                        <  2:1058445/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)
                                                                                    agtagaaaaCCCGCTGATGAGCGGGTTTTGTGCATTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAGAAAATTTTATCCATATTTTTTGAAAAAt                  >  1:2038915/6‑99 (MQ=25)
                                                                                    agtagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaaa                   >  2:2772221/6‑98 (MQ=38)
                                                                                     gtagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                >  2:2920045/5‑100 (MQ=39)
                                                                                     gCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCCTCGTTGAACAAATTTTATCCATATTTTTTGAAAAATGc                >  1:2992425/1‑100 (MQ=255)
                                                                                     gCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                >  2:1294526/1‑100 (MQ=255)
                                                                                     gCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                <  1:2448437/100‑1 (MQ=255)
                                                                                      tagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCa               >  2:1859448/4‑100 (MQ=39)
                                                                                      cacaAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCa               >  1:2114983/1‑100 (MQ=255)
                                                                                       acaAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                <  1:1066263/98‑1 (MQ=255)
                                                                                       acaAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaa              <  2:2182014/100‑1 (MQ=255)
                                                                                       acaAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaa              <  2:2181810/100‑1 (MQ=255)
                                                                                        gaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaaa             >  1:519006/2‑100 (MQ=255)
                                                                                          aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAca           >  1:2003303/1‑100 (MQ=255)
                                                                                          aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAca           >  1:1516094/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)
                                                                                               cGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATg      >  2:1899336/1‑100 (MQ=255)
                                                                                               cGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATg      >  1:422010/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)
                                                                                                   |                                                                                                   
CAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGC  >  NC_000913/1642271‑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: