Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 37,311 T→A 9.4% 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 9.4% 261.6 / 12.2 92T172S (ACT→TCT) ‡caiCcarnitine‑CoA ligase
Reads supporting (aligned to +/- strand):  ref base T (40/43);  new base A (3/6);  total (43/49)
Fisher's exact test for biased strand distribution p-value = 4.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATA  >  NC_000913/37162‑37458
                                                                                                                                                      |                                                                                                                                                   
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:787504/150‑4 (MQ=255)
               gACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCGCTATCGACTGACGATA‑cggcggaaa                                                                                                                                      <  2:781631/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:904613/150‑1 (MQ=255)
                                                     aaTATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAg                                                                                                <  2:239856/150‑1 (MQ=255)
                                                      atatCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGt                                                                                               <  2:1626690/150‑1 (MQ=255)
                                                      atatCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGt                                                                                               >  1:1475917/1‑150 (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                                                                                                >  2:466596/1‑144 (MQ=255)
                                                           cAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAg                                                                                                <  1:466596/144‑1 (MQ=255)
                                                           cAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACg                                                                                          >  2:61841/1‑150 (MQ=255)
                                                           cAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACg                                                                                          <  2:895412/150‑1 (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:810672/150‑1 (MQ=255)
                                                                                gggTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCgg                                                                     <  1:299991/150‑1 (MQ=255)
                                                                                gggTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCgg                                                                     >  2:89033/1‑150 (MQ=255)
                                                                                gggTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCgg                                                                     <  1:749850/150‑1 (MQ=255)
                                                                                  gTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCa                                                                                                      <  2:933109/115‑1 (MQ=255)
                                                                                  gTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCa                                                                                                      >  1:933109/1‑115 (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                                      >  1:1630297/1‑150 (MQ=255)
                                                                                                               ggtgCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCg                                      <  2:658770/150‑1 (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)
                                                                                                                                                    aGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATa  <  2:73531/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
CATTACCGTCAGGTAGACGTCATCGTCACGCAGTGCACACTGCCAGGCGGAGTAATATCCAGCGAAGCGCAGGTTGTAATGGGTAATCACCACACCTTTCGGTCGGGAGGTGGTGCCGGAGGTGAAGAGAATTTCCGCCG‑TATCGTCAGTCGATAGCGGCGGTGCATAGCACAAGGTGGCAGGTTGTTGATTTTTCAGTTGAGTAAACGAACTCACGCCATCATCAGCGGGAAGTGCCACATCTGTCAGGCAAATGTGCCGCAATTGAGTGGCATCTTCCTGCTGAATCTGTTGATA  >  NC_000913/37162‑37458

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