Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 37,307 T→A 5.4% D173V (GAC→GTC)  caiC ← carnitine‑CoA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091337,3070TA 5.4% 300.6 / 5.3 93D173V (GAC→GTC) caiCcarnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (43/45);  new base A (3/2);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.02e-01

AGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTG  >  NC_000913/37159‑37455
                                                                                                                                                     |                                                                                                                                                    
aGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTc                                                                                                                                                     <  2:1005772/150‑1 (MQ=255)
aGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTc                                                                                                                                                     >  2:1587583/1‑150 (MQ=255)
         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                                                                                                                >  2:192365/1‑150 (MQ=255)
                                     gCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgt                                                                                                                >  1:700947/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                                                     <  1:1093974/150‑1 (MQ=255)
                                                                                                acCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATc                                                     <  2:812055/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              >  1:1033167/1‑150 (MQ=255)
                                                                                                                                       ttCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCctg              >  2:1271244/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:1385957/6‑150 (MQ=255)
                                                                                                                                           gccgctATCGACTGACGATA‑CGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGa          >  1:1391623/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:1568660/150‑1 (MQ=255)
                                                                                                                                                   cGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGtt   >  2:1423633/1‑150 (MQ=255)
                                                                                                                                                    gTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTg  <  1:1242168/150‑1 (MQ=255)
                                                                                                                                                    gTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTg  <  2:1545145/150‑1 (MQ=255)
                                                                                                                                                    gTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTg  >  2:1308776/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTG  >  NC_000913/37159‑37455

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