Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,610,685 C→T 6.9% V88I (GTC→ATC)  era ← maturation of 16S RNA and assembly of 30S ribosomal subunit GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,610,6850CT 6.9% 186.0 / 3.0 58V88I (GTC→ATC) eramaturation of 16S RNA and assembly of 30S ribosomal subunit GTPase
Reads supporting (aligned to +/- strand):  ref base C (25/29);  new base T (2/2);  total (27/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

ATTGGTCGGGATGAATTTTATCAATTTTGTTTACGATTAAGAATACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTACTCCTTGCACCTTGTTTCTCGT  >  BS168_reference_genome/2610547‑2610820
                                                                                                                                          |                                                                                                                                       
aTTGGTCGGTATGAATTTTATCAATTTTGTTTACGATTAAGACTACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGAc                                                                                                                                         <  7:271207/139‑1 (MQ=255)
aTTGGTCGGGATGAATTTTATCAATTTTGTTTACGATTAAGAATACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGAc                                                                                                                                         <  6:96021/139‑1 (MQ=255)
aTTGGTCGGGATGAATTTTATCAATTTTGTTTACGATTAAGAATACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGAc                                                                                                                                         <  2:296909/139‑1 (MQ=255)
 ttGGTCGGGATGAATTTTATCAATTTTGTTTACGATTAAGAATACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGAct                                                                                                                                        >  8:62187/1‑139 (MQ=255)
       gggATGAATTTTATCAATTTTGTTTACGATTAAGAATACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGACTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCttt                                                                                                                                   >  7:299191/1‑138 (MQ=255)
         gATGAATTTTATCAATTTTGTTTACGATTAAGAATACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGAc                                                                                                                                         >  7:260547/1‑130 (MQ=255)
         gATGAATTTTATCAATTTTGTTTACGATTAAGAATACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGAc                                                                                                                                         >  1:306850/1‑130 (MQ=255)
         gATGAATTTTATCAATTTTGTTTACGATTAAGAATACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGAc                                                                                                                                         <  2:306850/130‑1 (MQ=255)
         gATGAATTTTATCAATTTTGTTTACGATTAAGAATACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGAc                                                                                                                                         <  8:260547/130‑1 (MQ=255)
                    tCAATTTTGTTTACGATTAAGAATACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTgcg                                                                                                                     <  8:23329/139‑1 (MQ=255)
                       aTTTTGTTTACGATTAAGAATACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCa                                                                                                                  >  3:52093/1‑139 (MQ=255)
                                  gATTAAGAATACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTcatcat                                                                                                        >  5:254617/1‑138 (MQ=255)
                                  gATTAAGAATACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATa                                                                                                       <  4:168548/139‑1 (MQ=255)
                                         aaTACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTccc                                                                                                <  2:156643/139‑1 (MQ=255)
                                          aTACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTcccc                                                                                               <  1:87848/139‑1 (MQ=255)
                                          aTACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTccc                                                                                                <  8:145433/138‑1 (MQ=255)
                                              aGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACt                                                                                                                <  2:301250/118‑1 (MQ=255)
                                              aGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACt                                                                                                                >  1:301250/1‑118 (MQ=255)
                                              aGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGt                                                                                           <  1:88843/139‑1 (MQ=255)
                                              aGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTAAGCATAGATCATAAACAAGATCAAATCGACTACTTTTAATGTATTTTGCTCAAGTTTCATCATAAAGTCCCCaact                                                                                           >  7:10698/1‑137 (MQ=255)
                                                   tAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGAc                                                                                                                                         <  8:143328/88‑1 (MQ=255)
                                                   tAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGAc                                                                                                                                         >  7:143328/1‑88 (MQ=255)
                                                   tAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATg                                                                                      >  2:62163/1‑139 (MQ=255)
                                                   tAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATg                                                                                      >  8:311602/1‑139 (MQ=255)
                                                   tAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGAGCCACTTTCATCATAAATTCCCCAACTTTATg                                                                                      >  4:155157/1‑139 (MQ=255)
                                                     gACATCGTCTGGAATTTCTCAATGTTGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAATATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCt                                                                                    <  7:47264/139‑1 (MQ=255)
                                                     gACATCGTCTGCAATTTCTCAATGATGACTTCGTCGCCCGTCCCGTAGCCGTCCTCAGCAATGAGCATCAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCt                                                                                    <  6:168686/139‑1 (MQ=255)
                                                     gACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGATTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCt                                                                                    >  6:81722/1‑139 (MQ=255)
                                                     gACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGATTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCt                                                                                    >  6:75863/1‑139 (MQ=255)
                                                          cGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAATATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATa                                                                                                       >  3:244668/1‑115 (MQ=255)
                                                          cGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGt                                                                               >  7:275812/1‑139 (MQ=255)
                                                          cGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGt                                                                               >  5:36468/1‑139 (MQ=255)
                                                          cGTCTGCAATTTCTCAATGATGAATTCGTAGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATa                                                                                                       <  4:244668/115‑1 (MQ=255)
                                                                cAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATg                                                                          <  6:36468/138‑1 (MQ=255)
                                                                          atgatgAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGgtg                                                               >  2:358063/1‑139 (MQ=255)
                                                                          atgatgAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGgtg                                                               <  1:89837/139‑1 (MQ=255)
                                                                              tgaATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGgtgtgt                                                            >  7:141049/1‑138 (MQ=255)
                                                                              tgaATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTc                                                           <  1:235985/139‑1 (MQ=255)
                                                                                     gtcgCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATaaaa                                                    <  7:62187/139‑1 (MQ=255)
                                                                                     gtcgCCCTTCCCGTACCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACGTCTTTTAATGAATTTTGCGACACTTTCTTCTTAAAGTCACCAAGATTATGCTTAGGTTTATCACTCCCTGGTGTGTCtatata                                                     >  8:96657/1‑136 (MQ=255)
                                                                                     ggcgCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATaaaa                                                    <  4:54801/137‑1 (MQ=255)
                                                                                          ccTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGt                                               <  8:141049/139‑1 (MQ=255)
                                                                                                cGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGa                                         <  1:358063/139‑1 (MQ=255)
                                                                                                         ccTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTgt                                >  3:346101/1‑139 (MQ=255)
                                                                                                         ccTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCACTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCttg                                >  7:264193/1‑137 (MQ=255)
                                                                                                          cTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTgt                                <  8:275812/138‑1 (MQ=255)
                                                                                                                aTTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGt                         <  4:346101/139‑1 (MQ=255)
                                                                                                                 ttGATCATAAACAAGATCAAATCGATTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTa                        <  5:81722/139‑1 (MQ=255)
                                                                                                                 ttGATCATAAACAAGATCAAATCGATTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTa                        <  5:75863/139‑1 (MQ=255)
                                                                                                                 ttGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTa                        <  7:311602/139‑1 (MQ=255)
                                                                                                                  tGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATATAGTCCCGAAGTTTTTGGTTAGGTTTATTATTCCCTGGGGTGTCATTTAAAATGGTTTGTGATGTGCCTGTCGTCTGTAc                       >  2:307394/1‑139 (MQ=255)
                                                                                                                  tGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTAc                       >  5:79954/1‑139 (MQ=255)
                                                                                                                    aTCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTccc                                                                                                >  4:44196/1‑64 (MQ=255)
                                                                                                                    aTCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTccc                                                                                                <  3:44196/64‑1 (MQ=255)
                                                                                                                    aTCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTACTc                     >  4:218241/1‑139 (MQ=255)
                                                                                                                     tCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTACTcc                    <  1:62163/139‑1 (MQ=255)
                                                                                                                        tAAACAATATCAAATCGACTTCTTTTAATGTATTATGCGCCACATTCATCATAAAGGCCCCAAGTTTATGCTTATGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTACTCCTTg                 <  4:28112/139‑1 (MQ=255)
                                                                                                                         aaaCAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTa                        >  2:4978/1‑131 (MQ=255)
                                                                                                                         aaaCAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTa                        <  1:4978/131‑1 (MQ=255)
                                                                                                                         aaaCAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTACTCCTTg                 >  5:252963/1‑138 (MQ=255)
                                                                                                                         aaaCAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTACTCCTTGc                >  8:276544/1‑139 (MQ=255)
                                                                                                                               gATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTAc                       <  2:116609/126‑1 (MQ=255)
                                                                                                                               gATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTAc                       >  1:116609/1‑126 (MQ=255)
                                                                                                                               gATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTACTCCTTGCACCTTg          >  7:123481/1‑139 (MQ=255)
                                                                                                                                      tCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTACTCCTTGCACCTTGTTTCtcg   <  3:155157/139‑1 (MQ=255)
                                                                                                                                      tCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTACTCCTTGCACCTTGTTTCtcg   <  8:123481/139‑1 (MQ=255)
                                                                                                                                      tCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTACTCCTTGCACCTTGTTTCtcg   <  6:254617/139‑1 (MQ=255)
                                                                                                                                        gACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTACTCCTTGCACCTTGTTTCtcgt  <  7:276544/138‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                       
ATTGGTCGGGATGAATTTTATCAATTTTGTTTACGATTAAGAATACAGGCGTAGACATCGTCTGCAATTTCTCAATGATGAATTCGTCGCCCTTCCCGTAGCCTTCCTCAGCATTGATCATAAACAAGATCAAATCGACTTCTTTTAATGTATTTTGCGCCACTTTCATCATAAAGTCCCCAAGTTTATGCTTAGGTTTATGAATCCCTGGTGTGTCAATAAAAATGGTTTGTGATGTGCCTGTCGTCAGTACTCCTTGCACCTTGTTTCTCGT  >  BS168_reference_genome/2610547‑2610820

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: