Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,898,609 G→A *196* (TAG→TAA cbrC → UPF0167 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,898,6090GA98.7% 236.3 / NA 76*196* (TAG→TAAcbrCUPF0167 family protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (38/37);  total (39/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 = 1.00e+00

AAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACG  >  NC_000913/3898476‑3898743
                                                                                                                                     |                                                                                                                                      
aagaagATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAattatt                                                                                                                                   <  1:308274/139‑1 (MQ=255)
      aTTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaata                                                                                                                             >  1:490140/1‑139 (MQ=255)
        tgtgAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaataat                                                                                                                           >  1:850598/1‑139 (MQ=255)
          tgAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTAGTTAAATAATga                                                                                                                         >  2:1189118/1‑139 (MQ=255)
             gaATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGaa                                                                                                                      <  2:850598/139‑1 (MQ=255)
                    taagagacagaTTCTGATCTAGCTAAATGCCTGCAAAAGGCTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGGCg                                                                                                               <  2:1068848/129‑1 (MQ=255)
                         tAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCg                                                                                                          <  1:140411/139‑1 (MQ=255)
                           gAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTa                                                                                                        <  1:1189118/139‑1 (MQ=255)
                           gAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTa                                                                                                        >  2:1183829/1‑139 (MQ=255)
                           gAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTa                                                                                                        <  1:248871/139‑1 (MQ=255)
                                tCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTca                                                                                                   <  2:875836/139‑1 (MQ=255)
                                  tGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                                <  1:713763/108‑1 (MQ=255)
                                  tGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                                >  2:713763/1‑108 (MQ=255)
                                        aGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCg                                                                                           <  2:124884/139‑1 (MQ=255)
                                         gCTAAATTCCTGCAAAAGGGTGGTCATTGTCAGTGTTTTCTCTTCCGCTGTCTCCACTGCGGCAAGCAGAGATCTTGGGGTGATCTGTCGTAGTCATTGAAATAATGAGAACAGACCAGAGCGTAATTcctgtctctta                                                                                          >  1:1068848/1‑129 (MQ=255)
                                              aTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaataat                                                                                                                           <  1:70504/101‑1 (MQ=255)
                                              aTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAaataat                                                                                                                           >  2:70504/1‑101 (MQ=255)
                                              aTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTa                                                                                     <  2:1319547/139‑1 (MQ=255)
                                              aTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTa                                                                                     >  2:1665505/1‑139 (MQ=255)
                                              aTGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTa                                                                                     >  1:1319547/1‑139 (MQ=255)
                                                   tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                                <  1:1437020/91‑1 (MQ=255)
                                                   tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                                >  2:1437020/1‑91 (MQ=255)
                                                   tGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCt                                                                                >  2:91823/1‑139 (MQ=255)
                                                             tggtCATTGTCAGGGTTATCTCTTCCGCCGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCg                                                                                                          <  2:237148/103‑1 (MQ=255)
                                                             tggtCATTGTCAGGGTTATCTCTTCCGCCGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCg                                                                                                          >  1:237148/1‑103 (MQ=255)
                                                               gtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCg                                                                    >  2:526750/1‑139 (MQ=255)
                                                               gtCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCg                                                                    <  1:526750/139‑1 (MQ=255)
                                                                  aTTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCtt                                                                               <  1:98463/125‑1 (MQ=255)
                                                                  aTTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCtt                                                                               >  2:98463/1‑125 (MQ=255)
                                                                              aTCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCtt                                                                               <  2:434658/113‑1 (MQ=255)
                                                                              aTCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCtt                                                                               >  1:434658/1‑113 (MQ=255)
                                                                                    tCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                                <  1:1053559/58‑1 (MQ=255)
                                                                                    tCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                                >  2:1053559/1‑58 (MQ=255)
                                                                                    tCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGc                                                                                                           >  1:1050451/1‑79 (MQ=255)
                                                                                    tCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGc                                                                                                           <  2:1050451/79‑1 (MQ=255)
                                                                                    tCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATc                                                                                             >  1:107273/1‑93 (MQ=255)
                                                                                    tCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATc                                                                                             <  2:107273/93‑1 (MQ=255)
                                                                                    tCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCg                                               >  2:1113838/1‑139 (MQ=255)
                                                                                     ccGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCgg                                              >  2:1430439/1‑139 (MQ=255)
                                                                                      cGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGc                                             <  1:1665505/139‑1 (MQ=255)
                                                                                           tctcCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATc                                                                                             >  1:1771877/1‑86 (MQ=255)
                                                                                           tctcCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATc                                                                                             <  2:1771877/86‑1 (MQ=255)
                                                                                           tctcCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTAt                                        >  2:1700572/1‑139 (MQ=255)
                                                                                             tcCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                                <  1:6390/49‑1 (MQ=255)
                                                                                             tcCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTaa                                                                                                                                >  2:6390/1‑49 (MQ=255)
                                                                                             tcCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCa                                      <  1:1700572/139‑1 (MQ=255)
                                                                                             tcCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCa                                      >  2:1788416/1‑139 (MQ=255)
                                                                                              ccACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAg                                     <  1:1788416/139‑1 (MQ=255)
                                                                                                     ggCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCg                                                                                                          <  1:130398/63‑1 (MQ=255)
                                                                                                     ggCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCg                                                                                                          >  2:130398/1‑63 (MQ=255)
                                                                                                     ggCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGc                                                                                         <  2:841848/80‑1 (MQ=255)
                                                                                                     ggCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGc                                                                                         <  2:985251/80‑1 (MQ=255)
                                                                                                     ggCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGc                                                                                         >  1:985251/1‑80 (MQ=255)
                                                                                                     ggCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGc                                                                                         >  1:841848/1‑80 (MQ=255)
                                                                                                     ggCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCa                                                  >  2:1557952/1‑119 (MQ=255)
                                                                                                     ggCAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCa                                                  <  1:1557952/119‑1 (MQ=255)
                                                                                                       gaaGCTGAGACTGTGGGGTGATTTGTCGTAATTATTTAAAGAATGAGAACGGGCCGGCTCGTACTTCAGACATCCGGCCTTATTTCTTAAGCTTACCCGAGTGCCAATCCTGCGGCAGCGGCCTTATCAGCGATTATCa                            >  1:390540/2‑139 (MQ=255)
                                                                                                       cAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACa                                                                                                                    <  2:965505/51‑1 (MQ=255)
                                                                                                       cAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACa                                                                                                                    >  1:965505/1‑51 (MQ=255)
                                                                                                       cAAGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCa                            >  1:236987/1‑139 (MQ=255)
                                                                                                         aGCTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATc                          >  2:1380662/1‑139 (MQ=255)
                                                                                                           cTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCg                                                                                                               >  2:1109504/1‑52 (MQ=255)
                                                                                                           cTGAGACTGTGGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCg                                                                                                               <  1:1109504/52‑1 (MQ=255)
                                                                                                                    tgGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTAcc                             >  1:1054191/1‑125 (MQ=255)
                                                                                                                    tgGGGTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTAcc                             <  2:1054191/125‑1 (MQ=255)
                                                                                                                     ggggTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCg                                                                                                          <  1:128019/47‑1 (MQ=255)
                                                                                                                     ggggTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCg                                                                                                          >  2:128019/1‑47 (MQ=255)
                                                                                                                     ggggTGATTTTTCGTAATTATTTA‑‑‑AATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGc           >  1:1468849/1‑139 (MQ=255)
                                                                                                                        gTGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGc           <  1:1113838/139‑1 (MQ=255)
                                                                                                                         tGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTa                                                                                     >  1:894299/1‑64 (MQ=255)
                                                                                                                         tGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTa                                                                                     <  2:894299/64‑1 (MQ=255)
                                                                                                                         tGATTTTTCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCt          >  1:1192081/1‑139 (MQ=255)
                                                                                                                              tttCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCg                      <  1:1464822/122‑1 (MQ=255)
                                                                                                                              tttCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCg                      >  2:1464822/1‑122 (MQ=255)
                                                                                                                              tttCGTAATTATTTA‑‑‑AATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACg  <  2:1468849/139‑1 (MQ=255)
                                                                                                                               ttCGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAaca    <  1:1430439/139‑1 (MQ=255)
                                                                                                                                 cGTAATTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACg  <  2:853485/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                      
AAGAAGATTGTGAGAATTTCGGTATTAGAAATTCTGATCTAGCTAAATGCCTGCAAAAGGGTGGTCATTGTCAGGGTTATCTCTTCCGCTGTCTCCACTGCGGCAAGCTGAGACTGTGGGGTGATTTTTCGTAGTTATTTAAATAATGAGAACAGGCCGGAGCGTAATTCACACATCCGGCCTTATTTCTTAAGCTTAACCGAGTGCTAATTCTGCGGCAGCGGCTTTATCAGCGATTACCATCAGCGAAGGATGCAGCTGCAACACG  >  NC_000913/3898476‑3898743

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: