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,2760GT95.9% 969.2 / 24.7 279H491N (CAC→AAC) gadCglutamate:gamma‑aminobutyric acid antiporter
Reads supporting (aligned to +/- strand):  ref base G (8/4);  new base T (156/111);  total (164/115)
Fisher's exact test for biased strand distribution p-value = 7.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 279 total aligned reads displayed.

CCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCGG  >  CP009273/1563130‑1563422
                                                                                                                                                  |                                                                                                                                                  
 cGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCt                                                                                                                                                <  1:4693933/150‑1 (MQ=255)
    tAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCtttt                                                                                                                                             <  1:1880620/150‑1 (MQ=255)
        tCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTg                                                                                                                                         >  2:3002201/1‑150 (MQ=255)
           aTACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGt                                                                                                                                      <  1:1970971/150‑1 (MQ=255)
                aCAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCGTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAAGCTAGTGTGGTGAACGTGCACGCGGGTGCAGGACGAAGTTACCTTTTGGTGCGGTCTGa                                                                                                                                 >  2:356801/1‑150 (MQ=255)
                  agcaagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACt                                                                                                                               >  2:3489936/1‑150 (MQ=255)
                   gcaagcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTg                                                                                                                              <  2:3639727/150‑1 (MQ=255)
                      agcaCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTg                                                                                                                           <  2:2969792/150‑1 (MQ=255)
                        caCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGAt                                                                                                                         >  2:2493988/1‑150 (MQ=255)
                        caCGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGAt                                                                                                                         >  1:1469408/1‑150 (MQ=255)
                           gAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTgg                                                                                                                      <  1:1033502/150‑1 (MQ=255)
                           gAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTgg                                                                                                                      >  1:328889/1‑150 (MQ=255)
                           gAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTgg                                                                                                                      >  1:3169135/1‑150 (MQ=255)
                             aaaaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCt                                                                                                                    >  1:946916/1‑150 (MQ=255)
                             aaaaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCt                                                                                                                    >  1:571148/1‑150 (MQ=255)
                             aaaaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCt                                                                                                                    >  2:2821061/1‑150 (MQ=255)
                              aaaGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTc                                                                                                                   <  1:709049/150‑1 (MQ=255)
                                aGGGAGCGCTGAATTATCGCTCCCTTGTCGTATACCCATTCCGACATGGTTAGTGTGTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCa                                                                                                                 <  2:166256/150‑1 (MQ=255)
                                aGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGACGTTACCTTTTGGTGCTGTCTGACTGTTGATTGGCTCCa                                                                                                                 >  2:533524/1‑150 (MQ=255)
                                aGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCa                                                                                                                 >  2:3858415/1‑150 (MQ=255)
                                aGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCa                                                                                                                 >  1:3450951/1‑150 (MQ=255)
                                  ggAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCaga                                                                                                               >  1:2612648/1‑150 (MQ=255)
                                    aGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGt                                                                                                             >  2:327348/1‑150 (MQ=255)
                                         tGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGc                                                                                                        >  2:3907065/1‑150 (MQ=255)
                                            aTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCgg                                                                                                     <  2:3964902/150‑1 (MQ=255)
                                            aTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCgg                                                                                                     <  2:1578384/150‑1 (MQ=255)
                                            aTTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCgg                                                                                                     >  2:4498438/1‑150 (MQ=255)
                                                 cGCTCCCTTGTCTTATAACCAGTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAttt                                                                                                <  2:1608923/150‑1 (MQ=255)
                                                  gCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAAGATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtttg                                                                                               >  2:4406971/1‑150 (MQ=255)
                                                   cTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtttgc                                                                                              <  1:3831069/150‑1 (MQ=255)
                                                    tCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtttgct                                                                                             >  2:2885203/1‑150 (MQ=255)
                                                    tCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtttgct                                                                                             >  2:4078002/1‑150 (MQ=255)
                                                      ccTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTAtttgcttt                                                                                           <  1:4036378/150‑1 (MQ=255)
                                                           tCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCtt                                                                                      >  2:1505405/1‑150 (MQ=255)
                                                            cTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCttt                                                                                     <  2:4124138/150‑1 (MQ=255)
                                                            cTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATAGAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCttt                                                                                     <  1:510936/150‑1 (MQ=255)
                                                             ttATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTa                                                                                    >  1:802634/1‑150 (MQ=255)
                                                              tataACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTAc                                                                                   >  1:3527176/1‑150 (MQ=255)
                                                                taACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGa                                                                                 >  1:841424/1‑150 (MQ=255)
                                                                  aCCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATc                                                                               >  1:256261/1‑150 (MQ=255)
                                                                     aTTCAGACATGGTGAGTGTTTCTTGTCAGTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATg                                                                            <  1:1202219/150‑1 (MQ=255)
                                                                      ttCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGa                                                                           >  1:558084/1‑150 (MQ=255)
                                                                      ttCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGa                                                                           >  1:3355543/1‑150 (MQ=255)
                                                                         aGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACa                                                                        <  1:1775849/150‑1 (MQ=255)
                                                                         aGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACa                                                                        <  1:1775466/150‑1 (MQ=255)
                                                                         aGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACa                                                                        >  1:1632794/1‑150 (MQ=255)
                                                                           aCATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGc                                                                      >  1:2721840/1‑150 (MQ=255)
                                                                               ggTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATag                                                                  >  2:3540971/1‑150 (MQ=255)
                                                                                 ttAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATagag                                                                >  2:1420330/1‑150 (MQ=255)
                                                                                     tgtTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATa                                                            >  2:4436410/1‑150 (MQ=255)
                                                                                     tgtTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATa                                                            <  1:1079409/150‑1 (MQ=255)
                                                                                     tgtTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATa                                                            <  1:4498438/150‑1 (MQ=255)
                                                                                     tgtTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATa                                                            >  1:2071440/1‑150 (MQ=255)
                                                                                      gtTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATaa                                                           <  1:4575472/150‑1 (MQ=255)
                                                                                       tttCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATaaa                                                          >  2:4170574/1‑150 (MQ=255)
                                                                                          cTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAaggg                                                       >  1:3642855/1‑150 (MQ=255)
                                                                                          cTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAaggg                                                       >  2:182945/1‑150 (MQ=255)
                                                                                          cTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCGTTGCCTTGACGATCATGAACAGCATCGAGACTAAaggg                                                       >  2:3189403/1‑150 (MQ=255)
                                                                                           ttGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggc                                                      >  2:1100788/1‑150 (MQ=255)
                                                                                            tGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggca                                                     >  2:3280888/1‑150 (MQ=255)
                                                                                            tGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggca                                                     <  1:3899734/150‑1 (MQ=255)
                                                                                              tcattcatCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggcggg                                                   >  2:2325483/1‑150 (MQ=255)
                                                                                              tcattcatCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAagggcagg                                                   >  1:231124/1‑150 (MQ=255)
                                                                                                 ttcatCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAgggcagggca                                                <  2:148739/150‑1 (MQ=255)
                                                                                                 ttcatCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAgggcagggca                                                <  1:127381/150‑1 (MQ=255)
                                                                                                  tcatcaCAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCaa                                               >  2:2659162/1‑150 (MQ=255)
                                                                                                   catcaCAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAg                                              >  1:3781911/1‑150 (MQ=255)
                                                                                                       acaATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTacc                                          >  1:1305102/1‑150 (MQ=255)
                                                                                                       acaATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTacc                                          >  2:840906/1‑150 (MQ=255)
                                                                                                         aaTATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTaccac                                        >  2:461903/1‑150 (MQ=255)
                                                                                                           tataGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTAccacca                                      >  2:4651012/1‑150 (MQ=255)
                                                                                                            atcgtgtGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAg                                     <  1:2176265/147‑1 (MQ=255)
                                                                                                            ataGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAg                                     <  2:2847559/150‑1 (MQ=255)
                                                                                                            ataGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAg                                     <  1:2804483/150‑1 (MQ=255)
                                                                                                            ataGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAg                                     >  2:3349437/1‑150 (MQ=255)
                                                                                                                tgtgGTGAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGaaa                                 <  1:2997285/150‑1 (MQ=255)
                                                                                                                    gtgAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTaa                             <  1:3002201/150‑1 (MQ=255)
                                                                                                                     tgAACGTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTaac                            >  1:3104192/1‑150 (MQ=255)
                                                                                                                          gTGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATCTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCa                       <  1:540951/150‑1 (MQ=255)
                                                                                                                           tGCACGCGGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAg                      >  1:3236860/1‑150 (MQ=255)
                                                                                                                               cgcgGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAAt                  >  1:1524169/1‑150 (MQ=255)
                                                                                                                               cgcgGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAAt                  <  2:4375467/150‑1 (MQ=255)
                                                                                                                               cgcgGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAAt                  >  1:2361091/1‑150 (MQ=255)
                                                                                                                                gcgGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAAtt                 <  2:4226145/150‑1 (MQ=255)
                                                                                                                                 cgGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTCCCACCAGGAAACTAACAACCAGTAATTc                >  1:226322/1‑150 (MQ=255)
                                                                                                                                 cgGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTc                <  1:3151130/150‑1 (MQ=255)
                                                                                                                                 cgGGTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTc                <  1:4334322/150‑1 (MQ=255)
                                                                                                                                  gggTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCa               <  2:2612648/150‑1 (MQ=255)
                                                                                                                                   ggTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCaa              <  2:1837822/150‑1 (MQ=255)
                                                                                                                                   ggTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCaa              >  2:1906728/1‑150 (MQ=255)
                                                                                                                                    gTGCAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAac             <  2:2256567/150‑1 (MQ=255)
                                                                                                                                       cAGGAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacata          >  1:3295799/1‑150 (MQ=255)
                                                                                                                                         ggAAGAAGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacataca        >  2:376959/1‑150 (MQ=255)
                                                                                                                                          gaagaaGTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAacatacat       >  2:2490791/1‑150 (MQ=255)
                                                                                                                                           aagaagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAAcatacata      <  2:1786719/150‑1 (MQ=255)
                                                                                                                                           aagaagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAAcatacata      >  2:2109379/1‑150 (MQ=255)
                                                                                                                                             gaagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATc    >  1:4343036/1‑150 (MQ=255)
                                                                                                                                              aagTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCg   <  2:4646108/150‑1 (MQ=255)
                                                                                                                                               agTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCgg  <  1:1357321/150‑1 (MQ=255)
                                                                                                                                               agTTACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCgg  >  2:1308628/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATTATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGTATTTGCTTTGCCTTTACGATCATGAACAGCATAGAGAATAAAGGGCAGGGCAAGTACCACCAGGAAACTAACAACCAGTAATTCAACATACATATCGG  >  CP009273/1563130‑1563422

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: