Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,407,294 T→A 11.1% intergenic (+28/‑81) accC → / → yhdT acetyl‑CoA carboxylase, biotin carboxylase subunit/DUF997 family putative inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,407,2940TA11.1% 197.1 / 14.7 63intergenic (+28/‑81)accC/yhdTacetyl‑CoA carboxylase, biotin carboxylase subunit/DUF997 family putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (27/29);  new base A (2/5);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.76e-01

CGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGT  >  NC_000913/3407149‑3407438
                                                                                                                                                 |                                                                                                                                                
cGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAtttt                                                                                                                                                <  1:914095/148‑1 (MQ=255)
cGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTc                                                                                                                                               <  4:828382/149‑1 (MQ=255)
 gACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTcc                                                                                                                                              <  4:338023/149‑1 (MQ=255)
    ggTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCgg                                                                                                                                            <  1:818860/148‑1 (MQ=255)
    ggTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGc                                                                                                                                           <  4:896586/149‑1 (MQ=255)
        tCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCttt                                                                                                                                       <  3:241249/149‑1 (MQ=255)
        tCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGAAAATCCGGCCtt                                                                                                                                        <  2:298296/148‑1 (MQ=255)
         cAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCttt                                                                                                                                       >  1:796543/1‑148 (MQ=255)
           aaaCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGCACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGCCTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCcttctt                                                                                                                                     >  1:301805/1‑148 (MQ=255)
            aaCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                                                                                   <  3:745526/149‑1 (MQ=255)
             aCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                                                                                   >  2:742152/1‑148 (MQ=255)
             aCCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTtta                                                                                                                                   <  2:706981/148‑1 (MQ=255)
                   gTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTgg                                                                                                                            <  4:11338/149‑1 (MQ=255)
                        tCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGAt                                                                                                                        >  2:67861/1‑148 (MQ=255)
                               aTCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAAcc                                                                                                                >  4:22747/1‑149 (MQ=255)
                               aTCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAAc                                                                                                                 <  1:701430/148‑1 (MQ=255)
                                      tcatGAATGACGAGAACTTCCAGCATCGTAGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAATAATAACACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTCTTTATTACTGGGGATCGACAACCCCCATaa                                                                                                         <  4:758652/149‑1 (MQ=255)
                                       catGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATaa                                                                                                         >  1:140734/1‑148 (MQ=255)
                                       catGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATaa                                                                                                         <  2:140734/148‑1 (MQ=255)
                                       catGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATaa                                                                                                         <  1:28876/148‑1 (MQ=255)
                                          gaatgaCGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCCGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTa                                                                                                     <  4:512816/149‑1 (MQ=255)
                                                   gaACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCCTAAGGTACAATCCCCg                                                                                            >  4:529968/1‑149 (MQ=255)
                                                    aaCTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCg                                                                                            >  2:523920/1‑148 (MQ=255)
                                                            gCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTttcttc                                                                                    >  1:35016/1‑148 (MQ=255)
                                                             cATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCAGCTTTCTTCAc                                                                                  <  3:529968/149‑1 (MQ=255)
                                                                       cTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATaa                                                                                                         <  1:258708/116‑1 (MQ=255)
                                                                       cTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATaa                                                                                                         >  2:258708/1‑116 (MQ=255)
                                                                       cTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAggg                                                                         <  2:35016/148‑1 (MQ=255)
                                                                       cTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAggg                                                                         >  1:463149/1‑148 (MQ=255)
                                                                         aaCATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTttctt                                                                                     <  2:132989/134‑1 (MQ=255)
                                                                         aaCATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTttctt                                                                                     >  1:132989/1‑134 (MQ=255)
                                                                         aaCATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAg                                                                           <  1:708998/144‑1 (MQ=255)
                                                                         aaCATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAg                                                                           >  2:708998/1‑144 (MQ=255)
                                                                         aaCATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGAc                                                                       >  1:853488/1‑148 (MQ=255)
                                                                           cATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaa                                                                     >  1:876514/1‑148 (MQ=255)
                                                                             tCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaaa                                                                   >  2:6584/1‑148 (MQ=255)
                                                                             tCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaaa                                                                   >  2:8343/1‑148 (MQ=255)
                                                                               cACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACaaaaaa                                                                 <  1:440184/148‑1 (MQ=255)
                                                                                          gaAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAg                                                                           >  4:711708/1‑127 (MQ=255)
                                                                                          gaAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAg                                                                           <  3:711708/127‑1 (MQ=255)
                                                                                             aaaaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAgg                                                                          >  4:64958/1‑125 (MQ=255)
                                                                                             aaaaCTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAgg                                                                          <  3:64958/125‑1 (MQ=255)
                                                                                                  tCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCGTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCa                                              <  2:301805/148‑1 (MQ=255)
                                                                                                    ggTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGt                                                 <  4:760385/143‑1 (MQ=255)
                                                                                                    ggTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGt                                                 >  3:760385/1‑143 (MQ=255)
                                                                                                         tCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCTGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCAt                                       >  2:434925/1‑148 (MQ=255)
                                                                                                           aGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATaa                                     <  1:6584/148‑1 (MQ=255)
                                                                                                             gAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGa                                                             >  3:506223/1‑122 (MQ=255)
                                                                                                             gAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGa                                                             <  4:506223/122‑1 (MQ=255)
                                                                                                             gAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAag                                   <  2:876514/148‑1 (MQ=255)
                                                                                                             gAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAag                                   >  1:143215/1‑148 (MQ=255)
                                                                                                                    aaGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGgcgcg                            <  2:463149/148‑1 (MQ=255)
                                                                                                                            tAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCgcgctgggcgc                    >  2:400960/1‑148 (MQ=255)
                                                                                                                              aaGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCgcgctgggcgctg                  <  1:742152/148‑1 (MQ=255)
                                                                                                                                 cGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAag                                   <  2:424650/128‑1 (MQ=255)
                                                                                                                                 cGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAag                                   >  1:424650/1‑128 (MQ=255)
                                                                                                                                      aaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCgcgctgggcgc                    >  4:389771/1‑138 (MQ=255)
                                                                                                                                      aaaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCgcgctgggcgc                    <  3:389771/138‑1 (MQ=255)
                                                                                                                                       aaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGc                                           >  4:183588/1‑114 (MQ=255)
                                                                                                                                       aaGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGc                                           <  3:183588/114‑1 (MQ=255)
                                                                                                                                        aGGCCGGAAAATCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGaca                                                           <  3:478352/97‑1 (MQ=255)
                                                                                                                                         ggCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCt       <  1:8343/148‑1 (MQ=255)
                                                                                                                                           ccGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCtttt    >  3:240621/1‑149 (MQ=255)
                                                                                                                                             ggATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGt  >  4:147659/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CGACGGTATCAAAACCAACGTTGATCTGCAGATCCGCATCATGAATGACGAGAACTTCCAGCATGGTGGCACTAACATCCACTATCTGGAGAAAAAACTCGGTCTTCAGGAAAAATAAGACTGCTAAAGCGTCAAAAGGCCGGATTTTCCGGCCTTTTTTATTACTGGGGATCGACAACCCCCATAAGGTACAATCCCCGCTTTCTTCACCCATCAGGGACAAAAAATGGACACTCGTTTTGTTCAGGCCCATAAAGAGGCGCGCTGGGCGCTGGGGCTGACCCTTTTGT  >  NC_000913/3407149‑3407438

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