Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,371,363 C→T 5.6% E45E (GAG→GAA hisC ← histidinol‑phosphate aminotransferase; tyrosine/phenylalanine aminotransferase (promiscuous)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,371,3630CT 5.6% 242.7 / 2.5 71E45E (GAG→GAAhisChistidinol‑phosphate aminotransferase; tyrosine/phenylalanine aminotransferase (promiscuous)
Reads supporting (aligned to +/- strand):  ref base C (32/35);  new base T (2/2);  total (34/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

AGAACCGTTGCCGAAAATAAGTGATGTTTCGCTGACATTGAGGTGCTTACTGAGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAGGCTGATAAGGCTTCAGCTGTTTTAAATGTTCTTTGATACGCAAA  >  BS168_reference_genome/2371228‑2371498
                                                                                                                                       |                                                                                                                                       
aGAACCGGTGCCTACATTACGTTAGGTTTCGCAGACCTTGAGGTGCTTAGTGAGCCTGGCCGCCAATGCGGCGCTATCCCCGTCCGGGTAATGGTCATGCTGTTGAAACTCATCATGAAGCGCCTCTTTTGCAGCct                                                                                                                                        <  4:80265/137‑1 (MQ=255)
 gAACCGTTGCCGAAAATAAGCGATGTTTCGCTGACATTGAGGGGCTTACTGAGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTg                                                                                                                                     <  3:21677/139‑1 (MQ=255)
   aCCGTTGCCTAAAATAAGTGATGTTTCGCTGACATTGAGGTGCTTACTGAGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGaa                                                                                                                                   <  1:292019/139‑1 (MQ=255)
   aCCGTTGCCGAAAATAAGTGATGTTTCGCTGACATTGAGGTGCTTACTGAGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGaa                                                                                                                                   <  5:335919/139‑1 (MQ=255)
    ccGTTGCCGAAAATAAGTGATGTTTCGCTGACATTGAGGTGCTTACTGCGCCTTGTCCGAAACGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAAc                                                                                                                                  <  7:219615/139‑1 (MQ=255)
    ccGTTGCCGAAAATAAGTGATGTTTCGCTGACATTGAGGTGCTTACTGAGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAAc                                                                                                                                  <  4:321898/139‑1 (MQ=255)
     cGTTGCCGAAAATAAGTGATGTTTCGCTGACATTGAGGTGCTTACTGAGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACa                                                                                                                                 <  6:362923/139‑1 (MQ=255)
       ttGCCGAAAATAAGTGATGTTTCGCTGACATTGAGGTGCTTACTGAGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGc                                                                                                                               >  1:316573/1‑139 (MQ=255)
       ttGCCGAAAATAAGTGATGTTTCGCTGACATTGAGGTGCTTACTGAGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGc                                                                                                                               >  2:149508/1‑139 (MQ=255)
                        tGTTTCGCTGACATTGAGGTGCTTACTGAGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGACCAGCCGTAAGGATTCTCGttt                                                                                                              >  5:417149/1‑139 (MQ=255)
                          tttCGCTGACATTGAGGTGCTTACTGAGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTg                                                                                                             >  1:378698/1‑138 (MQ=255)
                                     ttGAGGTGCTTACTGAGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTca                                                                                                 <  3:339371/139‑1 (MQ=255)
                                             cTTACTGAGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCGTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTGCACAACTTTa                                                                                         >  3:252680/1‑139 (MQ=255)
                                              tcaCTGAGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGATAAAGCGCAAGCTGTTGAATCTCATGTTGAAGCGCCTCTTTTGCAGCTTCTGAACAGCCGTACGGATTTTCGTTTGAAGCAAGCTTCACAACTTTAt                                                                                        >  3:71254/3‑139 (MQ=255)
                                                aCTGAGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATc                                                                                       >  1:74132/1‑138 (MQ=255)
                                                 cTGAGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATcc                                                                                      <  2:393423/138‑1 (MQ=255)
                                                    aGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGaa                                                                                                                                   <  7:181474/90‑1 (MQ=255)
                                                    aGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGaa                                                                                                                                   >  8:181474/1‑90 (MQ=255)
                                                    aGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGc                                                                                  >  5:343869/1‑139 (MQ=255)
                                                       cTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCa                                                                                >  6:375403/1‑138 (MQ=255)
                                                         tGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCatat                                                                             >  7:5769/1‑139 (MQ=255)
                                                          gcccGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGccctgt                                                                             <  3:243991/136‑5 (MQ=255)
                                                          gTCCGCCCAGCGGCGCTATCCCCGACCCGGTAATGGGCAAGCTGTTGAATCTCATGATGAAGCGCCTCTATTGCAGCCTCTGAACAGCCGTAAGGATTGTCGTTTGAAGCTAGCTTCTCAACTTTATCCAAGCCATAtt                                                                            <  1:392868/139‑1 (MQ=255)
                                                          gTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATcc                                                                                      <  3:156131/129‑1 (MQ=255)
                                                          gTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATcc                                                                                      >  4:156131/1‑129 (MQ=255)
                                                           tCCGCAAAGCGGCGCTATACCCGGCCGGTTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATAtt                                                                            <  7:402591/138‑1 (MQ=255)
                                                           tCCGAAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTc                                                                           <  4:265370/139‑1 (MQ=255)
                                                               cAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGAt                                                                       <  6:343869/139‑1 (MQ=255)
                                                                aaaGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATcc                                                                                      <  4:92500/123‑1 (MQ=255)
                                                                aaaGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATcc                                                                                      >  3:92500/1‑123 (MQ=255)
                                                                   gcggcgCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCa                                                                   <  4:320531/139‑1 (MQ=255)
                                                                      gcgcTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTg                                                                <  2:271590/139‑1 (MQ=255)
                                                                      gcgcTATACCCGTCCGGGTAAAGGGCTAGCTGGTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTg                                                                <  6:287443/139‑1 (MQ=255)
                                                                        gcTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCt                                                              <  2:378698/139‑1 (MQ=255)
                                                                            taCCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGaa                                                                                                                                   <  6:246070/66‑1 (MQ=255)
                                                                            taCCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGaa                                                                                                                                   >  5:251565/1‑66 (MQ=255)
                                                                            taCCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGaa                                                                                                                                   <  6:251565/66‑1 (MQ=255)
                                                                            taCCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGaa                                                                                                                                   >  5:246070/1‑66 (MQ=255)
                                                                                 gtgcGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGc                                                     >  5:396423/4‑139 (MQ=255)
                                                                                 gTCCGGGTAAAGGGCTAGCTGTTGATTCTCATGATGAAGCGACTCTTTTGCAGCCTCTGCACAGCCGTAAGGATTCTCGTTTGAAGGAAGCTTCACAACTTTTGCCCAGCCATATTCTGATTTCACTGCTTCCAACggg                                                     >  2:217505/1‑138 (MQ=255)
                                                                                 gTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGaa                                                                                                                                   <  1:40688/61‑1 (MQ=255)
                                                                                 gTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGaa                                                                                                                                   >  2:40688/1‑61 (MQ=255)
                                                                                 gTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGc                                                     >  8:367868/1‑139 (MQ=255)
                                                                                 gTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGc                                                     >  1:20758/1‑139 (MQ=255)
                                                                                 gTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGc                                                     >  3:259809/1‑139 (MQ=255)
                                                                                 gTCCGGATAAAGCGCAAGCTGTTGAATCTCATGTTGAAGCGCCTCTTTTGCAGCTTCTGAACAGCCGTACGGATTTTCGTTTGAAGCAAGCTTca                                                                                                 <  5:215995/95‑1 (MQ=255)
                                                                                 gTCCGGATAAAGCGCAAGCTGTTGAATCTCATGTTGAAGCGCCTCTTTTGCAGCTTCTGAACAGCCGTACGGATTTTCGTTTGAAGCAAGCTTca                                                                                                 >  6:215995/1‑95 (MQ=255)
                                                                                  tCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCt                                                    >  3:247356/1‑139 (MQ=255)
                                                                                      ggTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGTATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTc                                                 >  7:59647/1‑138 (MQ=255)
                                                                                      ggTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTc                                                 >  5:307911/1‑138 (MQ=255)
                                                                                       gTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCa                                               <  2:316573/139‑1 (MQ=255)
                                                                                              gCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAGGCTGAt                                        <  5:174208/139‑1 (MQ=255)
                                                                                                       tGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGaagcaag                                                                                                      >  6:296093/1‑68 (MQ=255)
                                                                                                       tGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGaagcaag                                                                                                      <  5:296093/68‑1 (MQ=255)
                                                                                                       tGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTATGATTTCAATGCTTCTATCGGCTTTCCAGGCTGTTAAGTCTTcg                               >  3:16889/1‑138 (MQ=255)
                                                                                                        gAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAGGCTGATAAGGCTTCAg                              >  2:38160/1‑139 (MQ=255)
                                                                                                         atcctcATTTTGAAGGGGCTCTTTTGGAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCTAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAGGCTGATATGGCTTCAg                              <  3:401936/135‑1 (MQ=255)
                                                                                                         aaTCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAGGCTGATAAGGCTTCAg                              <  6:58172/138‑1 (MQ=255)
                                                                                                         aaTCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAGGCTGATAAGGCTTCAg                              <  8:5769/138‑1 (MQ=255)
                                                                                                         aaTCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAGGCTGATAAGGCTTCAg                              <  6:396423/138‑1 (MQ=255)
                                                                                                               atgatgAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAGGCTGATAAGGCTTCAGCTGt                          <  6:307911/136‑1 (MQ=255)
                                                                                                                  atgaAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGAAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAGGCTCATAAGGCTTCAGCTGTTTTaaa                    >  1:311927/1‑139 (MQ=255)
                                                                                                                     aaGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAGGCTGATAAGGCTTCAGCTGTTTTAAATGt                 >  8:120335/1‑139 (MQ=255)
                                                                                                                      aGCGCCTCTTTTGCAGCTTCTGAACAGCCGTACGGATTTTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAGGCTGATAAGGCTTCAGTTGTTTTAAATGtt                <  4:71254/139‑1 (MQ=255)
                                                                                                                      aGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAgg                                             >  2:323219/1‑110 (MQ=255)
                                                                                                                      aGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAgg                                             <  1:323219/110‑1 (MQ=255)
                                                                                                                           ctctTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAGGCTGATAAGGCTTCAGCTGTTTTAAATGTTCTTTg           <  4:247356/139‑1 (MQ=255)
                                                                                                                           ctctTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAGGCTGATAAGGCTTCAGCTGTTTTAAATGTTCTTTg           >  2:208565/1‑139 (MQ=255)
                                                                                                                                ttGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAGGCTGATAAGGCTTCAGCTGTTTTAAATGTTCTTTGAt         <  4:259809/136‑1 (MQ=255)
                                                                                                                                     gCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTc                                                                           <  5:129172/65‑1 (MQ=255)
                                                                                                                                     gCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTc                                                                           >  6:129172/1‑65 (MQ=255)
                                                                                                                                     gCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAGGCTGATAAGGCTTCAGCTGTTTTAAATGTTCTTTGATACGCaaa  >  3:21907/1‑138 (MQ=255)
                                                                                                                                       |                                                                                                                                       
AGAACCGTTGCCGAAAATAAGTGATGTTTCGCTGACATTGAGGTGCTTACTGAGCCTTGTCCGCAAAGCGGCGCTATACCCGTCCGGGTAAAGGGCTAGCTGTTGAATCTCATGATGAAGCGCCTCTTTTGCAGCCTCTGAACAGCCGTAAGGATTCTCGTTTGAAGCAAGCTTCACAACTTTATCCAAGCCATATTCTGATTTCACTGCTTCAATCGGCTTTCCAGGCTGATAAGGCTTCAGCTGTTTTAAATGTTCTTTGATACGCAAA  >  BS168_reference_genome/2371228‑2371498

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: