Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 3,931,996 A→G 5.0% D118G (GAT→GGT)  CJ019_03714 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1223,931,9960AG 5.0% 262.9 / 3.7 81D118G (GAT→GGT) CJ019_03714hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (16/60);  new base G (2/2);  total (18/63)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

CAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAGGGACCACACACAGATATCATGATCCGAGAGTCGGCACGTGATGCTTTCATTTCGGCTAGAAGC  >  ACB122/3931862‑3932126
                                                                                                                                      |                                                                                                                                  
cAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTtgatga                                                                                                                                    >  3:239640/1‑135 (MQ=255)
       gttgATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACa                                                                                                                               <  6:46196/133‑1 (MQ=255)
              cAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTAtt                                                                                                                     <  7:131208/136‑1 (MQ=255)
               aaGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTAtt                                                                                                                     <  3:232730/135‑1 (MQ=255)
               aaGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTAtt                                                                                                                     <  2:217711/135‑1 (MQ=255)
                   ggCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTaaa                                                                                                                  >  5:225931/1‑134 (MQ=255)
                   ggCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTaaa                                                                                                                  >  4:341184/1‑134 (MQ=255)
                     ccATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAaacaa                                                                                                               <  2:228358/135‑1 (MQ=255)
                     ccATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAaacaa                                                                                                               <  3:341184/135‑1 (MQ=255)
                     ccATTCCGGCGTTCTGACACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAaacaa                                                                                                               <  5:3261/135‑1 (MQ=255)
                           cgtcGTTCTTCAACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATTATTACATAAAAATATTAAACAACAAAAAg                                                                                                        <  8:189556/133‑1 (MQ=255)
                           cgtcGTTCTGCCACTTCAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTAATGATTACATAAAACTATTAAACAACAAAAAg                                                                                                        <  8:253722/133‑1 (MQ=255)
                           cGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGTCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAg                                                                                                        <  7:17304/136‑1 (MQ=255)
                           cGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAg                                                                                                        <  8:67541/136‑1 (MQ=255)
                           cGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAg                                                                                                        <  7:134741/136‑1 (MQ=255)
                            ggCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAg                                                                                                        <  1:289073/135‑1 (MQ=255)
                            ggCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAg                                                                                                        <  5:85303/135‑1 (MQ=255)
                            ggCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAg                                                                                                        <  5:192489/135‑1 (MQ=255)
                            ggCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAg                                                                                                        <  3:16732/135‑1 (MQ=255)
                            ggCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAg                                                                                                        <  3:185626/135‑1 (MQ=255)
                            ggCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAg                                                                                                        <  2:234268/135‑1 (MQ=255)
                            ggCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATAAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAg                                                                                                        <  3:228264/135‑1 (MQ=255)
                            ggCGTGCTGCCACTAAAAATTCAGGCCGCGCAAAAATTACCGCCTTCGGCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCAAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAg                                                                                                        <  1:395779/135‑1 (MQ=255)
                              cGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCACAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCa                                                                                                     <  7:33207/136‑1 (MQ=255)
                              cGTTCTACCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCa                                                                                                     <  8:390829/136‑1 (MQ=255)
                               gTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCa                                                                                                     <  1:81780/135‑1 (MQ=255)
                               gTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCACAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCa                                                                                                     <  6:384916/135‑1 (MQ=255)
                               gTTCTGCCACTAAAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCATTTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCa                                                                                                     <  5:327273/135‑1 (MQ=255)
                                 tCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGTGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCaaa                                                                                                   <  3:308772/135‑1 (MQ=255)
                                     ccACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAAcc                                                                                                 <  8:350015/133‑1 (MQ=255)
                                     ccACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCt                                                                                                <  8:174888/134‑1 (MQ=255)
                                      cACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCt                                                                                                <  6:137575/133‑1 (MQ=255)
                                      cACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCt                                                                                                <  4:239640/133‑1 (MQ=255)
                                      cACTACAAATTCTGGCCGCCCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGAATCCAACATGCGGGCATCCCAATTAATGTTTGATTATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGt                                                                                              <  4:12753/135‑1 (MQ=255)
                                      cACTACAAATTCTGGCAGCGCAAAAATTAACGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGCTCATACCAAATAATGTATGATGATTACATAAAACTATTAAACAAAAAAAAGCCACACCTGt                                                                                              >  5:48871/1‑135 (MQ=255)
                                      cACTACAAATTCTGGACGCGAAAAAATTACAGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATTATTACATAAAACTATTAAACAACAAAAAGCCAAACCt                                                                                                <  2:262388/133‑1 (MQ=255)
                                             aaTTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGTAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTATTGTTTGATGATTACATAAAAATATTAAACAACAAAAAGCCAAACCTGTATACCAt                                                                                       <  2:166855/135‑1 (MQ=255)
                                             aaTTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCAt                                                                                       <  1:327799/135‑1 (MQ=255)
                                             aaTTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCAt                                                                                       <  5:398526/135‑1 (MQ=255)
                                             aaTTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCAt                                                                                       <  6:154559/135‑1 (MQ=255)
                                                     ccGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGActc                                                                               <  3:351909/135‑1 (MQ=255)
                                                      cgcgCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCAt                                                                                       >  8:29307/1‑126 (MQ=255)
                                                      cgcgCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCAt                                                                                       <  7:29307/126‑1 (MQ=255)
                                                      cgcgCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACtctc                                                                             <  7:116317/136‑1 (MQ=255)
                                                      cgcgCAAAAATTACAGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATACAACCAGCGGTCATCCCAAATAATCTATGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTAAACCATGGAGActct                                                                              >  1:87893/1‑135 (MQ=255)
                                                       gcgcAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATTTTTTATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACtctc                                                                             <  4:245220/135‑1 (MQ=255)
                                                       gcgcAAAAATTACCGCCTTCGTCAAGGACAGTAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACtctc                                                                             <  5:156001/135‑1 (MQ=255)
                                                         gcAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCga                                                                           <  6:225931/135‑1 (MQ=255)
                                                             aaaTTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGata                                                                       >  5:371057/1‑135 (MQ=255)
                                                                ttACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTc                                                                    <  5:371705/135‑1 (MQ=255)
                                                                    cGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAgg                                                                <  4:210608/135‑1 (MQ=255)
                                                                        ttCGTCAAGGCCAGGATATTTAATAGCATCCACCCTGCTGTCATCCCAATTAATGTTTGATGATTACTTAAAACTATTAAACCACAAAATGCCATACATGTATACCATGGAGACTCTCGAGATATTCAAGGGAcc                                                            <  6:60180/135‑1 (MQ=255)
                                                                        ttCGTCAAGGCCAGGAAATTTAAGATCATCCTACCTGCGCTGATCCCAATTAATGTTTTATGTTTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATAGCATGGAGACTCTCGAGATATTCAAGGGAcc                                                            <  6:6910/135‑1 (MQ=255)
                                                                        ttCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAGGGAcc                                                            <  5:262638/135‑1 (MQ=255)
                                                                        ttCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAGGGAcc                                                            <  5:332492/135‑1 (MQ=255)
                                                                            tCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAGGGACcacaca                                                       <  8:193215/136‑1 (MQ=255)
                                                                             cAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGCTGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAGGGACcacaca                                                       <  4:43582/135‑1 (MQ=255)
                                                                              aaGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTAATCCCAATTTATGTTTGATGATTAAATAAAACTTTTAAAAAACAAAAAGCCAAACCTGTATACGATGGAGACTCTCGAGATATTCGAGGGACcactca                                                       <  3:38004/134‑1 (MQ=255)
                                                                                   cAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAAGCTGTATACCATGGAGACTCTCGAGATATTCAAGGGACCACACACAGata                                                 <  5:74868/135‑1 (MQ=255)
                                                                                   cAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAGGGACCACACACAGata                                                 <  2:43718/135‑1 (MQ=255)
                                                                                        aaTTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAg                                                                                                        <  8:2159/75‑1 (MQ=255)
                                                                                        aaTTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAg                                                                                                        >  7:2159/1‑75 (MQ=255)
                                                                                             aaGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAg                                                                                                        >  6:334482/1‑70 (MQ=255)
                                                                                             aaGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAg                                                                                                        <  5:334482/70‑1 (MQ=255)
                                                                                                  cATCAAAACTGCGTTCAACCCAATTTATGTTTGATAATAACATAAAACTATTAAAAAACAAAAAGCAAAACCTGTATACCATTGAGACTCTCGAGATATTCAAGGGACCACACACAGATATCATGATCCGTGAGt                                  <  2:237009/135‑1 (MQ=255)
                                                                                                      cAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAACACGCCAAACCAGTATACCATGGAGACTCTCGAGATATTCAAGGGACCACACACAGATATCATGATCCGAGAGTCGGc                              <  8:362275/135‑1 (MQ=255)
                                                                                                       aaCCTGCGGTCATCCCAATTAATGTTTGATGGTTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCAt                                                                                       <  5:136090/77‑1 (MQ=255)
                                                                                                       aaCCTGCGGTCATCCCAATTAATGTTTGATGGTTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCAt                                                                                       <  5:136096/77‑1 (MQ=255)
                                                                                                       aaCCTGCGGTCATCCCAATTAATGTTTGATGGTTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCAt                                                                                       >  6:136090/1‑77 (MQ=255)
                                                                                                       aaCCTGCGGTCATCCCAATTAATGTTTGATGGTTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCAt                                                                                       >  6:136096/1‑77 (MQ=255)
                                                                                                       aaCCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAGGGACCACACACAGATATCATGATCCGAGAGTCGGCa                             >  4:223871/1‑135 (MQ=255)
                                                                                                       aaCCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAGGGACCACACACAGATATAATGATCCGAGAGTCGGCAc                            >  8:389717/1‑136 (MQ=255)
                                                                                                              ggTCATCGCAATTTATGTTTGATGATTACATAAAACTATTTATATACAAAAAGCCAAAACTTTTTACCATGGAGACTCTCGTGATATTCAAGGGACCACACACTGATTTCATGATCCTAGAGTCGGCACGTGAt                       <  2:87893/134‑1 (MQ=255)
                                                                                                                tCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAGGGACCACACACAGATATCATGATCCGAGAGTCGGCACGTGATGCt                    >  5:379559/1‑135 (MQ=255)
                                                                                                                     ccAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGatat                                                                      <  2:34178/80‑1 (MQ=255)
                                                                                                                     ccAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGatat                                                                      >  1:34178/1‑80 (MQ=255)
                                                                                                                      cAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAGGGACCACACACAGATATCATGATCCGAGAGTCGGCACGTGATGCTTTCAttt             >  7:236476/1‑136 (MQ=255)
                                                                                                                      cAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAGGGACCACACACAGATATCATGATCCGAGAGTCGGCACGTGATGCTTTCAtt              >  2:340137/1‑135 (MQ=255)
                                                                                                                          tAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAgg                                                                >  7:98173/1‑81 (MQ=255)
                                                                                                                          tAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAgg                                                                <  8:98173/81‑1 (MQ=255)
                                                                                                                            aTGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAGGGACCACACACAGATATCATGATCCGAGAGTCGGCACGTGATGCTTTCATTTCgg          >  4:189276/1‑133 (MQ=255)
                                                                                                                                ttGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAGGGACCACACACAGATATCATGATCCGAGAGTCGGCACGTGATGCTTTCATTTCGGCTAGaa    <  3:189276/135‑1 (MQ=255)
                                                                                                                                    tgatTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAGGGACCACACACAGATATCATGATCCGAGAGTCGGCACGTGATGCTTTCATTTCGGCTAGAAGc  >  1:82803/1‑133 (MQ=255)
                                                                                                                                      |                                                                                                                                  
CAACTTTGTTGATCCAAGTGGCCATTCCGGCGTTCTGCCACTACAAATTCTGGCCGCGCAAAAATTACCGCCTTCGTCAAGGCCAGGAAATTTAAGAGCATCCAACCTGCGGTCATCCCAATTAATGTTTGATGATTACATAAAACTATTAAACAACAAAAAGCCAAACCTGTATACCATGGAGACTCTCGAGATATTCAAGGGACCACACACAGATATCATGATCCGAGAGTCGGCACGTGATGCTTTCATTTCGGCTAGAAGC  >  ACB122/3931862‑3932126

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