Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,344,172 T→C S59G (AGT→GGT)  yciT ← putative DNA‑binding transcriptional regulator YciT

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,344,1720TC98.5% 228.4 / NA 65S59G (AGT→GGT) yciTputative DNA‑binding transcriptional regulator YciT
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (31/33);  minor base G (1/0);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTGACCC  >  NC_000913/1344037‑1344302
                                                                                                                                       |                                                                                                                                  
tCCGGGCGAGTAGGGCATTGCTTCTGCCATTGTCGCTAACGATGGTTTAGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCa                                                                                                                                <  2:325215/140‑1 (MQ=255)
     gtgAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGa                                                                                                                           <  2:616113/138‑1 (MQ=255)
     gCGAGTAGGGCATTGTTGCTGCAATTTTCGATAAAGAAGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGa                                                                                                                           <  1:366022/140‑1 (MQ=255)
       gAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGaaa                                                                                                                         <  2:117883/140‑1 (MQ=255)
        aGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAAc                                                                                                                        >  1:717848/1‑140 (MQ=255)
        aGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAAc                                                                                                                        >  1:92517/1‑140 (MQ=255)
        aGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAAc                                                                                                                        >  2:40140/1‑140 (MQ=255)
          tAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAg                                                                                                                      >  1:1005718/1‑140 (MQ=255)
            gggCATTGCTGCTGCCATTTTCGATAAAGATGGTTACGCACGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCa                                                                                                                    <  2:382105/140‑1 (MQ=255)
             ggCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCGTACGGGTGTCGACGTCGTCACGATCAAGCGAAACAGCaa                                                                                                                   >  2:295990/1‑140 (MQ=255)
              gCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCaaa                                                                                                                  <  2:1005718/140‑1 (MQ=255)
                 ttGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGcc                                                                                                               >  1:550119/1‑140 (MQ=255)
                  tgctgctgcCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCa                                                                                                              >  1:209670/1‑140 (MQ=255)
                    ctgctgcCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCAtg                                                                                                            >  1:511995/1‑140 (MQ=255)
                         gcCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCAc                                                                                                       >  1:539417/1‑140 (MQ=255)
                         gcCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCAc                                                                                                       <  2:198972/140‑1 (MQ=255)
                           cATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCCTTAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCAcgg                                                                                                     <  2:685682/140‑1 (MQ=255)
                              tttCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCAcggcgg                                                                                                  >  2:675664/1‑140 (MQ=255)
                                  gATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATTAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGt                                                                                              >  2:435084/1‑140 (MQ=255)
                                       aGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTc                                                                                         <  2:245648/140‑1 (MQ=255)
                                       aGATGGTTTCGCCCGGATGAACCAGCGACGTGGAATACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTc                                                                                         <  2:484512/140‑1 (MQ=255)
                                             tttCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGtttt                                                                                   >  2:827192/1‑140 (MQ=255)
                                                cGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCg                                                                                <  1:123893/140‑1 (MQ=255)
                                                cGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCg                                                                                <  1:827192/140‑1 (MQ=255)
                                                cGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCg                                                                                <  1:41063/140‑1 (MQ=255)
                                                cGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCg                                                                                >  1:555662/1‑140 (MQ=255)
                                                 gCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGa                                                                               <  1:850681/140‑1 (MQ=255)
                                                 gCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGa                                                                               <  2:550119/140‑1 (MQ=255)
                                                 gCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGa                                                                               <  2:209670/140‑1 (MQ=255)
                                                 gCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGa                                                                               <  1:897053/140‑1 (MQ=255)
                                                 gCCCGGTTGAACCAGCGACGCGCCACACTTGGCGAGTCCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGa                                                                               <  1:295990/140‑1 (MQ=255)
                                                  cccGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAg                                                                              >  2:614539/1‑140 (MQ=255)
                                                  cccGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAg                                                                              >  2:97297/1‑140 (MQ=255)
                                                  cccGGTAGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGAATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAg                                                                              >  2:616733/1‑140 (MQ=255)
                                                   ccGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAgt                                                                             <  2:92517/140‑1 (MQ=255)
                                                   ccGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAgt                                                                             <  1:40140/140‑1 (MQ=255)
                                                                cGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCt                                                                >  1:232794/1‑140 (MQ=255)
                                                                     cgGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGa                                                           <  2:675187/140‑1 (MQ=255)
                                                                             tCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAAc                                                   >  2:683590/1‑140 (MQ=255)
                                                                              cGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGTGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACt                                                  <  2:511995/140‑1 (MQ=255)
                                                                                              cAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTgg                                  >  2:226520/1‑140 (MQ=255)
                                                                                               aGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGc                                 >  2:611159/1‑140 (MQ=255)
                                                                                                    ataGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTg                            >  1:870872/1‑140 (MQ=255)
                                                                                                     taGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGc                           <  2:232794/140‑1 (MQ=255)
                                                                                                      aGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGcc                          >  1:969650/1‑140 (MQ=255)
                                                                                                      aGTTGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGcc                          >  1:979056/1‑140 (MQ=255)
                                                                                                        ttGCTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAg                        <  1:614539/140‑1 (MQ=255)
                                                                                                          gCTCATCATACGGGTATCGACGTCAGCACCATCAAGCCAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGAt                      <  1:435084/140‑1 (MQ=255)
                                                                                                           cTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATc                     >  2:914656/1‑140 (MQ=255)
                                                                                                           cTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATc                     >  1:167739/1‑140 (MQ=255)
                                                                                                           cTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATc                     >  1:593562/1‑140 (MQ=255)
                                                                                                           cTCATCATACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATc                     >  2:165022/1‑140 (MQ=255)
                                                                                                             catcatACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGATACTCCAGTGACTTTTGCCAGATCgg                   >  1:483306/1‑140 (MQ=255)
                                                                                                             catcatACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCgg                   >  2:663637/1‑140 (MQ=255)
                                                                                                              atcatACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGt                  <  1:675664/140‑1 (MQ=255)
                                                                                                                  tACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACg              <  2:979056/140‑1 (MQ=255)
                                                                                                                  tACGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACg              <  1:97297/140‑1 (MQ=255)
                                                                                                                    cGGGTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCt            <  1:683590/140‑1 (MQ=255)
                                                                                                                       gTTTCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAAc         >  1:654445/1‑140 (MQ=255)
                                                                                                                        tttCGACGTCATCACCAGCAAGGGAAACAGCAAAGGCATCTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAAcc        <  2:339091/140‑1 (MQ=255)
                                                                                                                          tCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTATCTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTg      <  1:165022/140‑1 (MQ=255)
                                                                                                                          tCGACGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGAAACTTCAGAAACGCCAGCGGCTTTTACCAGGTCGGTTACCCTAACCTg      >  2:478600/1‑140 (MQ=255)
                                                                                                                              cGTCATCACCATCAAGCGAAACTGCCAAGCCATGTGCCCGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTGAccc  <  1:314120/140‑1 (MQ=255)
                                                                                                                              cGTCATCACCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTGAccc  <  2:167739/140‑1 (MQ=255)
                                                                                                                                 catcaCCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACg              >  1:788783/1‑125 (MQ=255)
                                                                                                                                 catcaCCATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACg              <  2:788783/125‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                  
TCCGGGCGAGTAGGGCATTGCTGCTGCCATTTTCGATAAAGATGGTTTCGCCCGGTTGAACCAGCGACGCGGCAAACTCGGCGAGTTCACGCTTCAGCGTATAGTTGCTCATCATACGGGTTTCGACGTCATCACTATCAAGCGAAACAGCAAAGCCATGTGCACGGCGGAGGTAACTCAGTTTTTCGAGTGTGTTGAGATCCTGGCGAATGGTAACTTCAGAAACTCCAGTGGCTTTTGCCAGATCGGTTACGCTAACCTGACCC  >  NC_000913/1344037‑1344302

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑