Predicted mutation
evidence seq id position mutation freq annotation gene description
RA ACB122 1,010,006 G→A 5.9% Y45Y (TAC→TAT CJ019_00951 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*ACB1221,010,0060GA 5.9% 230.4 / 2.5 72Y45Y (TAC→TATCJ019_00951hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (26/39);  new base A (2/2);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.14e-01

TGAATGCGGTGCTGCCATCGACGCGGATGCGGTCGCTGTGGACGATGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCGTCGGGGAGGGTGGTGGCGATCA  >  ACB122/1009873‑1010139
                                                                                                                                     |                                                                                                                                     
tGAATGCGGTGCTGCCATCGACGCGGATGCGGTCGCTGTGGACGATGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATTGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGt                                                                                                                                      <  8:99190/135‑1 (MQ=255)
 gAATGCGGTGCTGCCATCGACGCGGATGCGGTCGCTGTGGACGATGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTAGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTa                                                                                                                                     >  2:45174/1‑135 (MQ=255)
     gCGGTGCTGCCATCGACGCGGATGCGGTCGCTGTGGACGATGTCCTCGCCGGCTTCATGGGGCAGGCCGGCGCGTACGTACTGGCGGCGTTACTGATGGCCGCGCCCGCGGTTGTGAAAGAAGTCGTGGTACAgg                                                                                                                                 >  2:354799/1‑135 (MQ=255)
     gCGGTGCTGCCATCGACGCGGATGCGGTCGCTGTGGACGATGTCCTCGCCGGCTGGATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTgg                                                                                                                                 >  2:87910/1‑135 (MQ=255)
      cGGTGCTGCCATCGACGCGGATGCGGTCGCTGTGGACGATGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGa                                                                                                                                >  3:278769/1‑135 (MQ=255)
        gTGCTGCCATCGACGCGGATGCGGTCGCTGTGGACGATGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGata                                                                                                                              <  8:365096/135‑1 (MQ=255)
        gTGCTGCCATCGACGCGGATGCGGTCGATGTGGACGATGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTGGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGata                                                                                                                              <  5:327360/135‑1 (MQ=255)
           ctgcCATCGACGCGGATGCGGTCGCTGTGGACGATGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCg                                                                                                                           <  3:300507/135‑1 (MQ=255)
           ctgcCATCGACGCGGATGCGGTCGCTGTGGACGAGGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCg                                                                                                                           <  3:209952/135‑1 (MQ=255)
                 tCGACGCTGATGCGGTCGCTGTGGACGATGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGggtgg                                                                                                                     <  6:84271/135‑1 (MQ=255)
                 tCGACGCGGATGGGGGCGCTGTGGACGATGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTGCTGAAAGAAGTCGTGGTACTGGATATCGGggtgg                                                                                                                     <  6:89694/135‑1 (MQ=255)
                 tCGACGCGGATGCGGTCGCTGTGGACGATGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGCGCACGCGGGTCTGAAAGCAGTCGGGGTACTGGATATCGGggtgg                                                                                                                     >  6:260925/1‑135 (MQ=255)
                    aCGCGGATGCGGTCGCTGTGGACGATGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTa                                                                                                                   <  3:345261/134‑1 (MQ=255)
                      gcgGATGCGGTCGCTGTGGACGATGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGGAGTTCCTGATGGCCGAGCCCGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTagag                                                                                                                >  7:6026/1‑135 (MQ=255)
                      gcgGATGCGGTCGCTGTGGACGATGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGCGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTagag                                                                                                                >  4:190651/1‑135 (MQ=255)
                                           gATGTCCTCGCGGGCCTCATGGGGCAGGCAGGCGCGTACGGAGGCGCGCAGTTCCTGATGGCCGAGCCCGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGGGGTAGAGCGCCATGACACCGTCGAGGTc                                                                                           <  2:256358/135‑1 (MQ=255)
                                           gATGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTc                                                                                           <  3:170677/135‑1 (MQ=255)
                                           gATGTCCGCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTc                                                                                           <  6:237441/135‑1 (MQ=255)
                                             tGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTcgc                                                                                         <  2:65490/135‑1 (MQ=255)
                                             tGTCCTCGCCGGCTGCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTcgc                                                                                         <  6:329912/135‑1 (MQ=255)
                                             ggTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTcgc                                                                                         <  5:316310/134‑1 (MQ=255)
                                                ccTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGAt                                                                                      >  2:318683/1‑135 (MQ=255)
                                                ccTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGAGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGAt                                                                                      >  1:270640/1‑135 (MQ=255)
                                                 cTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTGCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGGACTGGATATCGGGGTGGTTGCGCGCCATGACACCGGCGAGGTCGCGATg                                                                                     >  5:223122/1‑135 (MQ=255)
                                                   cGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCt                                                                                   >  2:68739/1‑135 (MQ=255)
                                                         cTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCcagc                                                                             <  4:32162/135‑1 (MQ=255)
                                                         cTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCcagc                                                                             <  4:278769/135‑1 (MQ=255)
                                                                 ggcaggcaGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAggtgg                                                                     <  8:44093/135‑1 (MQ=255)
                                                                      gcaggcGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAgc                                                                 >  4:133474/1‑134 (MQ=255)
                                                                              gtacgtaGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTcgc                                                                                         >  8:187794/1‑102 (MQ=255)
                                                                              gtacgtaGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTcgc                                                                                         <  7:187794/102‑1 (MQ=255)
                                                                                acgtaGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGAtt                                                      <  4:163273/135‑1 (MQ=255)
                                                                                acgtaGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGAtt                                                      <  3:133474/135‑1 (MQ=255)
                                                                                acgtaGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAACTCGTGCTACTGGATATCCGGGTGGTACAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGCTGGTAGCGCAGCACCAtt                                                      >  5:123683/1‑135 (MQ=255)
                                                                                   taGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATTCTTCCAGCAGCGGTGGTAGCGCAGCACGATTTcg                                                   >  3:214042/1‑135 (MQ=255)
                                                                                   taGTCGCGCAGTTCCTCATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGCCACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTcg                                                   >  7:275263/1‑135 (MQ=255)
                                                                                        gcgcAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGt                                              <  4:251059/135‑1 (MQ=255)
                                                                                         cgcAGTTCCTGATGGCCGATCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGTGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGt                                              <  4:151341/134‑1 (MQ=255)
                                                                                            aGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGCTACTGGATATCGGGGTGGTAGCGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGc                                          >  4:345090/1‑135 (MQ=255)
                                                                                              ttCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCgg                                        <  7:78170/135‑1 (MQ=255)
                                                                                               tCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTTGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCgg                                        <  1:354799/134‑1 (MQ=255)
                                                                                                     tGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGttt                                   <  3:52917/133‑1 (MQ=255)
                                                                                                      ggCCGAGCACGCGGGTCTGAATGAAGTCGTGGTACTGGATATCGGGGGGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTACAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCt                                <  6:342187/135‑1 (MQ=255)
                                                                                                        ccGAGCCCGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCg                              <  5:113336/135‑1 (MQ=255)
                                                                                                        ccGAGCACGCGGTTCTGAAAGACGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCg                              <  5:289449/135‑1 (MQ=255)
                                                                                                           aGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTcgc                                                                                         <  4:15943/73‑1 (MQ=255)
                                                                                                           aGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTcgc                                                                                         >  3:15943/1‑73 (MQ=255)
                                                                                                           aGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCgggg                           >  4:351530/1‑135 (MQ=255)
                                                                                                           aGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCgggg                           >  5:106790/1‑135 (MQ=255)
                                                                                                              aCGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCg                        >  3:330064/1‑135 (MQ=255)
                                                                                                              aCGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCg                        >  2:25166/1‑135 (MQ=255)
                                                                                                                gcgGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCGTc                      <  4:330064/135‑1 (MQ=255)
                                                                                                                gcgGTTCTGAAAGAAGGCGTGTTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCggc                      <  2:190369/135‑3 (MQ=255)
                                                                                                                     tCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGa                                                                                               >  4:355333/1‑57 (MQ=255)
                                                                                                                     tCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGa                                                                                               <  3:355333/57‑1 (MQ=255)
                                                                                                                         aaaGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTcc                                                                                <  2:358451/68‑1 (MQ=255)
                                                                                                                         aaaGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTcc                                                                                <  8:336604/68‑1 (MQ=255)
                                                                                                                         aaaGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTcc                                                                                >  7:336604/1‑68 (MQ=255)
                                                                                                                         aaaGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTcc                                                                                >  1:358451/1‑68 (MQ=255)
                                                                                                                         aaaGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCGTCGGGGAGggt             >  1:126151/1‑135 (MQ=255)
                                                                                                                          aagaagTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTATGCGGGGTTTGCTCGGGGCCGTCGGGGAggggg            <  2:321347/135‑3 (MQ=255)
                                                                                                                           agaagTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCGTCGGGGAGggtgg           <  5:128544/135‑1 (MQ=255)
                                                                                                                               gTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCGTCGGGGAGGGTGGTGGc       <  3:190651/135‑1 (MQ=255)
                                                                                                                                 cGTTGTACTGGATATCGGGGTGGTAGAGCGCCATGACCCCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCGTCGGGGAGGGTGGTGGCGa     <  1:71450/135‑1 (MQ=255)
                                                                                                                                 cGTGGTACTGGATCTATGTGTGTTAGAGCGCCATGACACCTTCGAGGTCGCGATGCTTCCAGCAGAGGTTGTAGCTCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCGTCGGGGAGGGTGGTTGCGa     <  4:119662/135‑1 (MQ=255)
                                                                                                                                 cGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCgg                                        <  1:357015/100‑1 (MQ=255)
                                                                                                                                 cGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCgg                                        >  2:357015/1‑100 (MQ=255)
                                                                                                                                 cGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCGTCGGGGAGGGTGGTGGCGa     >  4:23830/1‑135 (MQ=255)
                                                                                                                                 cGTGATACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCGTCGGGGAGGGTGGTGGc       >  7:269596/1‑133 (MQ=255)
                                                                                                                                 cGTGATACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCGTCGGGGAGGGTGGTGGc       <  6:381099/133‑1 (MQ=255)
                                                                                                                                 cGTGATACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCGTCGGGGAGGGTGGTGGc       <  8:269596/133‑1 (MQ=255)
                                                                                                                                 cGTGATACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCGTCGGGGAGGGTGGTGGc       >  5:381099/1‑133 (MQ=255)
                                                                                                                                    ggTAGTGGATTTCGGGGTGGTAGAGCTCCATGACACAGTCGATGTCGCGATGCTTCCAGCAGTGGTGGTAGCGCAGCACGTTTTCGGCGGTAAGGGGGGTTTGCTCGGGGCCGTCGGGGAGGGTGGTGGCGATCa  <  1:68739/135‑1 (MQ=255)
                                                                                                                                    ggTACTGGATAGCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCGTCGGGGAGGGTGGTGGCGATCa  <  2:283755/135‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                     
TGAATGCGGTGCTGCCATCGACGCGGATGCGGTCGCTGTGGACGATGTCCTCGCCGGCTTCATGGGGCAGGCAGGCGCGTACGTAGTCGCGCAGTTCCTGATGGCCGAGCACGCGGTTCTGAAAGAAGTCGTGGTACTGGATATCGGGGTGGTAGAGCGCCATGACACCGTCGAGGTCGCGATGCTTCCAGCAGAGGTGGTAGCGCAGCACGATTTCGGCGGTAAGCGGGGTTTGCTCGGGGCCGTCGGGGAGGGTGGTGGCGATCA  >  ACB122/1009873‑1010139

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: