Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,642,373 C→G 20.9% 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,3730CG20.9% 227.7 / 52.9 95intergenic (+306/+116)cspF/quuQQin prophage; cold shock protein/Qin prophage; putative antitermination protein Q
Reads supporting (aligned to +/- strand):  ref base C (48/27);  new base G (7/13);  total (55/40)
Fisher's exact test for biased strand distribution p-value = 2.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

ATAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCAT  >  NC_000913/1642274‑1642472
                                                                                                   |                                                                                                   
aTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatc                                                                                                     <  2:571591/100‑1 (MQ=255)
aTAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatc                                                                                                     <  2:2207091/100‑1 (MQ=255)
 tAACCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatca                                                                                                    >  1:2232946/1‑100 (MQ=255)
  atccAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCtcatc                                                                                                     <  1:2795408/96‑1 (MQ=255)
   aCCAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGc                                                                                                  >  2:1511417/1‑100 (MQ=255)
    ccAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCg                                                                                                 <  1:1044676/100‑1 (MQ=255)
    ccAAGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCg                                                                                                 >  1:1494341/1‑100 (MQ=255)
      aaGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCggg                                                                                               >  2:2089204/1‑100 (MQ=255)
       aGTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGt                                                                                              <  2:805196/100‑1 (MQ=255)
        gTCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGtt                                                                                             <  2:1588263/100‑1 (MQ=255)
         tCACCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttt                                                                                            <  2:387254/100‑1 (MQ=255)
           aCCAGCTGAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTc                                                                                          <  1:918713/100‑1 (MQ=255)
           aCCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttg                                                                                          <  2:2533172/100‑2 (MQ=255)
           aCCAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATGCAACGTTACCATTGCCTCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttg                                                                                          <  1:69020/100‑2 (MQ=38)
            ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgt                                                                                         <  1:1110601/100‑3 (MQ=255)
            ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgt                                                                                         <  2:2502524/100‑3 (MQ=255)
            ccAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAAAATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGAGTTTTc                                                                                          >  1:2261462/1‑99 (MQ=255)
             cAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTa                                                                                        >  2:508222/1‑100 (MQ=255)
             cAGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTa                                                                                        >  1:147398/1‑100 (MQ=255)
             cAGCTAAAGAATAAGTTAACACACATCAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgt                                                                                         <  1:2147649/99‑3 (MQ=38)
              aGCTAAATAATAAGTTAACAGATATGAGTACCGGAATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgc                                                                                       <  1:1290233/100‑5 (MQ=37)
              aGCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgc                                                                                       <  2:2513307/100‑5 (MQ=39)
               gCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGttttgtgct                                                                                      <  2:2223543/100‑6 (MQ=38)
               gCTAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACt                                                                                      >  2:637528/1‑100 (MQ=255)
                 tAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACttt                                                                                    >  2:2092152/1‑100 (MQ=255)
                 tAAATAATAAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACttt                                                                                    >  1:71991/1‑100 (MQ=255)
                   aataatAAGTTAACAGACATGAGTCCCGGGATGATATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGTTTTGTGCtttaa                                                                                  <  2:2063351/100‑3 (MQ=37)
                    ataataAGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTc                                                                                 <  2:1210085/100‑1 (MQ=255)
                    ataataAATTATCAGACTGGAGTCCCGGAAGGAAATTCAACACTACTCTTGCCCCATTTAAAGCACAAAACCCGCTGATGAGCGGGTTTTGTGCtttaa                                                                                  <  1:1267088/99‑3 (MQ=12)
                     taataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCTTTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                                <  2:1678324/100‑1 (MQ=255)
                     taataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                                <  1:1684148/100‑1 (MQ=255)
                     taataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCt                                                                                <  1:1950130/100‑1 (MQ=255)
                       ataaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTa                                                                              <  1:2775787/100‑1 (MQ=255)
                         aaGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAAc                                                                            <  1:1467074/100‑1 (MQ=255)
                          aGTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAAc                                                                            <  2:1790123/99‑1 (MQ=255)
                           gTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGt                                                                          >  2:2477859/1‑100 (MQ=255)
                           gTTAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGt                                                                          >  1:1341581/1‑100 (MQ=255)
                            ttAACAGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTc                                                                         >  1:1665099/1‑100 (MQ=255)
                                 aGACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGta                                                                    >  2:2167338/1‑100 (MQ=255)
                                  gACATGAGTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGtat                                                                   >  2:70798/1‑100 (MQ=255)
                                    cATGAGTCCTGGGATGAGATTCATCATTATCCTTGCCCCATTTAATGCATACAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGta                                                                    <  1:2264739/97‑1 (MQ=255)
                                         gTCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGc                                                            >  2:383096/1‑100 (MQ=255)
                                          tCCCGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGcc                                                           <  1:1972134/100‑1 (MQ=255)
                                            ccGGGATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCa                                                         <  2:1304410/100‑1 (MQ=255)
                                              gggATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATc                                                       <  2:2020205/100‑1 (MQ=255)
                                               ggATGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCg                                                      >  2:2569085/1‑100 (MQ=255)
                                                 aTGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGtt                                                    >  1:122400/1‑100 (MQ=255)
                                                 aTGAGATTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGtt                                                    >  1:1250626/1‑100 (MQ=255)
                                                 aTGAGATTAAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGtt                                                    >  2:1341056/1‑100 (MQ=255)
                                                    agaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaa                                                 >  1:1660383/1‑100 (MQ=255)
                                                    agaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaa                                                 >  2:609124/1‑100 (MQ=255)
                                                     gaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaaa                                                >  2:2784177/1‑100 (MQ=255)
                                                     gaTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGATATACAAAGCCCATCGTTGaaa                                                >  1:1420379/1‑100 (MQ=255)
                                                      aTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaaaa                                               >  1:1590404/1‑100 (MQ=255)
                                                      aTTCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaaaa                                               >  2:643063/1‑100 (MQ=255)
                                                       ttCAACATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGaaaaa                                              >  1:444164/1‑100 (MQ=255)
                                                           aCATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAttt                                          >  1:2419804/1‑100 (MQ=255)
                                                           aCATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAAttt                                          <  1:2322394/100‑1 (MQ=255)
                                                            cATTACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTCACGTCAGGTATACAAAGCCCATCGTTGAAAAAAtttt                                         >  1:1956279/1‑100 (MQ=255)
                                                              ttACCATTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTAt                                       <  2:1914506/100‑1 (MQ=255)
                                                               tACCCTTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTAt                                       <  1:2586798/99‑1 (MQ=255)
                                                                   aTTGCCCCATTTAAAGCACAAAACCCGCTCATCGGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCata                                  <  2:2399125/100‑1 (MQ=255)
                                                                   aTTGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCata                                  >  2:620296/1‑100 (MQ=255)
                                                                     tGCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGCTCTCTACTTTTTCTTACCGTCGGGTATACAAAGCCCATAGTTGAAACAATTTTATCCatat                                 >  1:672147/1‑99 (MQ=255)
                                                                      gCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttt                               >  1:969713/1‑100 (MQ=255)
                                                                      gCCCCATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttt                               >  1:2714234/1‑100 (MQ=255)
                                                                        cccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttttt                             >  2:1042038/1‑100 (MQ=255)
                                                                        cccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAttttt                             >  2:1402674/1‑100 (MQ=255)
                                                                         ccATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATAtttttt                            >  2:1537897/1‑100 (MQ=255)
                                                                          cATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTg                           <  1:146001/100‑1 (MQ=255)
                                                                          cATTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTg                           <  2:2384627/100‑1 (MQ=255)
                                                                           aTTTAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGa                          <  2:2714244/100‑1 (MQ=255)
                                                                              tAAAGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaa                       >  2:1087350/1‑100 (MQ=255)
                                                                               aaagtagaaaaCCCGCTGATGAGCGGGTTTTGTGCTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaaa                      >  2:508484/8‑100 (MQ=25)
                                                                               aaaGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGaaaaa                      >  1:1115335/1‑100 (MQ=255)
                                                                                aaGCACAAATCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAAt                     <  1:1226698/100‑1 (MQ=255)
                                                                                aaGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAAt                     >  2:1949990/1‑100 (MQ=255)
                                                                                 tgCACAAAACCCGCTCATCAGCGGGTTCTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATg                    >  1:718095/2‑100 (MQ=255)
                                                                                 aGCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATg                    >  2:63818/1‑100 (MQ=255)
                                                                                  gtagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                   >  2:505594/5‑100 (MQ=39)
                                                                                  gtagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGc                   >  2:691514/5‑100 (MQ=255)
                                                                                  gCACAAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATg                    >  1:2541392/1‑99 (MQ=255)
                                                                                   tagaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCa                  >  2:1025664/4‑100 (MQ=39)
                                                                                    agaaaaCCCGCTGATGAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATAAAAAGTCCATCGTTGAAAAAATTTCATGCATATTGTTTGAAAAATGc                   >  1:591877/3‑98 (MQ=37)
                                                                                    agaaaaCCCGCTCAAGAGCGGGTGTTCTACTTTTTCTTAACGTCGGGTATAGAACGACCCACGTTGAAAAACTTTTATCTATATTTTTTGAAAAATGCaa                 >  1:935746/3‑100 (MQ=255)
                                                                                    acaAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCa                  <  1:611021/99‑1 (MQ=255)
                                                                                     caAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaaa                >  2:2477517/1‑100 (MQ=255)
                                                                                     caAAACCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCaaa                >  2:25509/1‑100 (MQ=255)
                                                                                       aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGTGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAc               <  1:1778388/99‑1 (MQ=255)
                                                                                       aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAca              >  1:336830/1‑100 (MQ=255)
                                                                                       aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAca              >  2:2684124/1‑100 (MQ=255)
                                                                                       aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAca              >  2:1687227/1‑100 (MQ=255)
                                                                                       aaaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAca              <  2:2232956/100‑1 (MQ=255)
                                                                                        aaCCCGCTCATCAGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAAcat             <  2:2906677/100‑1 (MQ=255)
                                                                                               tcatcaGCGGGTTTTCTACTTTTTCTTATCGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCg      <  2:2679664/100‑1 (MQ=255)
                                                                                               tcatcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCg      >  2:2252329/1‑100 (MQ=255)
                                                                                                catcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGc     <  2:1034571/100‑1 (MQ=255)
                                                                                                  tcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCa   <  1:2393105/100‑1 (MQ=255)
                                                                                                  tcaGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCa   <  1:2748274/100‑1 (MQ=255)
                                                                                                   caGCGGGTTTTCTACTTTTTCTTAACGTCGGGTATACAAAGCCCATCGTTGAAAAAATTTTATCCATATTTTTTGAAAAATGCAAACATCATGTCGCCAt  >  2:1355932/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: