Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,825,775 T→A 7.5% S155S (TCA→TCT astC ← succinylornithine transaminase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,825,7750TA 7.5% 619.0 / 31.2 201S155S (TCA→TCTastCsuccinylornithine transaminase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (95/91);  new base A (10/5);  total (105/96)
Fisher's exact test for biased strand distribution p-value = 2.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 205 total aligned reads displayed.

CGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTC  >  CP009273/1825628‑1825921
                                                                                                                                                   |                                                                                                                                                  
cgcgTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCtgt                                                                                                                                                    <  1:1055674/148‑1 (MQ=255)
     tGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATa                                                                                                                                               >  2:759421/1‑148 (MQ=255)
     tGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATa                                                                                                                                               <  2:1144794/148‑1 (MQ=255)
     tGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATa                                                                                                                                               >  1:1481575/1‑148 (MQ=255)
         ggctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAggct                                                                                                                                           >  2:3227522/1‑148 (MQ=255)
          gctggcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAggctg                                                                                                                                          >  2:532603/1‑148 (MQ=255)
           ctggcACCACACCGCCTTGCCGCTGGATGGGTTCGACAATCAGCGCACTGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAATTCCTGTGAATAggctgg                                                                                                                                         <  1:2184822/148‑1 (MQ=255)
              gcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGgctggctg                                                                                                                                      >  1:115443/1‑148 (MQ=255)
              gcACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGgctggctg                                                                                                                                      >  1:3214974/1‑148 (MQ=255)
                   acacCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGTCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGATTAGTCTGGCTGccgtc                                                                                                                                 <  2:1802438/148‑4 (MQ=255)
                     acCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCATCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCAccc                                                                                                                               >  1:1316222/1‑148 (MQ=255)
                     acCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCAccc                                                                                                                               >  2:116451/1‑148 (MQ=255)
                     acCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCAccc                                                                                                                               >  1:2000717/1‑148 (MQ=255)
                      ccGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCg                                                                                                                              <  2:2065748/148‑1 (MQ=255)
                       cGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAAAATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGc                                                                                                                             >  2:1297505/1‑148 (MQ=255)
                        gCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCa                                                                                                                            <  1:2247664/148‑1 (MQ=255)
                         ccTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTCTATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCAc                                                                                                                           >  2:2875594/1‑148 (MQ=255)
                         ccTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCAc                                                                                                                           >  1:2463040/1‑148 (MQ=255)
                                 cTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTa                                                                                                                   <  2:921087/148‑1 (MQ=255)
                                  tGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTaa                                                                                                                  <  1:1855063/148‑1 (MQ=255)
                                    gATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAAc                                                                                                                <  1:1782907/148‑1 (MQ=255)
                                      tGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAg                                                                                                              >  1:234705/1‑148 (MQ=255)
                                          ttCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTg                                                                                                          <  2:1357566/148‑1 (MQ=255)
                                           tCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGc                                                                                                         >  2:2167696/1‑148 (MQ=255)
                                               cAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGAcc                                                                                                     >  1:2859005/1‑148 (MQ=255)
                                                aaTCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCCGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCa                                                                                                    >  1:3074389/1‑148 (MQ=255)
                                                  tCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATg                                                                                                  >  2:863307/1‑148 (MQ=255)
                                                   cACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGa                                                                                                 >  1:2147612/1‑148 (MQ=255)
                                                   cACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGa                                                                                                 >  1:3492697/1‑148 (MQ=255)
                                                   cACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGa                                                                                                 >  2:246195/1‑148 (MQ=255)
                                                   cACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAAACCTGTGAATAGGCAGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGa                                                                                                 >  1:842157/1‑148 (MQ=255)
                                                         acaGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCa                                                                                           <  1:2453271/148‑1 (MQ=255)
                                                          caGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAt                                                                                          >  2:255077/1‑148 (MQ=255)
                                                           aGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAtt                                                                                         <  1:996732/148‑1 (MQ=255)
                                                           aGGGAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAtt                                                                                         <  1:3041925/148‑1 (MQ=255)
                                                            ggTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAttt                                                                                        >  2:3256179/1‑148 (MQ=255)
                                                             gTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGAGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCAtttt                                                                                       >  2:290990/1‑148 (MQ=255)
                                                               agagTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTg                                                                                     <  1:660894/148‑1 (MQ=255)
                                                                gagTCGTCAATCAGCGCGCTGTCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGa                                                                                    <  1:1801369/148‑1 (MQ=255)
                                                                   tcgtcAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACg                                                                                 >  2:2746595/1‑148 (MQ=255)
                                                                      tcAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCa                                                                              <  2:1765778/148‑1 (MQ=255)
                                                                         aTCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGa                                                                           >  2:785937/1‑148 (MQ=255)
                                                                           cAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATg                                                                         <  1:3328884/148‑1 (MQ=255)
                                                                            aGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGc                                                                        <  2:2436081/148‑1 (MQ=255)
                                                                            aGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGc                                                                        >  1:22254/1‑148 (MQ=255)
                                                                             gcgcgcTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGAGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGcc                                                                       >  1:1632404/1‑148 (MQ=255)
                                                                             gcgcgcTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGAGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGcc                                                                       >  1:2923616/1‑148 (MQ=255)
                                                                              cgcgcTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCg                                                                      <  1:1105782/148‑1 (MQ=255)
                                                                                 gcTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGctc                                                                   <  2:1481575/148‑1 (MQ=255)
                                                                                  cTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGctct                                                                  <  1:1029140/148‑1 (MQ=255)
                                                                                  cTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGctct                                                                  <  2:3331882/148‑1 (MQ=255)
                                                                                    ggCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTa                                                                >  2:2011425/1‑148 (MQ=255)
                                                                                     tgagagTTTATATCGTTTTATGCAGCATGACTAATATACGACGGCAGTGGAGCAAACTCCTGAGATTAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCTTTTTTGTACGCCACTATTCCGCTCTTAt                                                               <  2:2193257/146‑1 (MQ=255)
                                                                                     gCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTAt                                                               <  1:2072499/148‑1 (MQ=255)
                                                                                         agTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCt                                                           >  1:3388775/1‑148 (MQ=255)
                                                                                         agTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCt                                                           >  1:1340266/1‑148 (MQ=255)
                                                                                          gTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTg                                                          <  1:2431989/148‑1 (MQ=255)
                                                                                              atatCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGt                                                      >  2:635672/1‑148 (MQ=255)
                                                                                               tatCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTa                                                     >  2:2925445/1‑148 (MQ=255)
                                                                                                atCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAg                                                    <  1:2025430/148‑1 (MQ=255)
                                                                                                atCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAg                                                    <  2:794653/148‑1 (MQ=255)
                                                                                                atCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTGATGGCTGCCGTAg                                                    >  1:1935167/1‑148 (MQ=255)
                                                                                                    ttATATGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGt                                                <  2:1979858/148‑1 (MQ=255)
                                                                                                     tataTGCAGCATGACGATTATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTc                                               >  2:2845888/1‑148 (MQ=255)
                                                                                                      atatGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCg                                              >  2:2952214/1‑148 (MQ=255)
                                                                                                      atatGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCg                                              <  1:3256179/148‑1 (MQ=255)
                                                                                                       tatGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGt                                             <  2:1823084/148‑1 (MQ=255)
                                                                                                        atGCAGCATGACGAATATCCGCCGGCAGTGGAGCGAAATCCTGTGAATAGGCTGGCTGCCCGCCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTg                                            >  1:3276387/1‑148 (MQ=255)
                                                                                                        atGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTg                                            >  1:1744623/1‑148 (MQ=255)
                                                                                                        atGCAGCATGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTg                                            <  2:1382372/148‑1 (MQ=255)
                                                                                                           cagcaTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGc                                         <  2:348271/148‑1 (MQ=255)
                                                                                                             gcaTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGAGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGa                                       <  1:531511/148‑1 (MQ=255)
                                                                                                              caTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGaa                                      <  1:421828/148‑1 (MQ=255)
                                                                                                              caTGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGaa                                      >  2:1435841/1‑148 (MQ=255)
                                                                                                                tGACGAATATCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAtt                                    <  2:1850099/148‑1 (MQ=255)
                                                                                                                 gACGAATATCCGCCGGCAGTGGAGCAAAATCCTGAGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAAttt                                   >  2:2401390/1‑148 (MQ=255)
                                                                                                                        atCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGcc                            <  1:255077/148‑1 (MQ=255)
                                                                                                                        atCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGcc                            <  2:3224261/148‑1 (MQ=255)
                                                                                                                         tCCGCCGGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCa                           <  2:165802/148‑1 (MQ=255)
                                                                                                                          ccgccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAg                          >  2:897210/1‑148 (MQ=255)
                                                                                                                          ccgccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAg                          >  2:1153444/1‑148 (MQ=255)
                                                                                                                          ccgccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAg                          >  1:2365305/1‑148 (MQ=255)
                                                                                                                          ccgccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAg                          >  2:1279504/1‑148 (MQ=255)
                                                                                                                          ccgccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAg                          >  1:1950506/1‑148 (MQ=255)
                                                                                                                             ccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGttt                       <  2:2505760/148‑1 (MQ=255)
                                                                                                                             ccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGttt                       <  1:2794094/148‑1 (MQ=255)
                                                                                                                             ccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGttt                       >  1:2620292/1‑148 (MQ=255)
                                                                                                                             ccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGttt                       <  2:2506102/148‑1 (MQ=255)
                                                                                                                             ccgGCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGttt                       <  2:2505755/148‑1 (MQ=255)
                                                                                                                               ggCAGTGGAGCAAAATCCTGAGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTa                     <  2:2130647/148‑1 (MQ=255)
                                                                                                                               ggCAGTGGAGCAAAATCCGGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAACCGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTa                     <  1:1548407/148‑1 (MQ=255)
                                                                                                                                gCAGTGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAg                    <  2:1207175/148‑1 (MQ=255)
                                                                                                                                 cAGTGGAGCAAAATCCTGAGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAgc                   >  1:2616702/1‑148 (MQ=255)
                                                                                                                                    tGGAGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGcgcc                >  2:354585/1‑148 (MQ=255)
                                                                                                                                     ggagcAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGcgccg               <  2:1702155/145‑1 (MQ=255)
                                                                                                                                      gagcAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGcgccgc              >  2:424211/3‑148 (MQ=255)
                                                                                                                                               ccTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGc     >  2:3414345/1‑148 (MQ=255)
                                                                                                                                               ccTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGc     <  2:262267/148‑1 (MQ=255)
                                                                                                                                               ccTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGc     <  1:814310/148‑1 (MQ=255)
                                                                                                                                                  gtgAATAGGCTGGCTGCCCACCCGCACTGACAGTACACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTc  <  2:3022386/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCGTTGCTGGCTGGCACCACACCGCCTTCCCCCTGGATGGGTTCGACAATCACCGCACAGGTAGAGTCGTCAATCAGCGCGCTGGCAGAGTTAATATCGTTATATGCAGCATGACGAATATCCGCCGGCAGTGGCGCAAAATCCTGTGAATAGGCTGGCTGCCCACCCGCACTGACAGTAAACAGCGTGCGACCATGAAACGCATTTTTGAACGCCACGATGCCGCTCTTATGGCTGCCGTAGCGGTCGTGAGCGAATTTACGCGCCAGTTTTAGCGCCGCTTCGTTGGCTTC  >  CP009273/1825628‑1825921

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