Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,825,763 C→A 91.2% A159A (GCG→GCT astC ← succinylornithine transaminase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,825,7630CA91.2% 826.0 / 66.8 274A159A (GCG→GCTastCsuccinylornithine transaminase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (20/4);  new base A (129/121);  total (149/125)
Fisher's exact test for biased strand distribution p-value = 2.56e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

Only 100 of 275 total aligned reads displayed.

ACCTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCG  >  CP009273/1825616‑1825908
                                                                                                                                                   |                                                                                                                                                 
     gTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGGTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAAt                                                                                                                                            >  1:4570818/1‑150 (MQ=255)
      tAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATc                                                                                                                                           >  2:3665709/1‑150 (MQ=255)
      tAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATc                                                                                                                                           >  1:2666793/1‑150 (MQ=255)
       aaaaaCGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATcc                                                                                                                                          <  1:1407415/150‑1 (MQ=255)
        aaaaCGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCt                                                                                                                                         >  2:2326626/1‑150 (MQ=255)
         aaaCGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCtg                                                                                                                                        <  2:1595975/150‑1 (MQ=255)
           aCGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCtgtg                                                                                                                                      >  2:3541892/1‑150 (MQ=255)
              cgTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAAt                                                                                                                                   >  1:3470327/1‑150 (MQ=255)
               gTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATa                                                                                                                                  >  2:2498294/1‑150 (MQ=255)
                     ggctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAggctgg                                                                                                                            <  2:1896612/150‑1 (MQ=255)
                       ctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAggctggct                                                                                                                          <  1:2404554/150‑1 (MQ=255)
                        tggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGgctggctg                                                                                                                         >  1:1741177/1‑150 (MQ=255)
                         ggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGc                                                                                                                        >  2:2814089/1‑150 (MQ=255)
                         ggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGc                                                                                                                        >  2:4426812/1‑150 (MQ=255)
                            accacACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCa                                                                                                                     <  2:2545201/150‑1 (MQ=255)
                              cacaCCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCAcc                                                                                                                   >  1:2625652/1‑150 (MQ=255)
                                cacCGCCGTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCg                                                                                                                 <  2:2119015/150‑1 (MQ=255)
                                 acagCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTAAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGAAGAATATCCGCCCGCAGTGGATCAAAATCCTGTGAATAGGCAGGCTTCCCAccccc                                                                                                                >  1:3536632/4‑148 (MQ=255)
                                   cGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCAc                                                                                                              >  2:885213/1‑150 (MQ=255)
                                      cTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGa                                                                                                           <  1:1284551/150‑1 (MQ=255)
                                           cccTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTa                                                                                                      <  1:135391/150‑1 (MQ=255)
                                              tGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAAc                                                                                                   >  2:1800713/1‑150 (MQ=255)
                                              tGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAAc                                                                                                   >  2:2359270/1‑150 (MQ=255)
                                                gATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATTGGCTGGCTGCCCACCCGCACTGACAGTAAACAg                                                                                                 >  1:1520519/1‑150 (MQ=255)
                                                gATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAg                                                                                                 >  2:2667516/1‑150 (MQ=255)
                                                gATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAg                                                                                                 >  2:274790/1‑150 (MQ=255)
                                                gATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAg                                                                                                 <  1:3254721/150‑1 (MQ=255)
                                                gATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAg                                                                                                 <  2:118988/150‑1 (MQ=255)
                                                gATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAg                                                                                                 >  1:2113293/1‑150 (MQ=255)
                                                 aTGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGc                                                                                                <  1:2359270/150‑1 (MQ=255)
                                                 aTGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGc                                                                                                >  2:937103/1‑150 (MQ=255)
                                                  tGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCg                                                                                               <  2:4113466/150‑1 (MQ=255)
                                                    ggTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTg                                                                                             >  1:4113649/1‑150 (MQ=255)
                                                      ttCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCg                                                                                           <  2:3586385/150‑1 (MQ=255)
                                                               cACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGaaa                                                                                  <  1:4527858/150‑1 (MQ=255)
                                                                aCCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAAc                                                                                 <  1:505989/150‑1 (MQ=255)
                                                                   gCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCa                                                                              >  2:1284250/1‑150 (MQ=255)
                                                                       aGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAtttt                                                                          <  1:754386/150‑1 (MQ=255)
                                                                          tAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGa                                                                       >  2:1507228/1‑150 (MQ=255)
                                                                           agagTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGaa                                                                      <  2:1387978/150‑1 (MQ=255)
                                                                                cgtcAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCa                                                                 <  2:2262038/150‑1 (MQ=255)
                                                                                 gtcAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCAc                                                                <  1:1098083/150‑1 (MQ=255)
                                                                                 gtcAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCAc                                                                >  1:3316501/1‑150 (MQ=255)
                                                                                  tcAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACg                                                               >  1:2291828/1‑150 (MQ=255)
                                                                                     aTCAGCGCGCTGGCAGAGTTAATATCGGTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGAAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATg                                                            >  2:4742769/1‑150 (MQ=255)
                                                                                        aGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCg                                                         >  1:1962810/1‑150 (MQ=255)
                                                                                          cgcgcTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGct                                                       >  1:1202838/1‑150 (MQ=255)
                                                                                             gcTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCtt                                                    <  2:2625652/150‑1 (MQ=255)
                                                                                                  cAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGc                                               <  2:1711663/150‑1 (MQ=255)
                                                                                                   agagTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCt                                              >  1:587336/1‑150 (MQ=255)
                                                                                                    gagTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTg                                             >  1:653765/1‑150 (MQ=255)
                                                                                                     agTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGc                                            >  2:3158874/1‑150 (MQ=255)
                                                                                                      gTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGcc                                           <  1:1800713/150‑1 (MQ=255)
                                                                                                      gTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGcc                                           <  1:1661616/150‑1 (MQ=255)
                                                                                                         aaTATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTa                                        >  1:2057269/1‑150 (MQ=255)
                                                                                                         aaTATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTa                                        >  1:3596229/1‑150 (MQ=255)
                                                                                                          atatCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAg                                       >  1:715002/1‑150 (MQ=255)
                                                                                                           tatCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGc                                      >  1:2687434/1‑150 (MQ=255)
                                                                                                            atCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCg                                     >  1:3802583/1‑150 (MQ=255)
                                                                                                              cGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGt                                   <  1:3013828/150‑1 (MQ=255)
                                                                                                              cGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGt                                   >  1:4000200/1‑150 (MQ=255)
                                                                                                                ttATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCg                                 <  1:811169/150‑1 (MQ=255)
                                                                                                                 tataTGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGt                                >  2:3875599/1‑150 (MQ=255)
                                                                                                                  atatGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTg                               <  1:1457345/150‑1 (MQ=255)
                                                                                                                   tatGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGa                              >  2:2075776/1‑150 (MQ=255)
                                                                                                                     tGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGc                            <  2:3737163/150‑1 (MQ=255)
                                                                                                                       cagcaTGACGAATATCCGCCGTCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGa                          <  1:680709/150‑1 (MQ=255)
                                                                                                                       cagcaTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGa                          <  2:4092581/150‑1 (MQ=255)
                                                                                                                       cagcaTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGa                          >  2:3460323/1‑150 (MQ=255)
                                                                                                                        agcaTGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGaa                         <  1:386299/150‑1 (MQ=255)
                                                                                                                        agcaTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGaa                         >  2:2857731/1‑150 (MQ=255)
                                                                                                                        agcaTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGaa                         >  1:78736/1‑150 (MQ=255)
                                                                                                                        agcaTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGaa                         >  1:2651565/1‑150 (MQ=255)
                                                                                                                          caTGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAtt                       >  2:2965802/1‑150 (MQ=255)
                                                                                                                          caTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAtt                       <  2:1073944/150‑1 (MQ=255)
                                                                                                                            tGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTa                     <  2:1429633/150‑1 (MQ=255)
                                                                                                                              aCGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTAcg                   >  1:2449797/1‑150 (MQ=255)
                                                                                                                              aCGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTAcg                   >  1:3637345/1‑150 (MQ=255)
                                                                                                                                gAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTAcgcg                 >  2:4256148/1‑150 (MQ=255)
                                                                                                                                gAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTAcgcg                 >  2:658216/1‑150 (MQ=255)
                                                                                                                                gAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTAcgcg                 >  1:4583972/1‑150 (MQ=255)
                                                                                                                                 aaTATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACgcgc                >  1:4304493/1‑150 (MQ=255)
                                                                                                                                   tatCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCa              >  2:803841/1‑150 (MQ=255)
                                                                                                                                     tCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGt            >  2:333636/1‑150 (MQ=255)
                                                                                                                                      ccgccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGtt           >  1:2133469/1‑150 (MQ=255)
                                                                                                                                       cgccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGttt          >  1:898000/1‑150 (MQ=255)
                                                                                                                                        gccgGCAGTGGAGCAAAATCCTGTGTATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGtttt         <  2:4657333/150‑1 (MQ=255)
                                                                                                                                        gccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGtttt         >  1:649391/1‑150 (MQ=255)
                                                                                                                                        gccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGtttt         <  1:2814089/150‑1 (MQ=255)
                                                                                                                                         ccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTa        <  1:369020/150‑1 (MQ=255)
                                                                                                                                          cgGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAg       >  2:94566/1‑150 (MQ=255)
                                                                                                                                          cgGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAg       >  2:288819/1‑150 (MQ=255)
                                                                                                                                          cgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAg       <  1:343101/150‑1 (MQ=255)
                                                                                                                                          cgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAg       >  2:4597102/1‑150 (MQ=255)
                                                                                                                                          cgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAg       >  1:2510160/1‑150 (MQ=255)
                                                                                                                                          cgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAg       >  2:2196546/1‑150 (MQ=255)
                                                                                                                                           ggCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAgc      <  1:4740657/150‑1 (MQ=255)
                                                                                                                                           ggCAGGGGAGCAAACTCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAgc      <  1:4410854/150‑1 (MQ=255)
                                                                                                                                            gCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAgcg     <  2:203615/150‑1 (MQ=255)
                                                                                                                                             cAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAgcgc    >  2:1758735/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
ACCTTGTAAAAACGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCG  >  CP009273/1825616‑1825908

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