Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 383,876 C→T 6.1% T20T (ACG→ACA)  yaiP ← putative family 2 glycosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913383,8760CT 6.1% 205.3 / 2.7 65T20T (ACG→ACA) yaiPputative family 2 glycosyltransferase
Reads supporting (aligned to +/- strand):  ref base C (25/36);  new base T (2/2);  total (27/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

GGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCAC  >  NC_000913/383732‑384024
                                                                                                                                                |                                                                                                                                                    
ggAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTAc                                                                                                                                                  >  2:315216/1‑149 (MQ=255)
   gCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTa                                                                                                                                               <  6:218639/149‑1 (MQ=255)
       tcgtTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTaaaaa                                                                                                                                           >  4:34050/1‑149 (MQ=255)
       tcgtTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTaaaaa                                                                                                                                           <  8:201152/149‑1 (MQ=255)
        cgtTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAAc                                                                                                                                          >  8:71962/1‑149 (MQ=255)
         gtTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAAcc                                                                                                                                         <  7:71962/149‑1 (MQ=255)
         gtTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAAc                                                                                                                                          >  8:88590/1‑148 (MQ=255)
            aaCGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCAt                                                                                                                                       <  7:88590/148‑1 (MQ=255)
             aCGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATa                                                                                                                                      <  4:328459/148‑1 (MQ=255)
             aCGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATa                                                                                                                                      <  3:151112/148‑1 (MQ=255)
             aCGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATa                                                                                                                                      <  1:342485/148‑1 (MQ=255)
             aCGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATa                                                                                                                                      <  8:342846/148‑1 (MQ=255)
                      aCCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGc                                                                                                                            >  4:282512/1‑149 (MQ=255)
                       ccaggCAAAAATCAGGGTTACGCAGTAGATTCTCCTGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATACTGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCt                                                                                                                           <  6:312397/146‑1 (MQ=255)
                        ccGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGTGTACTTAAAAACCATAGc                                                                                                                                    <  2:205138/139‑1 (MQ=255)
                        ccGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGTGTACTTAAAAACCATAGc                                                                                                                                    <  2:200049/139‑1 (MQ=255)
                        ccGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGTGTACTTAAAAACCATAGc                                                                                                                                    >  1:205138/1‑139 (MQ=255)
                        ccGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGTGTACTTAAAAACCATAGc                                                                                                                                    >  1:200049/1‑139 (MQ=255)
                          ggCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGTCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATg                                                                                                                        <  5:310775/149‑1 (MQ=255)
                          ggCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATg                                                                                                                        >  2:8451/1‑149 (MQ=255)
                           gCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATgg                                                                                                                       <  8:40361/149‑1 (MQ=255)
                             aaaaaTAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGAc                                                                                                                     <  3:242240/149‑1 (MQ=255)
                                aaTAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATa                                                                                                                  >  1:215781/1‑149 (MQ=255)
                                  tAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTAAGCGTACTTAAAAACCATAGCAATGTTGCTATggtttttaa                                                                                                                <  6:171444/149‑8 (MQ=255)
                                              aGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCa                                                                                                    <  1:18465/149‑1 (MQ=255)
                                                       tCCAGTGACTGCGCAAGACACGGGCCTTAGTTATACGCAGGTATAATGCCGTCAATACACCATTTCTTTTTACTCGGTTTAAGACTTACCGAAATTAAAAAACATAGCAACATTGCTATGGACATACAGATAAATATCAAGGTTATCAt                                                                                           >  8:321593/1‑149 (MQ=255)
                                                         cAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATg                                                                                          <  3:303054/148‑1 (MQ=255)
                                                         cAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATg                                                                                          >  6:292614/1‑148 (MQ=255)
                                                               cTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTAAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTg                                                                                   <  5:74007/149‑1 (MQ=255)
                                                               cTGCGCCAGACACGGGCCTGCGTTATACGCAGGTATAATGGCGTCAATCCAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTg                                                                                   <  2:338712/149‑1 (MQ=255)
                                                                  cgcCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCt                                                                                 <  7:160222/148‑1 (MQ=255)
                                                                     cAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCaaa                                                                             <  5:100517/149‑1 (MQ=255)
                                                                      aGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGc                                                                                         <  8:117905/136‑1 (MQ=255)
                                                                      aGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGc                                                                                         >  7:117905/1‑136 (MQ=255)
                                                                               ccTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAACACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAtttt                                                                   <  2:274754/149‑1 (MQ=255)
                                                                                 ttCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAttt                                                                    >  1:36514/1‑146 (MQ=255)
                                                                                 ttCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCAttt                                                                    <  2:36514/146‑1 (MQ=255)
                                                                                 ttCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTAt                                                                 >  5:300674/1‑149 (MQ=255)
                                                                                    gTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTATTGACTGTGTTTAAGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGaa                                                              <  1:104673/149‑1 (MQ=255)
                                                                                        taCGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTAc                                                          <  5:292614/149‑1 (MQ=255)
                                                                                                       gCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCAcc                                              >  5:254037/1‑146 (MQ=255)
                                                                                                          tCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATc                                        <  4:99900/149‑1 (MQ=255)
                                                                                                             aTACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGttt                                                                                               <  7:84774/91‑1 (MQ=255)
                                                                                                             aTACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGttt                                                                                               >  8:84774/1‑91 (MQ=255)
                                                                                                              tACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTa                                    <  4:78639/149‑1 (MQ=255)
                                                                                                                cAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTAt                                   <  6:254037/148‑1 (MQ=255)
                                                                                                                cAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCACGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTAtt                                  >  7:181989/1‑149 (MQ=255)
                                                                                                                 aGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTg                                 <  2:215781/149‑1 (MQ=255)
                                                                                                                 aGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTg                                 <  3:34050/149‑1 (MQ=255)
                                                                                                                    ctttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGAcc                              >  1:82995/1‑149 (MQ=255)
                                                                                                                     tttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCt                             >  5:219278/1‑149 (MQ=255)
                                                                                                                     tttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCt                             >  1:251100/1‑149 (MQ=255)
                                                                                                                     tttctttTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCt                             >  3:270374/1‑149 (MQ=255)
                                                                                                                         ttttGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCtgtt                          >  3:83455/1‑148 (MQ=255)
                                                                                                                           ttGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCtgtt                          <  4:83455/146‑1 (MQ=255)
                                                                                                                              aCTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATaaa                                                                                                          >  6:22185/1‑63 (MQ=255)
                                                                                                                              aCTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATaaa                                                                                                          <  5:22185/63‑1 (MQ=255)
                                                                                                                                  ggTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGaaa                <  3:155885/149‑1 (MQ=255)
                                                                                                                                  ggTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGa                  <  1:315216/147‑1 (MQ=255)
                                                                                                                                   gTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAAt               >  7:222135/1‑149 (MQ=255)
                                                                                                                                    tttACGTCTTAGCGTTCTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGaa                                                              <  6:199664/101‑1 (MQ=255)
                                                                                                                                    tttACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGaa                                                              >  5:199664/1‑101 (MQ=255)
                                                                                                                                    tttACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAAtt              <  1:8451/149‑1 (MQ=255)
                                                                                                                                           cTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGTTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAAc       >  2:118351/1‑149 (MQ=255)
                                                                                                                                           cTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAAc       >  6:20302/1‑149 (MQ=255)
                                                                                                                                                cGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCCCCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCAc  >  6:7981/1‑149 (MQ=255)
                                                                                                                                                cGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCACGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCAc  <  8:181989/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
GGAGCCGTCGTTAACGCAAATTACCCGGCAAAAATAAGGGTTACGCAGTAGATTATCCAGTGACTGCGCCAGACACGGGCCTTCGTTATACGCAGGTATAATGGCGTCAATACAGCCTTTCTTTTGACTGGGTTTACGTCTTAGCGTACTTAAAAACCATAGCAACATTGCTATGGACATACAGATAAATATCCAGGTTTTCATGCAATTTGCTCAAAGAATCATTTTATGAATTACAAAGCCTTCACCCAGATCGCTATTGACCTGTTGTCCGCGAAATTGTGCAACTGCAC  >  NC_000913/383732‑384024

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑