Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,411,036 A→G 40.7% intergenic (‑447/‑50) fmnP ← / → serA FMN permease/3‑phosphoglycerate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,411,0360AG40.7% 40.4 / 102.7 86intergenic (‑447/‑50)fmnP/serAFMN permease/3‑phosphoglycerate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (28/22);  new base G (23/11);  total (51/35)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.98e-02

ACGATAACTGGTTGGCTGTATAGAGAATTGCTTCAAATTAAATATATTACCTGCAAGCTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACTTATTGAATCAGACTTTATTGAAATCGTTCAAAAAAA  >  BS168_reference_genome/2410903‑2411174
                                                                                                                                     |                                                                                                                                          
aCGATAACTGGTTGGCTGTATAGAGAATTGCTTCAAATTAAATATATTACCTGCAAGCTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTa                                                                                                                                       <  6:221728/139‑1 (MQ=255)
      tgcgttcgtttgaaatcatcgctgacgtttggaaaaacagataacagctataaaaaaagcaatgtagcggggagacccatacattgcttTTTACTGCGTCAATACACGTTGACACTCTTTTGAAAGTGTGTTAAATTAt                                                                                                                                 <  5:106055‑M1/50‑1 (MQ=255)
          ttcgtttgaaatcatcgctgacgtttggaaaaacagataacagctataaaaaaagcaatgtagcggggagacccatacattgcttTTTACTGCGTCAATACACGTTGACACTCTTTTGAAAGTGTGTTAAATTATCagt                                                                                                                             <  6:53325‑M1/54‑2 (MQ=255)
          ttcgtttgaaatcatcgctgacgtttggaaaaacagataacagctataaaaaaagcaatgtagcggggagacccatacattgcttTTTACTGCGTCAATACACGTTGACACTCTTTTGAAAGAGTGTTAAATTATCagt                                                                                                                             <  1:55789‑M1/54‑2 (MQ=255)
            tGGCTGTATAGAGAATTGCTTCAAATTAAATATATTACCTGCAAGCTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGata                                                                                                                           >  8:92859/1‑139 (MQ=255)
                  tataGAGAATTGCTTCAAATTAAATATATTACCTGCAAGCTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGt                                                                                                                     >  1:270672/1‑139 (MQ=255)
                  tataGAGAATTGCTTCAAATTAAATATATTACCTGCAAGCTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGt                                                                                                                     >  6:366212/1‑139 (MQ=255)
                       agaATTGCTTCAAATTAAATATATTACCTGCAAGCTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCATTACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGt                                                                                                                 >  3:200605/1‑138 (MQ=255)
                       agaATTGCTTCAAATTAAATATATTACCTGCAAGCTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCATTACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGt                                                                                                                 >  3:193757/1‑138 (MQ=255)
                         aaTTGCTTCAAATTAAATATATTACCTGCAAGCTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCaa                                                                                                              >  5:165789/1‑139 (MQ=255)
                         aaTTGCTTCAAATTAAATATATTACCTGCAAGCTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCaa                                                                                                              >  5:186067/1‑139 (MQ=255)
                            tGCTTCAAATTAAATATATTACCTGCAAGCTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAAtt                                                                                                            >  3:277927/1‑138 (MQ=255)
                              cTTCAAATTAAATATATTACCTGCAAGCTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTag                                                                                                         >  3:6054/1‑139 (MQ=255)
                                   aaTTAAATATATTACCTGCAAGCTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTaggagg                                                                                                     >  1:226259/1‑138 (MQ=255)
                                              ttACCTGCAAGCTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGAt                                                                                         <  4:407366/139‑1 (MQ=255)
                                                    gCAAGCTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTc                                                                                    >  2:113983/1‑138 (MQ=255)
                                                         cTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTAt                                                                              <  4:200605/139‑1 (MQ=255)
                                                         cTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTAt                                                                              <  4:193757/139‑1 (MQ=255)
                                                         cTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTAt                                                                              <  2:270672/139‑1 (MQ=255)
                                                            tCAGATCAGTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTGTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAAGTTAGGAGGAAATCTATCGATGTTTCGAGTATTgg                                                                           >  3:245244/1‑139 (MQ=255)
                                                               gATCATTGATTTATTAGGCTTTACTTTTATACTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGtct                                                                        <  4:250278/139‑1 (MQ=255)
                                                                  cATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGtct                                                                        >  1:200810/1‑136 (MQ=255)
                                                                    ttGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGTAAATCTAACGATGTTTCGAGTATTGGTCTCAGa                                                                    >  1:151149/1‑138 (MQ=255)
                                                                            ttaGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAAc                                                                                            <  1:345484/106‑1 (MQ=255)
                                                                            ttaGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAAc                                                                                            >  2:345484/1‑106 (MQ=255)
                                                                              acgcgttACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGc                                                          <  2:367918/133‑1 (MQ=255)
                                                                              aGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGc                                                          <  6:186067/138‑1 (MQ=255)
                                                                                 cTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAac                                                       >  2:113804/1‑138 (MQ=255)
                                                                                  tttACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAacga                                                     >  3:54950/1‑139 (MQ=255)
                                                                                   ttACTTTTATCATTTACTGCGTCAATACACTTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGTAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAacga                                                     <  4:342543/138‑1 (MQ=255)
                                                                                    tACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAATTCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAacgacg                                                   >  6:378097/1‑139 (MQ=255)
                                                                                    tACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAacgacg                                                   >  7:394830/1‑139 (MQ=255)
                                                                                    tACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAacgacg                                                   >  8:276804/1‑139 (MQ=255)
                                                                                          tgctttttactgcGTCAATACACGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattc                                              >  7:249172‑M1/6‑58 (MQ=255)
                                                                                          tgctttttactgcGTCAATACACGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattc                                              >  5:383876‑M1/6‑58 (MQ=255)
                                                                                          tgctttttactgcGTCAATACACGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattc                                              >  7:64979‑M1/6‑58 (MQ=255)
                                                                                          tgctttttactgcGTCAATACACGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattc                                              >  7:344251‑M1/6‑58 (MQ=255)
                                                                                          tgctttttactgcGTCAATACACGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgt                                                             <  6:326251‑M1/118‑66 (MQ=255)
                                                                                          tgctttttactgcGTCAATACACGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgt                                                             >  5:326251‑M1/6‑58 (MQ=255)
                                                                                            tCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACa                                            <  4:408506/138‑1 (MQ=255)
                                                                                            tCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACa                                            <  7:92859/138‑1 (MQ=255)
                                                                                            ctttttactGCGTCAATACACGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaattgaagaatactccaactcaattagaacctaatgttcctgtaacaagaagccattcaat                                           >  6:47883‑M1/4‑56 (MQ=255)
                                                                                            ctttttactGCGTCAATACACGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaat                                           >  6:367979‑M1/4‑56 (MQ=255)
                                                                                              cTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTAAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTACCGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGc                                          <  6:88680/138‑1 (MQ=255)
                                                                                               tttACTGCGTCAATACACGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgg                                         >  5:48752‑M1/1‑53 (MQ=255)
                                                                                               tttACTGCGTCAATACACGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgg                                         >  6:166788‑M1/1‑53 (MQ=255)
                                                                                                ttACTGCGTCAATACACTTTGACAATCTTTTGTTAATTTTTAAAATTATAATATTTTTAGTTTTTCAATTTTGGAGGAAATCTAACGATGTTTGTAGTATTGGTCTCAGAAAAAATGAGCCACGACGGTTTACAGCCa                                        <  2:200810/138‑1 (MQ=255)
                                                                                                 TACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACTGCCACt                                      <  6:90900‑M1/139‑1 (MQ=255)
                                                                                                 TACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACt                                      <  2:13595‑M1/139‑1 (MQ=255)
                                                                                                 TACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACt                                      <  7:6015‑M1/139‑1 (MQ=255)
                                                                                                 TACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACt                                      <  6:204708‑M1/139‑1 (MQ=255)
                                                                                                   cTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAacga                                                     >  8:385414/1‑122 (MQ=255)
                                                                                                   cTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAacga                                                     <  7:385414/122‑1 (MQ=255)
                                                                                                   cTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACTTa                                    <  1:113804/139‑1 (MQ=255)
                                                                                                      cGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCa                                        <  3:300470/132‑1 (MQ=255)
                                                                                                      cGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCa                                        >  4:300470/1‑132 (MQ=255)
                                                                                                      cGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACTTAtt                                  >  8:321393/1‑138 (MQ=255)
                                                                                                      cGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGAAACGAAGGTTTACAGCCACTTAtt                                  >  6:205160/1‑138 (MQ=255)
                                                                                                      cGTCAATACACGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgggatttgt                                  >  1:183387‑M1/1‑46 (MQ=255)
                                                                                                      cGTCAATACACGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgggatttg                                   <  5:166788‑M1/137‑92 (MQ=255)
                                                                                                      cGTCAATACACGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatggcatttgt                                  >  3:313234‑M1/1‑46 (MQ=255)
                                                                                                          aaTACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACTTATTGAATc                             <  7:321393/139‑1 (MQ=255)
                                                                                                            tACACGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgggatttgtcattttg                           <  8:64979‑M1/139‑100 (MQ=255)
                                                                                                            tACACGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgggatttgtcattttg                           <  6:117336‑M1/139‑100 (MQ=255)
                                                                                                              ctcgtTGACACTGTTTTTATATTTTTTTTAATTATCAGtgaaatgaatgtaattattttaaaattaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgggatttgtcattttgat                         <  2:355952‑M1/137‑102 (MQ=255)
                                                                                                              cacGTTGACACTATTTTGAAAGTGTTTTAAATTATCAGtgacattaaaggaagtattttaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgggatttgtcattttgat                         <  6:48752‑M1/139‑102 (MQ=255)
                                                                                                                cGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgggatttgtcattttgatct                       >  5:140936‑M1/1‑36 (MQ=255)
                                                                                                                cGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgggatttgtcattttgatct                       >  5:170227‑M1/1‑36 (MQ=255)
                                                                                                                cGTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgggatttgtcattttgatc                        >  6:267350‑M1/1‑36 (MQ=255)
                                                                                                                 gTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgggatttgtcattttgatctc                      <  4:314622‑M1/139‑105 (MQ=255)
                                                                                                                 gTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgggatttgtcattttgatct                       >  6:18103‑M1/1‑35 (MQ=255)
                                                                                                                 gTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgggatttgtcattttgatct                       >  8:308788‑M1/1‑35 (MQ=255)
                                                                                                                 gTTGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgatcctgtaacaagaagccattcaatgggatttgtcattttgatct                       >  7:389135‑M1/1‑35 (MQ=255)
                                                                                                                   tGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctg                                                              >  8:315105‑M1/1‑33 (MQ=255)
                                                                                                                   tGACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctg                                                              <  7:315105‑M1/97‑65 (MQ=255)
                                                                                                                    gACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgggatttgtcattttgatctcatg                   >  4:358319‑M1/1‑32 (MQ=255)
                                                                                                                    gACACTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagacattcaatgggatttgtcattttgatctcatg                   >  2:70499‑M1/1‑32 (MQ=255)
                                                                                                                     acacTCTTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgggatttgtcattttgatctcatg                   <  5:267350‑M1/138‑108 (MQ=255)
                                                                                                                        ctctTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgggatttgtcattttgatctcatgtgcg               >  1:107901‑M1/1‑28 (MQ=255)
                                                                                                                         tctTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccacctcaattagaaccaaatgttcctgtaacaagaagccattcaaggggatttgtcattttgatctcatgtgcgg              >  2:45058‑M1/1‑27 (MQ=255)
                                                                                                                         tctTTTGAAAGTGTGTTAAATTATCAGtgacatgaaaggaagtatttgaaaatgaagaatactccaactcaattagaaccaaatgttcctgtaacaagaagccattcaatgggatttgtcattttgatctcatgtgcgg              >  2:94041‑M1/1‑27 (MQ=255)
                                                                                                                           ttttGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCa                                        <  7:43809/111‑1 (MQ=255)
                                                                                                                           ttttGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCa                                        >  8:43809/1‑111 (MQ=255)
                                                                                                                           ttttGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACTTATTGAATCAGACTTTATTGAATTCg            >  2:219734/1‑139 (MQ=255)
                                                                                                                           ttttGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACTTATTGAATCAGACTTTATTGAAATCg            >  7:294515/1‑139 (MQ=255)
                                                                                                                           ttttGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACTTATTGAATCAGACTTTATTGAAATCg            >  4:188728/1‑139 (MQ=255)
                                                                                                                                   atatGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACTTATTGAATCAGACTTTATTGAAATCGTTCaaaa     >  7:408032/1‑138 (MQ=255)
                                                                                                                                   atatGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACTTATTGAATCAGACTTTATTGAAATCGTTCaaaa     >  3:168676/1‑138 (MQ=255)
                                                                                                                                     atGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACTTATTGAATCAGAc                         <  8:316960/116‑1 (MQ=255)
                                                                                                                                     atGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACTTATTGAATCAGACTTTATTGAAATCGTTCaaaaaaa  >  7:235513/1‑139 (MQ=255)
                                                                                                                                     atGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACTTATTGAATCAGACTTTATTGAAATCGTTCaaaaaaa  >  7:241557/1‑139 (MQ=255)
                                                                                                                                     atGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACTTATTGAATCAGACTTTATTGAAATCGTTCaaaa     >  7:15006/1‑136 (MQ=255)
                                                                                                                                     atGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACTTATTGAATCAGACTTTATTGAAATCGTTCaaaa     <  8:15006/136‑1 (MQ=255)
                                                                                                                                     atGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTAGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACTTATTGAATCAGAc                         >  7:316960/1‑116 (MQ=255)
                                                                                                                                     |                                                                                                                                          
ACGATAACTGGTTGGCTGTATAGAGAATTGCTTCAAATTAAATATATTACCTGCAAGCTGTCAGATCATTGATTTATTAGGCTTTACTTTTATCCTTTACTGCGTCAATACACGTTGACACTCTTTTGAGAATATGTTAAATTATCAGATATTTAGTTTGTCAATTTAGGAGGAAATCTAACGATGTTTCGAGTATTGGTCTCAGACAAAATGAGCAACGACGGTTTACAGCCACTTATTGAATCAGACTTTATTGAAATCGTTCAAAAAAA  >  BS168_reference_genome/2410903‑2411174

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: