Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 37,311 T→A 5.3% T172S (ACT→TCT) ‡ caiC ← carnitine‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091337,3110TA 5.3% 294.9 / 5.6 94T172S (ACT→TCT) ‡caiCcarnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (43/46);  new base A (3/2);  total (46/48)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

CGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATAC  >  NC_000913/37168‑37459
                                                                                                                                                |                                                                                                                                                    
cGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTTGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAg                                                                                                                                                <  2:627196/150‑1 (MQ=255)
 gTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAgc                                                                                                                                               >  2:734167/1‑150 (MQ=255)
    aGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGCTATCGACTGACGATacggc                                                                                                                                             <  2:1316235/150‑5 (MQ=255)
    aGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAgcggc                                                                                                                                            >  1:795041/1‑150 (MQ=255)
     ggTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAgcggcg                                                                                                                                           >  1:951469/1‑150 (MQ=255)
     ggTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAgcggcg                                                                                                                                           >  2:1184068/1‑150 (MQ=255)
        aGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGCTATCGACTGACGATA‑cggcggaa                                                                                                                                        <  1:600889/150‑3 (MQ=255)
          aCGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                      <  2:1483894/150‑1 (MQ=255)
          aCGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                                      <  2:1552683/150‑1 (MQ=255)
           cGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCAt                                                                                                                                     >  1:1546721/1‑150 (MQ=255)
             tCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAg                                                                                                                                   <  2:951469/150‑1 (MQ=255)
              cATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGc                                                                                                                                  <  1:1407352/150‑1 (MQ=255)
                   tCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAg                                                                                                                             <  1:1099849/150‑1 (MQ=255)
                   cctggCAGTGTGCACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAg                                                                                                                             >  2:736251/5‑150 (MQ=255)
                    cACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAgg                                                                                                                            >  1:681701/1‑150 (MQ=255)
                     aCGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGt                                                                                                                           >  2:902036/1‑150 (MQ=255)
                     aCGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGt                                                                                                                           <  1:428514/150‑1 (MQ=255)
                       gCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTTGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTgg                                                                                                                         <  1:544810/150‑1 (MQ=255)
                         aGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCa                                                                                                                       >  2:1206096/1‑150 (MQ=255)
                         aGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCa                                                                                                                       <  2:1273025/150‑1 (MQ=255)
                          gTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAg                                                                                                                      >  2:213618/1‑150 (MQ=255)
                            gCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgt                                                                                                                    >  1:1568660/1‑150 (MQ=255)
                            gCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgt                                                                                                                    >  1:1378951/1‑150 (MQ=255)
                            gCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgt                                                                                                                    >  1:700947/1‑150 (MQ=255)
                            gCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgt                                                                                                                    >  2:192365/1‑150 (MQ=255)
                             cacaCTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgtt                                                                                                                   <  2:1141288/150‑1 (MQ=255)
                               cacTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgttgt                                                                                                                 >  1:1580152/1‑150 (MQ=255)
                                acTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgttgtt                                                                                                                <  2:397035/150‑1 (MQ=255)
                                 cTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGttgttg                                                                                                               >  2:1242168/1‑150 (MQ=255)
                                     cAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGAttt                                                                                                           >  1:903450/1‑150 (MQ=255)
                                       ggCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGAttttt                                                                                                         >  1:167973/1‑150 (MQ=255)
                                        gCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTc                                                                                                        <  1:1597252/150‑1 (MQ=255)
                                         cGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCa                                                                                                       <  2:167973/150‑1 (MQ=255)
                                          ggAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAg                                                                                                      <  2:143532/150‑1 (MQ=255)
                                             gTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTg                                                                                                   >  1:637973/1‑150 (MQ=255)
                                             gTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTg                                                                                                   <  1:428581/150‑1 (MQ=255)
                                               aaTATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAg                                                                                                 <  1:213618/150‑1 (MQ=255)
                                                atatCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGt                                                                                                <  1:734167/150‑1 (MQ=255)
                                                 tatCCAGCGAATCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGGGAAGAGAATTTCCGCCG‑TATCTTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTa                                                                                               <  2:708084/150‑1 (MQ=255)
                                                  atCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTaa                                                                                              <  2:681701/150‑1 (MQ=255)
                                                  atCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTaa                                                                                              >  2:1543734/1‑150 (MQ=255)
                                                  atCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTaa                                                                                              >  2:996157/1‑150 (MQ=255)
                                                  atCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTaa                                                                                              >  1:1613253/1‑150 (MQ=255)
                                                     cAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACg                                                                                           >  1:1032832/1‑150 (MQ=255)
                                                     cAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACg                                                                                           >  2:891467/1‑150 (MQ=255)
                                                     cAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACg                                                                                           <  1:1235142/150‑1 (MQ=255)
                                                     cAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACg                                                                                           <  2:1546721/150‑1 (MQ=255)
                                                          aaGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTc                                                                                      >  1:534452/1‑150 (MQ=255)
                                                             cgcAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACg                                                                                   <  1:829960/150‑1 (MQ=255)
                                                             cgcAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACg                                                                                   >  1:1234185/1‑150 (MQ=255)
                                                                  gTTGTAATGGGTAATGACCACTCGTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGATAATTTCCGCCG‑TATAGTCAGTCGATAGCGGCGGTTCATATCACAAGGTGGCAGGTTTTTGATTTTTCAGTTGAGTAAACGAACTCACGCcatc                                                                              <  2:530364/150‑1 (MQ=255)
                                                                  gTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCcatc                                                                              <  1:1391588/150‑1 (MQ=255)
                                                                   ttGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCcatca                                                                             <  2:903450/150‑1 (MQ=255)
                                                                         tGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCg                                                                       >  1:347415/1‑150 (MQ=255)
                                                                            gTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGa                                                                    >  1:693240/1‑150 (MQ=255)
                                                                               aTCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGt                                                                 <  1:1543734/150‑1 (MQ=255)
                                                                                tCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTg                                                                >  1:1008257/1‑150 (MQ=255)
                                                                                 caccacACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGc                                                               <  2:1580152/150‑1 (MQ=255)
                                                                                  accacACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGcc                                                              <  2:1085007/150‑1 (MQ=255)
                                                                                   ccacACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCca                                                             <  2:1032832/150‑1 (MQ=255)
                                                                                       acCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATc                                                         <  2:812055/150‑1 (MQ=255)
                                                                                       acCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATc                                                         <  1:1093974/150‑1 (MQ=255)
                                                                                        ccTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCt                                                        >  1:235908/1‑150 (MQ=255)
                                                                                        ccTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCt                                                        <  1:1642494/150‑1 (MQ=255)
                                                                                         cTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTg                                                       <  1:902036/150‑1 (MQ=255)
                                                                                              ggtcggGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAgg                                                  <  1:1029109/150‑1 (MQ=255)
                                                                                                  gggAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCaaa                                              <  1:1634801/150‑1 (MQ=255)
                                                                                                   ggAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAAt                                             <  2:637973/150‑1 (MQ=255)
                                                                                                      ggtggtGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAAtgtg                                          >  1:646876/1‑150 (MQ=255)
                                                                                                         ggtgCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCg                                       <  2:1631876/150‑1 (MQ=255)
                                                                                                         ggtgCCGGAGGTGAAGAGAATTTCCGCCG‑CATCGTCAGTCGATAGAGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTTAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCg                                       >  2:1049786/1‑150 (MQ=255)
                                                                                                             ccGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAAt                                   <  1:790594/150‑1 (MQ=255)
                                                                                                                          gaATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCtt                      >  2:621039/1‑150 (MQ=255)
                                                                                                                             tttCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCct                   >  1:1139532/1‑150 (MQ=255)
                                                                                                                             tttCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCct                   >  2:625869/1‑150 (MQ=255)
                                                                                                                              ttCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCctg                  >  2:1271244/1‑150 (MQ=255)
                                                                                                                              ttCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCctg                  >  1:1033167/1‑150 (MQ=255)
                                                                                                                              ttCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCctg                  >  2:1069302/1‑150 (MQ=255)
                                                                                                                               tCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCctgc                 <  1:736251/150‑1 (MQ=255)
                                                                                                                               tCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCctgc                 >  2:223997/1‑150 (MQ=255)
                                                                                                                                 cgccgctATCGACTGACGATA‑CGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCctgctg               >  2:1282585/7‑150 (MQ=255)
                                                                                                                                  gccgctATCGACTGACGATA‑CGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGa              >  1:1391623/6‑150 (MQ=255)
                                                                                                                                  gccgctATCGACTGACGATA‑CGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGa              >  1:1385957/6‑150 (MQ=255)
                                                                                                                                  gccg‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGa              >  1:737975/1‑150 (MQ=255)
                                                                                                                                          cGTCAGTCGATAGCGGCGGTGCTTTGCACAAGGTGGCTTGTTGTTGATTTTTGAGTTGAGTAAACTAACTCACGCCATGATGAGGGGTATTTGGCACTTGTGTCATTCAAATGTGCCGCATTTTAGTGGCTTGTTCGTTATGAATGTGtt       <  2:1033167/150‑1 (MQ=255)
                                                                                                                                          cGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGtt       <  2:646876/150‑1 (MQ=255)
                                                                                                                                          cGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGtt       >  2:1423633/1‑150 (MQ=255)
                                                                                                                                          cGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGtt       <  2:1568660/150‑1 (MQ=255)
                                                                                                                                           gTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTg      <  1:1242168/150‑1 (MQ=255)
                                                                                                                                           gTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTg      <  2:1545145/150‑1 (MQ=255)
                                                                                                                                           gTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTg      >  2:1308776/1‑150 (MQ=255)
                                                                                                                                             cAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGAt    <  1:891467/150‑1 (MQ=255)
                                                                                                                                              aGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATa   <  1:1069302/150‑1 (MQ=255)
                                                                                                                                               gTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATAc  >  2:972427/1‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATAC  >  NC_000913/37168‑37459

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: ‑