Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,605,889 A→C 6.0% L135R (CTG→CGG)  PP_1407 ← sulfate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,605,8890AC 6.0% 190.1 / 4.6 67L135R (CTG→CGG) PP_1407sulfate transporter
Reads supporting (aligned to +/- strand):  ref base A (31/32);  new base C (2/2);  total (33/34)
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.90e-01

TGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGCCACGGCACCGGC  >  NC_002947/1605740‑1606037
                                                                                                                                                     |                                                                                                                                                    
tgatCAATGTGACCCAGTGAATTTCCCCCAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCa                                                                                                                                                      <  2:1068670/150‑1 (MQ=255)
 gatCAATGTGACCCAGTGAATTTCCCCCAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAg                                                                                                                                                     <  2:413335/150‑1 (MQ=255)
     aaTGTGACCCAGTGAATTTCCCCCAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATc                                                                                                                                                 <  1:1110104/150‑1 (MQ=255)
      aTGTGACCCAGTGAATTTCCCCCAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTACCCGATCAGGATCg                                                                                                                                                >  2:606693/1‑150 (MQ=255)
       tgtgACCCAGTGAATTTCCCCCAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCgg                                                                                                                                               <  1:973647/150‑1 (MQ=255)
             ccAGTGAATTTCCCCCAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTTGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCga                                                                                                                                         >  1:817269/1‑150 (MQ=255)
              cAGTGAATTTCCCCCAGGCGATGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCga                                                                                                                                         >  1:1156411/1‑149 (MQ=255)
              cAGTGAATTTCCCCCAGGCGATGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCga                                                                                                                                         <  2:1156411/149‑1 (MQ=255)
                   aaTTTCCCCCAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAgaagaa                                                                                                                                   >  2:842646/1‑150 (MQ=255)
                          cccAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCCTGATCGGctacct                                                                                                                                         <  2:957841/137‑6 (MQ=255)
                          cccAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCCTGATCGGctacct                                                                                                                                         >  1:957841/1‑132 (MQ=255)
                           ccAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGa                                                                                                                           <  2:75733/150‑1 (MQ=255)
                            cAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGAt                                                                                                                          <  2:1095014/150‑1 (MQ=255)
                                      aagaagTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGcc                                                                                                                >  1:70376/1‑150 (MQ=255)
                                           gTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCgcgc                                                                                                           >  2:1179034/1‑150 (MQ=255)
                                                    gCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCGGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCcgc                                                                                                  >  1:364532/1‑150 (MQ=255)
                                                     cTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCcgcg                                                                                                 >  2:483726/1‑150 (MQ=255)
                                                          gATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCGGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAAt                                                                                            <  2:364532/150‑1 (MQ=255)
                                                          gATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAAt                                                                                            >  2:948745/1‑150 (MQ=255)
                                                            tGAAACCGTCGCCCCCGATCTTGAAGCCCACCACCTTGGACAGCTGCTCGGCGATCAGTCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGACCTTCGCCCGCGCCAGGCCCGCGGCAATCa                                                                                          <  1:484046/150‑1 (MQ=255)
                                                             gAAACCGTCGCCCTCGATCTTGAAGCCCACCCCCTAGCACTGCTACCCGGCGATCAGGCTCAGGCCGATACCGTGGATGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAg                                                                                         <  1:606693/150‑1 (MQ=255)
                                                              aaaCCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGc                                                                                        <  1:1179034/150‑1 (MQ=255)
                                                                     cGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGc                                                                                 <  1:1256732/150‑1 (MQ=255)
                                                                         cTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGAc                                                                             >  1:1271036/1‑150 (MQ=255)
                                                                          tCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGAcc                                                                            >  2:671095/1‑150 (MQ=255)
                                                                             aTCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGAc                                                                             >  1:749661/1‑146 (MQ=255)
                                                                             aTCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGAc                                                                             <  2:749661/146‑1 (MQ=255)
                                                                              tCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTgg                                                                                                                              >  1:857278/1‑96 (MQ=255)
                                                                              tCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTgg                                                                                                                              <  2:857278/96‑1 (MQ=255)
                                                                              tCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCa                                                                        >  2:1008300/1‑150 (MQ=255)
                                                                                ttGAAGCCCACCACCTTGTACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACg                                                                      <  1:686912/150‑1 (MQ=255)
                                                                                ttGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACg                                                                      >  1:900402/1‑150 (MQ=255)
                                                                                   aaGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGCAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTg                                                                   >  1:830951/1‑150 (MQ=255)
                                                                                    aGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGa                                                                  <  2:817269/150‑1 (MQ=255)
                                                                                     gCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGAc                                                                 >  1:933946/1‑150 (MQ=255)
                                                                                     gCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGAc                                                                 >  2:816794/1‑150 (MQ=255)
                                                                                     gCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGAc                                                                 >  1:1270399/1‑150 (MQ=255)
                                                                                        caccacCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGAt                                                              >  1:688002/1‑150 (MQ=255)
                                                                                         accaccTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATc                                                             >  1:63384/1‑150 (MQ=255)
                                                                                           caccTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCAc                                                           >  2:1021257/1‑150 (MQ=255)
                                                                                            accTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCAcc                                                          <  1:1008300/150‑1 (MQ=255)
                                                                                              cTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGa                                                        <  2:900402/150‑1 (MQ=255)
                                                                                              cTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGa                                                        <  1:671095/150‑1 (MQ=255)
                                                                                              cTTGGACAGCTGCCCCGAGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGa                                                        <  2:615643/150‑1 (MQ=255)
                                                                                                        tGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTcgacg                                              >  2:1282722/1‑150 (MQ=255)
                                                                                                          cccGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGAt                                            <  2:456736/150‑1 (MQ=255)
                                                                                                           ccGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATg                                           <  1:483726/150‑1 (MQ=255)
                                                                                                                  tcaggctcaggcCGATACTGTTGAGGTAGCCGATAACCCTAGAATGACAGAAGAAGCTCACGATGATCCCAGTCGTCGCCAGTCCATTTGGAATTATCATCACGCCGTCCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCg                                    <  2:786029/150‑1 (MQ=255)
                                                                                                                      gctcaggcCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTc                                >  1:1211103/1‑150 (MQ=255)
                                                                                                                           ggcCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCa                           >  2:935250/1‑150 (MQ=255)
                                                                                                                                    cGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCg                  <  1:816794/150‑1 (MQ=255)
                                                                                                                                     gTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCgg                 >  1:843606/1‑150 (MQ=255)
                                                                                                                                      ttGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGt                >  1:100540/1‑150 (MQ=255)
                                                                                                                                        gAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGc              <  1:842646/150‑1 (MQ=255)
                                                                                                                                          ggTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATc                                                             >  1:925726/1‑101 (MQ=255)
                                                                                                                                          ggTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATc                                                             <  2:925726/101‑1 (MQ=255)
                                                                                                                                           gcggCCGATCCTGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTg                                                                   >  2:243397/4‑94 (MQ=255)
                                                                                                                                           gcggCCGATCCTGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTg                                                                   <  1:243397/91‑1 (MQ=255)
                                                                                                                                           gTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGCCAc           >  1:1080558/1‑150 (MQ=255)
                                                                                                                                           gTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGCCAc           <  2:1211103/150‑1 (MQ=255)
                                                                                                                                             aGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTcga                                                >  2:217696/1‑111 (MQ=255)
                                                                                                                                             aGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTcga                                                <  1:217696/111‑1 (MQ=255)
                                                                                                                                               ccGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGCCACGGCa       <  2:70376/150‑1 (MQ=255)
                                                                                                                                                  aTCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGCCACGGCACCg    >  2:664399/1‑150 (MQ=255)
                                                                                                                                                  aTCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGCCACGGCACCg    <  2:1271036/150‑1 (MQ=255)
                                                                                                                                                  aTCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGCCACGGCACCg    <  1:948745/150‑1 (MQ=255)
                                                                                                                                                   tCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCCGAAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGCCACGGCACCgg   <  1:13037/150‑1 (MQ=255)
                                                                                                                                                    cAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGCCACGGCACCGGc  <  2:1080558/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TGATCAATGTGACCCAGTGAATTTCCCCCAGGCGCTGGAAGAAGTTGATCAGGCTGAGGATGAAACCGTCGCCCTCGATCTTGAAGCCCACCACCTTGGACAGCTGCCCGGCGATCAGGCTCAGGCCGATACCGTTGAGGTAGCCGATCAGGATCGGCCGCGAGAAGAAGCTGGCGATGAACCCGGCCCGCGCCAGGCCCGCGGCAATCAGCATCACGCCGACCAGCACGGTGACGATCACCGACAGTTCGACGATGCGCTGCGGGTCACCCATGGCCAGCGGTGCCACGGCACCGGC  >  NC_002947/1605740‑1606037

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑