Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 37,307 T→A 9.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 9.4% 275.8 / 15.6 96D173V (GAC→GTC) caiCcarnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (41/46);  new base A (3/6);  total (44/52)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.59e-01

CAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTG  >  NC_000913/37158‑37452
                                                                                                                                                      |                                                                                                                                                 
cAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGt                                                                                                                                                   >  1:658770/1‑150 (MQ=255)
cAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGt                                                                                                                                                   <  2:662429/150‑1 (MQ=255)
 aGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCGAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTc                                                                                                                                                  <  2:1153143/150‑1 (MQ=255)
  ggCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCa                                                                                                                                                 >  2:157692/1‑150 (MQ=255)
  ggCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCa                                                                                                                                                 <  2:1129475/150‑1 (MQ=255)
   gCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTATCTGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAg                                                                                                                                                <  2:1625743/150‑1 (MQ=255)
    cATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGt                                                                                                                                               <  1:1463427/150‑1 (MQ=255)
         ccGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATa                                                                                                                                          >  1:68740/1‑150 (MQ=255)
              aGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGCTATCGACTGACGATacggc                                                                                                                                      <  1:635114/150‑5 (MQ=255)
              aGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAgcggc                                                                                                                                     <  2:301574/150‑1 (MQ=255)
               ggTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAgcggcg                                                                                                                                    >  1:1618974/1‑150 (MQ=255)
               ggTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTGGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAgcggcg                                                                                                                                    >  2:1509160/1‑150 (MQ=255)
               ggTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGGTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGCTATCGACTGACGATa‑cggcg                                                                                                                                    <  2:810457/150‑1 (MQ=255)
                 tAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGCTATCGACTGACGATA‑cggcgga                                                                                                                                  <  1:1224942/150‑2 (MQ=255)
                 tAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGCTATCGACTGACGATA‑cggcgga                                                                                                                                  <  1:910887/150‑2 (MQ=255)
                   gACGTCATCGTCACGCAGTGCACACTGCCATGCGGAGTAATATCCTGCGAAGCGCAGGTTGTTATGGGTAATGACGACTCCTTTCTGTCGGGAGGTTGTGCCGGAGGTGATGAGAATTTCGGCCT‑TATCGTAATTCGATAGCGGCGGTGc                                                                                                                                <  2:734375/150‑1 (MQ=255)
                   gACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGCTATCGACTGACGATA‑cggcggaaa                                                                                                                                <  2:781631/150‑4 (MQ=255)
                   gACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGCTATCGACTGACGATA‑cggcggaaa                                                                                                                                <  2:787504/150‑4 (MQ=255)
                   gACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGc                                                                                                                                >  1:361975/1‑150 (MQ=255)
                    aCGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCa                                                                                                                               <  1:824861/150‑1 (MQ=255)
                     cGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCAt                                                                                                                              <  2:52359/150‑1 (MQ=255)
                      gTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATa                                                                                                                             <  1:490321/150‑1 (MQ=255)
                       tCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAg                                                                                                                            >  1:688570/1‑150 (MQ=255)
                        cATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGc                                                                                                                           >  1:904613/1‑150 (MQ=255)
                           cGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGcaca                                                                                                                        >  2:851301/1‑150 (MQ=255)
                           cGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGcaca                                                                                                                        >  2:1270344/1‑150 (MQ=255)
                              cACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAgg                                                                                                                     >  2:387920/1‑150 (MQ=255)
                               aCGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGt                                                                                                                    >  2:749850/1‑150 (MQ=255)
                               aCGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGt                                                                                                                    <  2:1355898/150‑1 (MQ=255)
                                  cAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGc                                                                                                                 >  1:640497/1‑150 (MQ=255)
                                      gCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgt                                                                                                             >  2:1554236/1‑150 (MQ=255)
                                      gCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgt                                                                                                             >  2:85009/1‑150 (MQ=255)
                                        acacTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATTGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgttg                                                                                                           <  1:207303/150‑1 (MQ=255)
                                         cacTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGgttgt                                                                                                          >  2:1058450/1‑150 (MQ=255)
                                           cTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGTTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGttgttg                                                                                                        <  1:1270344/150‑1 (MQ=255)
                                             gCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGAt                                                                                                      >  2:538960/1‑150 (MQ=255)
                                                aGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGAtttt                                                                                                   <  2:465422/150‑1 (MQ=255)
                                                 ggCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGAttttt                                                                                                  >  2:144605/1‑150 (MQ=255)
                                                  gCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTc                                                                                                 <  2:688570/150‑1 (MQ=255)
                                                   cGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCa                                                                                                <  1:638528/150‑1 (MQ=255)
                                                    ggAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAg                                                                                               <  1:567546/150‑1 (MQ=255)
                                                    ggAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAg                                                                                               <  2:68740/150‑1 (MQ=255)
                                                       gTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTg                                                                                            <  1:1501094/150‑1 (MQ=255)
                                                         aaTATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAg                                                                                          <  2:239856/150‑1 (MQ=255)
                                                         aaTATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAg                                                                                          <  2:904613/150‑1 (MQ=255)
                                                          atatCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGt                                                                                         >  1:1475917/1‑150 (MQ=255)
                                                          atatCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGt                                                                                         <  2:1626690/150‑1 (MQ=255)
                                                            atCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTaa                                                                                       >  2:34320/1‑150 (MQ=255)
                                                            atCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTaa                                                                                       >  2:592918/1‑150 (MQ=255)
                                                              ccAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAAc                                                                                     >  2:1347988/1‑150 (MQ=255)
                                                               cAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAg                                                                                          <  1:466596/144‑1 (MQ=255)
                                                               cAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAg                                                                                          >  2:466596/1‑144 (MQ=255)
                                                               cAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACg                                                                                    <  2:895412/150‑1 (MQ=255)
                                                               cAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACg                                                                                    >  2:61841/1‑150 (MQ=255)
                                                                    aaGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTc                                                                               <  2:195552/150‑1 (MQ=255)
                                                                        gcAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGc                                                                           >  1:1261668/1‑150 (MQ=255)
                                                                            gTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCcatc                                                                       <  2:908846/150‑1 (MQ=255)
                                                                             ttGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCcatca                                                                      >  2:965970/1‑150 (MQ=255)
                                                                               gTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGcc                                                                          >  1:35459/1‑144 (MQ=255)
                                                                               gTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGcc                                                                          <  2:35459/144‑1 (MQ=255)
                                                                               gTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCatcatc                                                                    <  2:640497/150‑1 (MQ=255)
                                                                                 aaTGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAg                                                                  >  1:500247/1‑150 (MQ=255)
                                                                                 aaTGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTCCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACTAACTCACGCCATCATCAg                                                                  >  2:1467343/1‑150 (MQ=255)
                                                                                   tGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCg                                                                <  2:697/150‑1 (MQ=255)
                                                                                    gggTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCgg                                                               >  2:89033/1‑150 (MQ=255)
                                                                                    gggTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCgg                                                               <  2:810672/150‑1 (MQ=255)
                                                                                    gggTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCgg                                                               <  1:299991/150‑1 (MQ=255)
                                                                                    gggTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCgg                                                               <  1:749850/150‑1 (MQ=255)
                                                                                      gTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCa                                                                                                >  1:933109/1‑115 (MQ=255)
                                                                                      gTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCa                                                                                                <  2:933109/115‑1 (MQ=255)
                                                                                             ccacACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCca                                                      >  1:1428211/1‑150 (MQ=255)
                                                                                              cacaCCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCcac                                                     >  1:251320/1‑150 (MQ=255)
                                                                                                      tcggtcggGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCa                                             <  1:538960/150‑1 (MQ=255)
                                                                                                          tcggGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCa                                         <  1:1554236/150‑1 (MQ=255)
                                                                                                             ggAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAAt                                      <  1:1007535/150‑1 (MQ=255)
                                                                                                                ggtggtGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAAtgtg                                   >  1:1463181/1‑150 (MQ=255)
                                                                                                                   ggtgCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCg                                <  2:658770/150‑1 (MQ=255)
                                                                                                                   ggtgCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCg                                >  1:1630297/1‑150 (MQ=255)
                                                                                                                      gCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCaa                             >  2:1213394/1‑150 (MQ=255)
                                                                                                                       ccGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGAGAAGTGCCACATCTGTCAGGCAACTGTGCCGCAAt                            <  1:1273892/150‑1 (MQ=255)
                                                                                                                        cGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAAtt                           >  1:515250/1‑150 (MQ=255)
                                                                                                                         ggAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTg                          <  2:500247/150‑1 (MQ=255)
                                                                                                                              tGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTgg                     >  2:627380/1‑150 (MQ=255)
                                                                                                                              tGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTgg                     <  1:486675/150‑1 (MQ=255)
                                                                                                                               gAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGc                    >  1:98606/1‑150 (MQ=255)
                                                                                                                                 agagAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCAt                  >  1:172805/1‑150 (MQ=255)
                                                                                                                                  cagaATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATc                 <  2:1041798/149‑1 (MQ=255)
                                                                                                                                   agaATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCGCATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCt                >  1:224746/1‑150 (MQ=255)
                                                                                                                                      aTTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTACCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTcc             <  1:495511/150‑1 (MQ=255)
                                                                                                                                       gcaccgccgcTATCGACTGACGATA‑CGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCct            >  2:151930/4‑150 (MQ=255)
                                                                                                                                        caccgccgcTATCGACTGACGATA‑CGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCctg           >  2:627507/3‑150 (MQ=255)
                                                                                                                                          ccgccg‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGAAGGTTGTTGATTTTTAAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATTTGCCGCAATTGAGTGGCATCTTCctgct         >  2:1605023/1‑150 (MQ=255)
                                                                                                                                            gccgctATCGACTGACGATA‑CGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGCGGCATCTTCCTGCTGa       >  1:1564528/6‑150 (MQ=255)
                                                                                                                                            gccg‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGa       <  1:89033/150‑1 (MQ=255)
                                                                                                                                               g‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATc    <  1:157692/150‑1 (MQ=255)
                                                                                                                                               g‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTCATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATc    <  1:85009/150‑1 (MQ=255)
                                                                                                                                               catATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAAt     <  2:41049/148‑1 (MQ=255)
                                                                                                                                                  aTCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTg  <  1:965970/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                 
CAGGCATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTG  >  NC_000913/37158‑37452

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