Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,642,3700CG22.6% 317.1 / 88.3 138intergenic (+303/+119)cspF/quuQQin prophage; cold shock protein/Qin prophage; putative antitermination protein Q
Reads supporting (aligned to +/- strand):  ref base C (57/49);  new base G (15/16);  total (72/66)
Fisher's exact test for biased strand distribution p-value = 6.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 149 total aligned reads displayed.

CAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGC  >  NC_000913/1642271‑1642469
                                                                                                   |                                                                                                   
cAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtc                                                                                                     >  1:1396855/1‑100 (MQ=255)
 aaaTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtca                                                                                                    <  2:2596458/100‑1 (MQ=255)
  aaTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcat                                                                                                   >  2:1455338/1‑100 (MQ=255)
   aTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatc                                                                                                  <  1:2761387/100‑1 (MQ=255)
   aTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatc                                                                                                  >  2:1963262/1‑100 (MQ=255)
      aCCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGc                                                                                               >  1:1112969/1‑100 (MQ=255)
      aCCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGc                                                                                               >  2:926782/1‑100 (MQ=255)
       ccAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCg                                                                                              <  1:2130091/100‑1 (MQ=255)
       ccAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCg                                                                                              >  1:1820168/1‑100 (MQ=255)
        cAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGAATCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCg                                                                                              >  2:2714922/1‑99 (MQ=255)
         aaGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCggg                                                                                            <  1:609247/100‑1 (MQ=255)
          aGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGt                                                                                           <  1:546759/100‑1 (MQ=255)
           gTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGtt                                                                                          >  1:1430430/1‑100 (MQ=255)
            tCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttt                                                                                         <  1:2422752/100‑1 (MQ=255)
            tCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttt                                                                                         <  1:1876791/100‑1 (MQ=255)
             cACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGTATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGtttt                                                                                        >  2:1780422/1‑100 (MQ=255)
             cACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGtttt                                                                                        >  2:2268910/1‑100 (MQ=255)
             cACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGtttt                                                                                        <  1:2429277/100‑1 (MQ=255)
               ccAGCTAAATAATAAGTTACCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgt                                                                                      <  1:2595143/100‑3 (MQ=39)
               ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgt                                                                                      <  1:2351270/100‑3 (MQ=255)
               ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgt                                                                                      <  2:410627/100‑3 (MQ=255)
               ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCt                                                                                      >  1:1413975/1‑100 (MQ=255)
               ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCt                                                                                      <  2:481689/100‑1 (MQ=255)
                cAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtg                                                                                     <  1:1378131/100‑4 (MQ=39)
                cAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtg                                                                                     <  1:783150/100‑4 (MQ=39)
                cAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTa                                                                                     <  1:10915/100‑1 (MQ=255)
                 aGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgc                                                                                    <  2:2283766/100‑5 (MQ=39)
                 aGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTAc                                                                                    <  2:2312589/100‑1 (MQ=255)
                  gCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgct                                                                                   <  2:464557/100‑6 (MQ=38)
                   cTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgctn                                                                                  <  2:2471742/100‑7 (MQ=38)
                    tAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACttt                                                                                 >  2:2472433/1‑100 (MQ=255)
                     naataatAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACtttt                                                                                >  2:80044/2‑100 (MQ=255)
                     aaaTAATAAGTTAACAGACATGAGTCCCGGGGTGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGATTTTCTACtttt                                                                                <  1:1896517/100‑1 (MQ=255)
                     aaaTAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACtttt                                                                                >  2:1931640/1‑100 (MQ=255)
                     aaaTAATAAGTTAACAGACATGAGTCCCGGGACGAGATTCAACATTACCAGTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGTTTTGTGCttt                                                                                 <  2:824837/99‑1 (MQ=37)
                       ataataAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTc                                                                              <  1:2258541/100‑1 (MQ=255)
                        taataaGTTAACAGACATGAGTCCCGGGATGAGATTCAATATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                             <  1:866578/100‑1 (MQ=255)
                        taataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                             <  2:1404287/100‑1 (MQ=255)
                        taataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                             >  2:1938135/1‑100 (MQ=255)
                        taataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                             <  2:1324049/100‑1 (MQ=255)
                         aataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCtt                                                                            >  2:741672/1‑100 (MQ=255)
                            aaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAAc                                                                         <  1:551847/100‑1 (MQ=255)
                             aGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGTACGAAACCCTCTCATCAGCGGGTTTTCTACTTTTTCTTAACg                                                                        >  2:2027679/1‑100 (MQ=255)
                              gTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGt                                                                       >  1:1650936/1‑100 (MQ=255)
                                  aCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCggg                                                                   <  1:987128/100‑1 (MQ=255)
                                  aCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCggg                                                                   <  1:112843/100‑1 (MQ=255)
                                  aCAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTATTTCTTAACGTCggg                                                                   >  1:2842808/1‑100 (MQ=255)
                                   cAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGt                                                                  >  1:1995614/1‑100 (MQ=255)
                                     gACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtat                                                                >  1:787258/1‑100 (MQ=255)
                                     gACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtat                                                                >  1:2628622/1‑100 (MQ=255)
                                     gACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtat                                                                >  2:1531219/1‑100 (MQ=255)
                                         tGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCTGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACaa                                                            <  2:1552644/100‑1 (MQ=255)
                                         tGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACaa                                                            >  2:1005122/1‑100 (MQ=255)
                                            gTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAGGc                                                         <  2:1820155/100‑1 (MQ=255)
                                            gTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGc                                                         >  1:2166398/1‑100 (MQ=255)
                                              cccGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGccc                                                       <  1:1317855/100‑1 (MQ=255)
                                              cccGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGccc                                                       <  2:1033189/100‑1 (MQ=255)
                                                 gggATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATc                                                    <  1:1955159/100‑1 (MQ=255)
                                                  ggATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCg                                                   <  2:2166385/100‑1 (MQ=255)
                                                   gATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGt                                                  >  2:2528445/1‑100 (MQ=255)
                                                   gATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGt                                                  <  2:605108/100‑1 (MQ=255)
                                                     tGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTg                                                <  1:965859/100‑1 (MQ=255)
                                                            cAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTCTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAt                                         <  2:165490/100‑1 (MQ=255)
                                                             aaCATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtt                                        <  1:1317264/100‑1 (MQ=255)
                                                              aCATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAttt                                       <  2:2809936/100‑1 (MQ=255)
                                                               cATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtttt                                      >  1:1479616/1‑100 (MQ=255)
                                                               cATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAtttt                                      <  1:1024233/100‑1 (MQ=255)
                                                                    ccATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCa                                 >  2:1967222/1‑100 (MQ=255)
                                                                      aTTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCata                               <  1:926784/100‑1 (MQ=255)
                                                                         gCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttt                            >  1:2702586/1‑100 (MQ=255)
                                                                           cccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttttt                          >  2:2861991/1‑100 (MQ=255)
                                                                            ccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttttt                         >  2:1164015/1‑100 (MQ=255)
                                                                            ccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttttt                         <  2:329000/100‑1 (MQ=255)
                                                                              aTTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGa                       >  2:378818/1‑100 (MQ=255)
                                                                               tccgatctagaaaaCCCGCTGATGATCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaa                      <  1:112854/90‑1 (MQ=255)
                                                                                ttAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaa                     <  2:2792965/100‑1 (MQ=255)
                                                                                aaaaagtagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaa                     >  2:1366153/10‑100 (MQ=37)
                                                                                  aaagtagaaaaCCCGCTGATGAGCGGGTTTTGTGCTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaaa                   >  1:1574664/8‑100 (MQ=25)
                                                                                  aaaGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGCAAAAATTTTATCCATATTTTTTGaaaaa                   >  2:2744051/1‑100 (MQ=255)
                                                                                  aaaGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaaa                   >  2:547809/1‑100 (MQ=255)
                                                                                  aaaGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaaa                   >  1:724435/1‑100 (MQ=255)
                                                                                  aaaGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaaa                   >  1:1599026/1‑100 (MQ=255)
                                                                                    agtagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATg                 >  1:1490901/6‑100 (MQ=38)
                                                                                     gtagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                >  1:634395/5‑100 (MQ=39)
                                                                                     gtagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                >  1:1400656/5‑100 (MQ=39)
                                                                                     gCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                >  1:598005/1‑100 (MQ=255)
                                                                                      tagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCa               >  2:945377/4‑100 (MQ=39)
                                                                                       agaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAGAtcg                >  2:112854/3‑96 (MQ=38)
                                                                                       agaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaa              >  2:414733/3‑100 (MQ=255)
                                                                                       acaAAACCCGCTTATCAGCGGGTTTTATACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                <  2:548390/98‑1 (MQ=255)
                                                                                        taaaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaaa             <  1:2832958/99‑1 (MQ=255)
                                                                                        caAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaaa             <  2:2842787/100‑1 (MQ=255)
                                                                                        caAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaaa             >  2:1175979/1‑100 (MQ=255)
                                                                                          aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAca           >  1:136381/1‑100 (MQ=255)
                                                                                          aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAca           <  1:1940971/100‑1 (MQ=255)
                                                                                            aCCCGCTGATGAGCGGGTTTTGTGCTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcatc         >  1:2289039/1‑100 (MQ=38)
                                                                                                  tcatcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCg   <  1:2495229/100‑1 (MQ=255)
                                                                                                  tcatcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCg   <  1:2027693/100‑1 (MQ=255)
                                                                                                   catcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGc  >  2:1639772/1‑100 (MQ=255)
                                                                                                   catcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAGTTTATCTATATTTTTTGAAAAATGCAAACATCATGTCGc  >  2:2031900/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
CAAATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGC  >  NC_000913/1642271‑1642469

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