Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,563,276 G→T 100% H491N (CAC→AAC)  gadC ← glutamate:gamma‑aminobutyric acid antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,563,2760GT98.4% 777.4 / 2.9 243H491N (CAC→AAC) gadCglutamate:gamma‑aminobutyric acid antiporter
Reads supporting (aligned to +/- strand):  ref base G (4/0);  new base T (132/107);  total (136/107)
Fisher's exact test for biased strand distribution p-value = 1.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 243 total aligned reads displayed.

CGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATA  >  CP009273/1563131‑1563418
                                                                                                                                                 |                                                                                                                                              
  ttAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCt                                                                                                                                            >  1:3092014/1‑148 (MQ=255)
  ttAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCt                                                                                                                                            >  2:376077/1‑148 (MQ=255)
   tAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCtt                                                                                                                                           >  1:2039142/1‑148 (MQ=255)
   tAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCtt                                                                                                                                           >  1:1985306/1‑148 (MQ=255)
   tAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCtt                                                                                                                                           <  2:2384565/148‑1 (MQ=255)
          aTACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGc                                                                                                                                    >  1:172416/1‑148 (MQ=255)
               aCAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCt                                                                                                                               >  2:880718/1‑148 (MQ=255)
                cAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTg                                                                                                                              >  2:393537/1‑148 (MQ=255)
                cAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTg                                                                                                                              >  2:3398871/1‑148 (MQ=255)
                 agcaagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGa                                                                                                                             >  1:328478/1‑148 (MQ=255)
                 agcaagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGa                                                                                                                             >  1:1652800/1‑148 (MQ=255)
                       caCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTg                                                                                                                       <  1:281666/148‑1 (MQ=255)
                           aaaaaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTg                                                                                                                   >  2:2968657/1‑148 (MQ=255)
                              aaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCt                                                                                                                >  1:1118613/1‑148 (MQ=255)
                                 ggAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCa                                                                                                             >  1:90746/1‑148 (MQ=255)
                                  gAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCag                                                                                                            <  2:611061/148‑1 (MQ=255)
                                     cGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGt                                                                                                         >  2:2250277/1‑148 (MQ=255)
                                     cGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGt                                                                                                         <  1:3099132/148‑1 (MQ=255)
                                      gATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTg                                                                                                        <  2:3076636/148‑1 (MQ=255)
                                       aTGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGa                                                                                                       >  2:1357541/1‑148 (MQ=255)
                                        tGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGAc                                                                                                      >  1:2634537/1‑148 (MQ=255)
                                         gAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACg                                                                                                     >  1:650530/1‑148 (MQ=255)
                                          aaTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGACGAAGTTACCTTTTGGTGCGTTATGACTGTTGATTGGCTCCCGAGTGACGc                                                                                                    >  2:1313161/1‑148 (MQ=255)
                                          aaTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGc                                                                                                    >  2:218012/1‑148 (MQ=255)
                                           aTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGcc                                                                                                   <  2:1747395/148‑1 (MQ=255)
                                           aTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGcc                                                                                                   >  1:486906/1‑148 (MQ=255)
                                            ttATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCg                                                                                                  <  2:1469050/148‑1 (MQ=255)
                                            ttATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCg                                                                                                  >  2:3054880/1‑148 (MQ=255)
                                             tATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCgg                                                                                                 >  2:2770187/1‑148 (MQ=255)
                                             tATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCgg                                                                                                 >  2:1639446/1‑148 (MQ=255)
                                                cGCTCCCTTGTCTTATAACCATTCATACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAt                                                                                              >  1:84317/1‑148 (MQ=255)
                                                 gCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtt                                                                                             >  1:1841146/1‑148 (MQ=255)
                                                    cccTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtttgc                                                                                          <  2:2912071/148‑1 (MQ=255)
                                                         gTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtttgctttgc                                                                                     <  1:2606667/148‑1 (MQ=255)
                                                           cTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCt                                                                                   <  1:2102752/148‑1 (MQ=255)
                                                              ataACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTa                                                                                <  1:2929859/148‑1 (MQ=255)
                                                              ataACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTa                                                                                <  1:1604706/148‑1 (MQ=255)
                                                              ataACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTa                                                                                <  1:561942/148‑1 (MQ=255)
                                                                 aCCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGa                                                                             <  1:1353075/148‑1 (MQ=255)
                                                                  ccATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGAt                                                                            <  2:2654097/148‑1 (MQ=255)
                                                                   cATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATc                                                                           <  1:1150185/148‑1 (MQ=255)
                                                                    aTTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCa                                                                          >  2:195986/1‑148 (MQ=255)
                                                                    aTTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCa                                                                          >  2:1663355/1‑148 (MQ=255)
                                                                     ttCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCAt                                                                         >  1:3537625/1‑148 (MQ=255)
                                                                        aGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGaa                                                                      >  2:2315728/1‑148 (MQ=255)
                                                                          aCATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACa                                                                    >  2:1388249/1‑148 (MQ=255)
                                                                              ggTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCAt                                                                <  2:1177187/148‑1 (MQ=255)
                                                                               gTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATa                                                               <  2:3377155/148‑1 (MQ=255)
                                                                                ttAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATag                                                              <  2:1134249/148‑1 (MQ=255)
                                                                                 tAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATaga                                                             <  1:2250277/148‑1 (MQ=255)
                                                                                  aGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATagag                                                            <  2:2976395/148‑1 (MQ=255)
                                                                                    tgtTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGaa                                                          <  2:3449680/148‑1 (MQ=255)
                                                                                    tgtTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGaa                                                          >  1:2220398/1‑148 (MQ=255)
                                                                                    tgtTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGaa                                                          >  2:232197/1‑148 (MQ=255)
                                                                                      tttCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATa                                                        >  2:2451323/1‑148 (MQ=255)
                                                                                      tttCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATa                                                        >  1:3109871/1‑148 (MQ=255)
                                                                                       ttCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATaa                                                       >  2:3150985/1‑148 (MQ=255)
                                                                                       ttCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATaa                                                       >  1:813051/1‑148 (MQ=255)
                                                                                        tCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATaaa                                                      <  1:3172442/148‑1 (MQ=255)
                                                                                         cTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAag                                                     >  2:2072325/1‑148 (MQ=255)
                                                                                          ttGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagg                                                    <  2:2347373/148‑1 (MQ=255)
                                                                                          ttGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagg                                                    >  1:3224384/1‑148 (MQ=255)
                                                                                          ttGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagg                                                    <  2:562114/148‑1 (MQ=255)
                                                                                           tGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAaggg                                                   <  2:1935880/148‑1 (MQ=255)
                                                                                                ttcatCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggcaggg                                              <  1:988436/148‑1 (MQ=255)
                                                                                                ttcatCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggcaggg                                              >  1:976297/1‑148 (MQ=255)
                                                                                                 tcatcaCAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggcagggc                                             >  2:794269/1‑148 (MQ=255)
                                                                                                    tcaCTATTTAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAg                                          <  2:3227174/148‑1 (MQ=255)
                                                                                                     cacaATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGt                                         >  1:2826752/1‑148 (MQ=255)
                                                                                                     cacaATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGt                                         >  2:846960/1‑148 (MQ=255)
                                                                                                     cacaATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGt                                         >  1:1007119/1‑148 (MQ=255)
                                                                                                          tataGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTaccac                                    >  1:3313998/1‑148 (MQ=255)
                                                                                                             aGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAg                                 >  1:2996125/1‑148 (MQ=255)
                                                                                                              gtgtGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTGTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAgg                                <  2:3323425/148‑1 (MQ=255)
                                                                                                               tgtgGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGa                               >  1:423932/1‑148 (MQ=255)
                                                                                                                 tggtgAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGaaa                             <  2:1974898/148‑1 (MQ=255)
                                                                                                                    tgAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTa                          <  2:2093528/148‑1 (MQ=255)
                                                                                                                    tgAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTa                          <  1:3336911/148‑1 (MQ=255)
                                                                                                                       aCGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTaaca                       >  1:2392955/1‑148 (MQ=255)
                                                                                                                        cGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTaacaa                      <  2:3344502/148‑1 (MQ=255)
                                                                                                                        cGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTaacaa                      <  2:1583519/148‑1 (MQ=255)
                                                                                                                          tGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAAcc                    <  1:1339738/148‑1 (MQ=255)
                                                                                                                           gCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCa                   <  2:1910092/148‑1 (MQ=255)
                                                                                                                           gCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCa                   <  1:1922361/148‑1 (MQ=255)
                                                                                                                            cACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAg                  <  2:328921/148‑1 (MQ=255)
                                                                                                                            cACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAg                  <  1:2058028/148‑1 (MQ=255)
                                                                                                                                cgGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAAt              <  1:406561/148‑1 (MQ=255)
                                                                                                                                  ggTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTc            >  2:1772263/1‑148 (MQ=255)
                                                                                                                                   gTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCa           <  2:90746/148‑1 (MQ=255)
                                                                                                                                   gTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCa           <  2:1891314/148‑1 (MQ=255)
                                                                                                                                   gTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCa           <  1:1808374/148‑1 (MQ=255)
                                                                                                                                   gTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCa           >  1:1926238/1‑148 (MQ=255)
                                                                                                                                    tGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGTTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCaa          <  2:298276/148‑1 (MQ=255)
                                                                                                                                     gCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAac         <  1:2968657/148‑1 (MQ=255)
                                                                                                                                       aGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacat       >  2:1266731/1‑148 (MQ=255)
                                                                                                                                         gaagaaGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacatac     <  2:2038586/148‑1 (MQ=255)
                                                                                                                                         gaagaaGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacatac     >  2:1829140/1‑148 (MQ=255)
                                                                                                                                         gaagaaGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacatac     <  1:2931544/148‑1 (MQ=255)
                                                                                                                                            gaagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAAcatacata  >  1:2133550/1‑148 (MQ=255)
                                                                                                                                            gaagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAAcatacata  >  2:2506435/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                              
CGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATA  >  CP009273/1563131‑1563418

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: