Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 3,290,481 C→T 5.5% M335I (ATG→ATA dhbC ← isochorismate synthase (siderophore‑specific)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome3,290,4810CT 5.5% 242.2 / 2.6 73M335I (ATG→ATAdhbCisochorismate synthase (siderophore‑specific)
Reads supporting (aligned to +/- strand):  ref base C (38/31);  new base T (2/2);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.39e-01

CTGAAGCTCGTCCTCAGGCTTTGAACCGGCCACAACACCAGCTCCAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCGGGTGGAGCGCCGCCGCCAATTCAAGTGCGGTTACAGACGGGTCGGA  >  BS168_reference_genome/3290346‑3290618
                                                                                                                                       |                                                                                                                                         
cTGAAGCTCGGCCTCATGCTTTGAACCGGCCACAACACCAGCGCCAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGTATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATg                                                                                                                                        <  5:5043/139‑1 (MQ=255)
          tCCTCAGGCTTTGAACCGGCCACAACACCAGCTCCAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAg                                                                                                                               >  7:142335/1‑138 (MQ=255)
          tCCTCAGGCTTTGAACCGGCCACAACACCAGCTCCAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAg                                                                                                                               >  5:190480/1‑138 (MQ=255)
          tCCTCAGGCTTTGAACCGGCCACAACACCAGCTCCAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAg                                                                                                                               >  6:125339/1‑138 (MQ=255)
          tCCTCAGGCTTTGAACCGGCCACAACACCAGCTCCAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAg                                                                                                                               >  1:163400/1‑138 (MQ=255)
          tCCTCAGGCTTTGAACCGGCCACAACACCAGCTCCAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCAACGGCATCGTCACACCATCCGACCATGCCGGTAAAg                                                                                                                               >  5:52295/1‑138 (MQ=255)
          tCCTCAGGCTTTGAACCGGCCACAACACCAGCTCCAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACCGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCCTCCGCCCATGCCGGTAAAg                                                                                                                               >  3:184454/1‑138 (MQ=255)
               aGGCTTTGAACCGGCCACAACACCAGCTCCAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACcg                                                                                                                         <  6:190480/139‑1 (MQ=255)
               aGGCTTTGAACCGGCCACAACACCAGCTCCAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACcg                                                                                                                         <  2:163400/139‑1 (MQ=255)
               aGGCTTTGAACCGGCCACAACACCAGCTCCAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACcg                                                                                                                         <  8:142335/139‑1 (MQ=255)
                       aaCCGGCCACAACACCAGCTCCAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCaaa                                                                                                                 <  5:125339/139‑1 (MQ=255)
                            gCCACAACACCAGCTCCAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGtt                                                                                                            >  8:226121/1‑139 (MQ=255)
                                      cAGCTCCAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCt                                                                                                     >  5:51152/1‑136 (MQ=255)
                                      cAGCTCCAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCt                                                                                                     >  2:238818/1‑136 (MQ=255)
                                           ccAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCaaa                                                                                                                 <  7:255585/119‑1 (MQ=255)
                                           ccAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCaaa                                                                                                                 >  8:255585/1‑119 (MQ=255)
                                                      gCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTGGCTTGc                                                                                   >  1:151647/1‑138 (MQ=255)
                                                      gCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGc                                                                                   >  1:247533/1‑138 (MQ=255)
                                                      gCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCa                                                                                  <  3:68369/139‑1 (MQ=255)
                                                         gAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGtt                                                                                                            <  1:279768/110‑1 (MQ=255)
                                                         gAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGtt                                                                                                            >  2:279768/1‑110 (MQ=255)
                                                         gAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGa                                                                               <  7:226121/139‑1 (MQ=255)
                                                           gTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATc                                                                             <  5:29208/139‑1 (MQ=255)
                                                            tGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGAt                                                                              >  3:183122/1‑137 (MQ=255)
                                                            tGAGCGTTCTTCTGCGTCTGCACAACGGATGGTCACGCTCCATTCTCCGTCACCGGCATCGTCACACCAGCCGACCATGCCGGTAAAGAACCCGCGGTCAAATGGTTCACTGCTGAGAATAGCTTCGAGTGCAAGCTcc                                                                            >  2:249946/1‑138 (MQ=255)
                                                              aGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGt                                                                          <  8:167737/139‑1 (MQ=255)
                                                              aGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGt                                                                          <  1:249946/139‑1 (MQ=255)
                                                                cGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCg                                                                        >  2:35760/1‑139 (MQ=255)
                                                                cGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCg                                                                        <  6:51152/139‑1 (MQ=255)
                                                                    cttctgcttctgCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGgtg                                                                     <  6:261625/138‑1 (MQ=255)
                                                                    cttctgcttctgCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGGAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGgtgt                                                                    >  8:162586/1‑139 (MQ=255)
                                                                       ctgcttctgcACAAGGGATGGGCACGATCCATGCTCCGTCACCGGCTGCGTCACACCATCCGACCATGCCGGTAAAGAACCCGCGGTAAAATGGTTCATTGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTc                                                                  <  2:125697/138‑1 (MQ=255)
                                                                          cttctgcACAACGGATGGTCACGATCCATTCTCCGTCACCGGCAACGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCa                                                              >  1:125697/1‑139 (MQ=255)
                                                                             ctgcACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGAc                                                           >  5:249058/1‑139 (MQ=255)
                                                                             ctgcACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGAc                                                           >  7:208844/1‑139 (MQ=255)
                                                                               gcACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCg                                                                        <  2:112997/124‑1 (MQ=255)
                                                                               gcACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCg                                                                        <  2:250125/124‑1 (MQ=255)
                                                                               gcACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCg                                                                        >  1:250125/1‑124 (MQ=255)
                                                                               gcACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCg                                                                        >  1:112997/1‑124 (MQ=255)
                                                                                     cGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCg                                                   >  5:125983/1‑139 (MQ=255)
                                                                                          ggTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCggg                                              >  8:264649/1‑139 (MQ=255)
                                                                                                gATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGTTGTTCCGCAGCCGGCTGGCGGCGGGTGGAgc                                        >  8:145841/1‑139 (MQ=255)
                                                                                                 aTCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCGGGTGGAgcg                                       <  1:238818/139‑1 (MQ=255)
                                                                                                  tCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCGGGTGGAgcgc                                      >  8:207974/1‑139 (MQ=255)
                                                                                                  tCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCGGGTGGAgcgc                                      >  6:131462/1‑139 (MQ=255)
                                                                                                   ccATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCGGGTGGAGcgcc                                     <  2:287344/139‑1 (MQ=255)
                                                                                                      ttCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGtt                                                                                                            <  8:117796/65‑1 (MQ=255)
                                                                                                      ttCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGtt                                                                                                            >  7:117796/1‑65 (MQ=255)
                                                                                                           cGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTACGCAGACGGCTGGCGACCGGGGGAGCCCCCCCCCCa                              >  2:99677/1‑138 (MQ=255)
                                                                                                              cACCGGCATCGTCACACCATCAGACCAAGCCGGTAAAGAAAACGCGGTCAAATGGTTCAATGCTGAGACTCGCTTCTCTTGCAACATCAGTCGGTGTTACGCAGACGGCAGGCGTCGGCTGGAGCGCCCCCGCCaaaa                           >  6:139987/1‑136 (MQ=255)
                                                                                                                ccGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCGGGTGGAGCGCCGCCGCCAATTc                          >  6:79171/1‑137 (MQ=255)
                                                                                                                   gCATCGTCACCCCATCCGACCATGCCGTTAAAGAAACCGCGGTCAAATGGTGCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCGGGTGGAGCGCCGCCGCCAATTCAAGTg                     <  2:202107/139‑1 (MQ=255)
                                                                                                                   gCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCGGGTGGAGCGCCGCCGCCAATTCAAGTg                     <  8:208844/139‑1 (MQ=255)
                                                                                                                   gCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCGGGTGGAGCGCCGCCGCCAATTCAAGTg                     <  6:52295/139‑1 (MQ=255)
                                                                                                                   gCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCGGGTGGAGCGCCGCCGCCAATACAAGTg                     <  3:130751/139‑1 (MQ=255)
                                                                                                                             acCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTg                                                                                    <  2:228829/66‑1 (MQ=255)
                                                                                                                             acCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTg                                                                                    >  1:228829/1‑66 (MQ=255)
                                                                                                                             acCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCGGGTGGAGCGCCGCCGCCAATTCAAGTGCGGTTACAGa           <  5:79171/139‑1 (MQ=255)
                                                                                                                             acCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCGGGTGGAGCGCCGCCGCCAATTCAAGTGCGGTTACAGa           <  1:62581/139‑1 (MQ=255)
                                                                                                                                 tCCGACTATGCCGGTAAAGAAACCGCGGTCAAATGGtt                                                                                                            >  6:74026/1‑38 (MQ=255)
                                                                                                                                 tCCGACTATGCCGGTAAAGAAACCGCGGTCAAATGGtt                                                                                                            >  6:78324/1‑38 (MQ=255)
                                                                                                                                 tCCGACTATGCCGGTAAAGAAACCGCGGTCAAATGGtt                                                                                                            <  5:78324/38‑1 (MQ=255)
                                                                                                                                 tCCGACTATGCCGGTAAAGAAACCGCGGTCAAATGGtt                                                                                                            <  5:74026/38‑1 (MQ=255)
                                                                                                                                 tCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGc                                                        <  4:213171/90‑1 (MQ=255)
                                                                                                                                 tCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGc                                                        >  3:213171/1‑90 (MQ=255)
                                                                                                                                 tCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCGGGTGGAGCGCCGCCGGCAAATCAAGTGCGGTTCCAGACggg       >  6:243286/1‑139 (MQ=255)
                                                                                                                                 tCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGGCGGGTGGCGCGCCGCCGCCCATTAAAATGCGGGTACACACcgg       >  2:62581/1‑139 (MQ=255)
                                                                                                                                 tCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACCGCTTGCGGCGGGTGGGGCGCCCCCGCCAATTCAAATGCGGTTTCACACggg       >  6:259215/1‑139 (MQ=255)
                                                                                                                                  ccGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGt                                                                          >  5:216836/1‑71 (MQ=255)
                                                                                                                                  ccGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGt                                                                          <  6:216836/71‑1 (MQ=255)
                                                                                                                                  ccGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCGGGTGGAGCGCCGCCGCCAATTCAAGTGCGGTTACAGACGGGt      >  2:100853/1‑139 (MQ=255)
                                                                                                                                    gACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCGGGTGGAGCGCCGCCGCCAATTCAAGTGCGGTTACAGACGGGTCg    <  1:99677/139‑1 (MQ=255)
                                                                                                                                    gACCATCCCGCAAAAGACCCCGCGGGCAACTGGTGCCATGCGGAGAATCGCTTCGCTCGCCCGAGCCGGCGGGGGGCCGCAGACGGCGGGCGTCGGGTGGAGCGCCGCCGCCAATTCCAGGGCGGTTACCGACGGGTCg    <  4:174157/139‑1 (MQ=255)
                                                                                                                                      ccATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGGCGGGTGGAGCGCCGCCGCCAATTCAAGTGCGGTTACAGACGGGTCGGa  <  8:150831/139‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                         
CTGAAGCTCGTCCTCAGGCTTTGAACCGGCCACAACACCAGCTCCAGCATACAGGCGGAGTGAGCGTTCTTCTGCTTCTGCACAACGGATGGTCACGATCCATTCTCCGTCACCGGCATCGTCACACCATCCGACCATGCCGGTAAAGAAACCGCGGTCAAATGGTTCAATGCTGAGAATCGCTTCTCTTGCAAGATCAGTCGGTGTTCCGCAGACGGCTGGCGTCGGGTGGAGCGCCGCCGCCAATTCAAGTGCGGTTACAGACGGGTCGGA  >  BS168_reference_genome/3290346‑3290618

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: