Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 35,462 T→C 6.1% Y234C (TAT→TGT)  caiD ← crotonobetainyl‑CoA hydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091335,4620TC 6.1% 194.1 / 2.9 65Y234C (TAT→TGT) caiDcrotonobetainyl‑CoA hydratase
Reads supporting (aligned to +/- strand):  ref base T (28/33);  new base C (2/2);  total (30/35)
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

TGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTAACCAGCTGCTGAGCCAGTTCGCGGGCGTTATCCATCAG  >  NC_000913/35313‑35610
                                                                                                                                                     |                                                                                                                                                    
tGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGAt                                                                                                                                                      >  1:874701/1‑150 (MQ=255)
  gACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAg                                                                                                                                                    >  1:522989/1‑150 (MQ=255)
         gCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCa                                                                                                                                             <  1:1045245/150‑1 (MQ=255)
         gCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCa                                                                                                                                             <  1:868103/150‑1 (MQ=255)
         gCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCa                                                                                                                                             <  1:1206442/150‑1 (MQ=255)
            gTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAaca                                                                                                                                          <  1:949355/150‑1 (MQ=255)
             tCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAacac                                                                                                                                         >  2:318068/1‑150 (MQ=255)
                   gAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCt                                                                                                                                   >  2:1272624/1‑150 (MQ=255)
                        accGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGaa                                                                                                                              >  2:492252/1‑150 (MQ=255)
                         ccGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAat                                                                                                                             <  2:650849/150‑1 (MQ=255)
                              aTTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGc                                                                                                                        >  2:137545/1‑150 (MQ=255)
                               ttAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCg                                                                                                                       <  1:648256/150‑1 (MQ=255)
                               ttAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCg                                                                                                                       <  2:21171/150‑1 (MQ=255)
                                tAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGa                                                                                                                      >  1:639648/1‑150 (MQ=255)
                                tAACCCCTCAAACGCGTAATAGCTAAAGTATCTAAACGTCCTTTCCACACCGGATACCGCATATCGGCAAACGCCAGCGACCCTACAATCCCATATTCCCAATACAGAACCGATGTATAGTGTATCAACAACCCGCTCCCAATAAAGCGa                                                                                                                      >  2:58820/1‑150 (MQ=255)
                                 aaCCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAac                                                                                                                                           <  1:38449/128‑1 (MQ=255)
                                 aaCCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAac                                                                                                                                           >  2:38449/1‑128 (MQ=255)
                                       tCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGctt                                                                                                               <  2:17473/150‑1 (MQ=255)
                                       tCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGctt                                                                                                               >  2:1051812/1‑150 (MQ=255)
                                       tCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGACAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGctt                                                                                                               >  2:38975/1‑150 (MQ=255)
                                       tCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGACAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGctt                                                                                                               >  2:44168/1‑150 (MQ=255)
                                        cAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGcttc                                                                                                              >  2:319808/1‑150 (MQ=255)
                                        cAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGcttc                                                                                                              <  1:894809/150‑1 (MQ=255)
                                            ggCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTtcttct                                                                                                          >  2:900855/1‑150 (MQ=255)
                                               gTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTAcc                                                                                                       <  2:71122/150‑1 (MQ=255)
                                               gTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTAcc                                                                                                       <  2:750824/150‑1 (MQ=255)
                                                  aTAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGc                                                                                                    >  2:1006181/1‑150 (MQ=255)
                                                                  aaCGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGc                                                                                    <  1:137545/150‑1 (MQ=255)
                                                                     gTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGt                                                                                 <  1:619902/150‑1 (MQ=255)
                                                                      tCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTa                                                                                <  1:889028/150‑1 (MQ=255)
                                                                          ttCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCTATGGCATCTTCCGAATGCAGAACCGATTGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCGGGTGGTGCGGTAGAtc                                                                            <  2:470347/150‑1 (MQ=255)
                                                                            ccACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGAtctc                                                                          >  2:1109310/1‑150 (MQ=255)
                                                                              acacCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCtt                                                                        <  1:318068/150‑1 (MQ=255)
                                                                                acCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTc                                                                      <  1:869447/150‑1 (MQ=255)
                                                                                 ccGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCa                                                                     <  1:855918/150‑1 (MQ=255)
                                                                                 ccGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCa                                                                     <  1:528075/150‑1 (MQ=255)
                                                                                 ccGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCa                                                                     >  2:1028135/1‑150 (MQ=255)
                                                                                      tCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGGGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGcgcc                                                                <  1:738979/150‑1 (MQ=255)
                                                                                      tCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGACAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGcgcc                                                                <  1:44168/150‑1 (MQ=255)
                                                                                      tCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGACAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGcgcc                                                                <  1:38975/150‑1 (MQ=255)
                                                                                      gcgcgCTTCTCGGCAAACGCCAGCGGCCCTTCCCGGGCATCTTCCGAATGCAGAACCGCTGGATAGGGTGTCAACACGCCGCTGCGAATATAGCGATACGCTTCTGCTACCGGCATTTCGCTGGGGGTGCGGTAGATCTCTTTCAGcgcc                                                                <  1:743930/149‑1 (MQ=255)
                                                                                          gcTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCaa                                                            <  1:58820/150‑1 (MQ=255)
                                                                                             tctcGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCg                                                         <  1:117057/150‑1 (MQ=255)
                                                                                               tcGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATAATATTCTTCTACCTGCATTTCGATGGTGTTGCTTTATATCTCTTTCAGTGCAGCAATCGcc                                                       <  2:639648/150‑1 (MQ=255)
                                                                                                    aaaCGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCgg                                                  >  2:44389/1‑150 (MQ=255)
                                                                                                     aaCGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCggg                                                 >  1:1413444/1‑150 (MQ=255)
                                                                                                           aGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTg                                           <  1:1051812/150‑1 (MQ=255)
                                                                                                           aGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTg                                           >  2:4617/1‑150 (MQ=255)
                                                                                                            gCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGt                                          <  1:341898/150‑1 (MQ=255)
                                                                                                            gCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGt                                          <  1:81414/150‑1 (MQ=255)
                                                                                                            gCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGt                                          <  2:874701/150‑1 (MQ=255)
                                                                                                              ggCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTa                                        >  1:771413/1‑150 (MQ=255)
                                                                                                               gCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTaa                                       <  1:492252/150‑1 (MQ=255)
                                                                                                                 ccTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTAAcc                                     <  1:319808/150‑1 (MQ=255)
                                                                                                                   ttCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTAACCAg                                   >  1:679641/1‑150 (MQ=255)
                                                                                                                         ggCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTAACCAGctgctg                             >  1:621609/1‑150 (MQ=255)
                                                                                                                               ttCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTAACCAGCTGCTGAGCCAg                       >  2:1385340/1‑150 (MQ=255)
                                                                                                                               ttCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTAACCAGCTGCTGAGCCAg                       >  2:1089932/1‑150 (MQ=255)
                                                                                                                                tCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTAACCAGCTGCTGAGCCAGt                      >  1:1189530/1‑150 (MQ=255)
                                                                                                                                tCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTAACCAGCTGCTGAGCCAGt                      >  1:263327/1‑150 (MQ=255)
                                                                                                                                      tGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTAACCAGCTGCTGAGCCAGTTCGCgg                >  1:103506/1‑150 (MQ=255)
                                                                                                                                       gCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTAACCAGCTGCTGAGCCAGTTCGCggg               >  2:1131324/1‑150 (MQ=255)
                                                                                                                                         aGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTAACCAGCTGCTGAGCCAGTTCGCGGGCg             >  2:102342/1‑150 (MQ=255)
                                                                                                                                            aCCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTAACCAGCTGCTGAGCCAGTTCGCGGGCGTTa          <  1:102342/150‑1 (MQ=255)
                                                                                                                                               gATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTAACCAGCTGCTGAGCCAGTTCGCGGGCGTTATcc       <  1:1272624/150‑1 (MQ=255)
                                                                                                                                                    aTAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTAACCAGCTGCTGAGCCAGTTCGCGGGCGTTATCCATCAg  >  2:1244322/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGGACAAACGCCGTCGGGTGAACCACCGGAATTAACCCCTCAAAGGCGTAATAGCTCACGGTTGTTAACGTCCTTTCCACACCGGATCGCGCTTCTCGGCAAACGCCAGCGGCCCTTCAATGGCATCTTCCGAATGCAGAACCGATGGATAGTGTTTCAACACGCCGCTGCGAATATAGCGATACGCTTCTTCTACCGGCATTTCGCTGGTGGTGCGGTAGATCTCTTTCAGCGCCGCAATCGCCAGCGGGGCGCTGTTAACCAGCTGCTGAGCCAGTTCGCGGGCGTTATCCATCAG  >  NC_000913/35313‑35610

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