Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 40,191 C→G 6.4% G66R (GGG→CGG)  caiA ← crotonobetaine reductase subunit II, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091340,1910CG 6.4% 201.0 / 5.6 63G66R (GGG→CGG) caiAcrotonobetaine reductase subunit II, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (33/26);  new base G (2/2);  total (35/28)
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.06e-01

CGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCG  >  NC_000913/40047‑40336
                                                                                                                                                |                                                                                                                                                 
cgcgGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCcgcg                                                                                                                                               >  2:249537/1‑149 (MQ=255)
 gcgGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTGACCCAGCCCGGCAACTGGAAAAGCCCATAGGTTGGTGCCCCCAGCCGTCCCAGCTCCATCCCCCCGTCGTCGAGCGTACCAAACCCCGCGt                                                                                                                                              >  6:179916/1‑149 (MQ=255)
     aaaGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTTGTGACCCCAGACGTACCAGCTACATCAACACGGCGGCGAGAGTAACAAACCCCGCGTCCAg                                                                                                                                          <  4:117968/149‑1 (MQ=255)
         ccATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGacca                                                                                                                                      >  7:76209/1‑149 (MQ=255)
         ccATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGacca                                                                                                                                      >  6:69906/1‑149 (MQ=255)
            tAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCg                                                                                                                                   >  5:205009/1‑149 (MQ=255)
            tAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCg                                                                                                                                   >  5:200260/1‑149 (MQ=255)
                       aTCGGCACTTCTGTGCCTTCGCGCAGCAAGGTGTTGAACCCGCCCGGCAACTCGTACAGGACATAGGATGGTGCCCCCAGCCGTCCCAGCTCCAGCCACACGGCGGCGAGAGTAACACACCCCGCGTCCAGACCAGCGTGCTCGTCAgg                                                                                                                        >  4:178023/1‑149 (MQ=255)
                            ctgttGTGTGCCTGCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTTGTACAGCACATAGGTTGGGGCCCCCAGCCGTCCCAGCGCCATCCACACGGCGGCGAGAGTAACACACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATc                                                                                                                    <  2:245932/145‑1 (MQ=255)
                            ctctTGTGTGCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAggg                                                                                                                       >  4:112538/1‑145 (MQ=255)
                            ctctTGTGTGCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAggg                                                                                                                       <  3:112538/145‑1 (MQ=255)
                            ctctTGTGTGCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGAt                                                                                                                     >  5:125616/1‑147 (MQ=255)
                            ctctTGTGTGCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATca                                                                                                                   <  3:134070/149‑1 (MQ=255)
                             tctTGTGTGCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATcag                                                                                                                  >  6:188195/1‑149 (MQ=255)
                                  tgtgCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGa                                                                                                              <  6:200260/148‑1 (MQ=255)
                                  tgtgCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGa                                                                                                              <  1:248971/148‑1 (MQ=255)
                                  tgtgCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGa                                                                                                              <  6:205009/148‑1 (MQ=255)
                                   gtgCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACt                                                                                                            >  8:117267/1‑149 (MQ=255)
                                       cTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGCCCACCGTGCTCTTCAGGGATCAGCAGACTGTc                                                                                                         >  5:340326/1‑148 (MQ=255)
                                        ttCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCGGCAGACTGTCGa                                                                                                       >  2:4489/1‑149 (MQ=255)
                                           gcgcAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATAc                                                                                                    <  7:117267/149‑1 (MQ=255)
                                                  aaGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATc                                                                                             >  1:327670/1‑149 (MQ=255)
                                                       gtTGAACCCGCCCGGCAACTGGTACAGCACTTAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGc                                                                                          >  7:282793/1‑147 (MQ=255)
                                                             cccgcccgGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAACCCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCtt                                                                                  >  3:176553/1‑149 (MQ=255)
                                                             cccgcccgGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCGCGCGTCCAGACCACCGTg                                                                                                                                 <  1:113092/102‑1 (MQ=255)
                                                             cccgcccgGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCGCGCGTCCAGACCACCGTg                                                                                                                                 >  2:108363/1‑102 (MQ=255)
                                                             cccgcccgGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCGCGCGTCCAGACCACCGTg                                                                                                                                 <  1:108363/102‑1 (MQ=255)
                                                             cccgcccgGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCGCGCGTCCAGACCACCGTg                                                                                                                                 >  2:113092/1‑102 (MQ=255)
                                                               cgcccgGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGCCCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTg                                                                                >  5:332653/1‑149 (MQ=255)
                                                                       aaCTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCa                                                                                        >  7:209138/1‑133 (MQ=255)
                                                                       aaCTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCa                                                                                        <  8:209138/133‑1 (MQ=255)
                                                                         cTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGtt                                                                      >  3:20064/1‑149 (MQ=255)
                                                                          tGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATcagc                                                                                                                 >  7:113897/1‑105 (MQ=255)
                                                                          tGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATcagc                                                                                                                 <  8:113897/105‑1 (MQ=255)
                                                                               cAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCggg                                                                 >  3:203372/1‑148 (MQ=255)
                                                                                 gCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAg                                                              <  1:249537/149‑1 (MQ=255)
                                                                                  cacaTAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGa                                                             <  4:176553/149‑1 (MQ=255)
                                                                                     aTAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCg                                                                   >  3:45752/1‑140 (MQ=255)
                                                                                     aTAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCg                                                                   <  4:45752/140‑1 (MQ=255)
                                                                                      tAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCAAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCt                                                         <  8:282793/149‑1 (MQ=255)
                                                                                          ttGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGGCa                                                     <  5:69906/149‑1 (MQ=255)
                                                                                            ggTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCAc                                                    <  5:188195/148‑1 (MQ=255)
                                                                                             gTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACgg                                                  <  6:125616/149‑1 (MQ=255)
                                                                                               gCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTc                                                <  2:327670/149‑1 (MQ=255)
                                                                                                cccccAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGCCTGTCGATCCCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCg                                               >  1:253508/1‑149 (MQ=255)
                                                                                                  cccAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCa                                             >  4:76889/1‑149 (MQ=255)
                                                                                                  cccAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCa                                             >  3:138364/1‑149 (MQ=255)
                                                                                                   cccgACGTCCCAGCTCCATCCACACGGCCGCGAGAGAAACAACCCCCGCGTCCAGACCACCGTGCGCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCAc                                            <  4:228028/146‑1 (MQ=255)
                                                                                                      gACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTc                                                                     <  4:260611/121‑1 (MQ=255)
                                                                                                      gACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTc                                                                     >  3:260611/1‑121 (MQ=255)
                                                                                                         gTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGc                                       >  3:138089/1‑148 (MQ=255)
                                                                                                              aGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCttt                                                                                 >  8:281223/1‑101 (MQ=255)
                                                                                                              aGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATa                                 <  8:76209/149‑1 (MQ=255)
                                                                                                                     tccaCACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCGGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTcc                          <  1:4489/149‑1 (MQ=255)
                                                                                                                     tccaCACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTcc                          >  6:156667/1‑149 (MQ=255)
                                                                                                                       cacaCGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTccc                         <  6:340326/148‑1 (MQ=255)
                                                                                                                                  gagaGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTgg             >  3:144702/1‑149 (MQ=255)
                                                                                                                                      gTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGCTCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCAt         <  3:76889/149‑1 (MQ=255)
                                                                                                                                      ctcaCAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGCTCAGCAGACTTTCGATACCCATATCCGCCAGTGCTGTGACAAAACGTTCCTGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCAt         <  5:179916/146‑1 (MQ=255)
                                                                                                                                         aCAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCgg                                                                  >  4:81473/1‑89 (MQ=255)
                                                                                                                                         aCAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCgg                                                                  <  3:81473/89‑1 (MQ=255)
                                                                                                                                            aaCCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTc   >  7:270501/1‑149 (MQ=255)
                                                                                                                                             aCCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTcg  >  1:12166/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CGCGGAAAGCCATAATTTTGTCGATCTGCTCTTGTGTGCCTTCGCGCAGGAAGGTGTTGAACCCGCCCGGCAACTGGTACAGCACATAGGTTGGTGCCCCCAGACGTCCCAGCTCCATCCACACGGCGGCGAGAGTAACAAACCCCGCGTCCAGACCACCGTGCTCTTCAGGGATCAGCAGACTGTCGATACCCATATCCGCCAGTGCTTTGACAAAACGTTCCGGGTAGACGCTGTCACGGTCGCACTCGGCAAAATAGGCCTCCCAGTTTTCGCTGGCCATCAGTTCG  >  NC_000913/40047‑40336

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