Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,385,199 A→G 6.5% *422Q (TAA→CAA)  oprD ← basic amino acid specific porin OprD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,385,1990AG 6.5% 174.7 / 2.3 61*422Q (TAA→CAA) oprDbasic amino acid specific porin OprD
Reads supporting (aligned to +/- strand):  ref base A (30/27);  new base G (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.36e-01

TTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGCTCAGGTCTTTCGCCGGGCCTTCCTGCAGCACGTACTTGATGT  >  NC_002947/1385052‑1385348
                                                                                                                                                   |                                                                                                                                                     
ttGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACa                                                                                                                                                     <  2:876934/150‑1 (MQ=255)
 tGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAg                                                                                                                                                    <  1:820040/150‑1 (MQ=255)
 tGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAg                                                                                                                                                    >  1:1065004/1‑150 (MQ=255)
  gCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAgg                                                                                                                                                   >  2:1166115/1‑150 (MQ=255)
    cTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGTAAGCCTGGCTTGCAgtgc                                                                                                                                                 <  2:1651873/150‑4 (MQ=255)
     tAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGGCGAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATg                                                                                                                                                >  1:1060567/1‑150 (MQ=255)
            ttGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCg                                                                                                                                         >  2:1828523/1‑150 (MQ=255)
                aaTGGCGCCGCGTTGGACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTa                                                                                                                                     >  2:1326520/1‑150 (MQ=255)
                 aTAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCAAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAAGTGCAGGGCACTGCACGCCAGGCTTACAGGATGCTCAGCGGGTAc                                                                                                                                    >  2:307057/1‑150 (MQ=255)
                     ccccGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGTAAGCCTGGCTTGCAgt                                                                                                                                                   <  1:558088/131‑2 (MQ=255)
                     ccccGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGTAAGCCTGGCTTGCAgt                                                                                                                                                   >  2:558088/1‑130 (MQ=255)
                      cccGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCAc                                                                                                                               >  1:1612617/1‑150 (MQ=255)
                          cATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATc                                                                                                                           <  2:1905309/150‑1 (MQ=255)
                          cATTGAACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATc                                                                                                                           <  1:1253597/150‑1 (MQ=255)
                               tACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGgcg                                                                                                                      <  2:1059954/150‑1 (MQ=255)
                                 cGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCa                                                                                                                    >  2:1612573/1‑150 (MQ=255)
                                    gggggAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGtt                                                                                                                 <  1:1166115/150‑1 (MQ=255)
                                          aaCGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGAATAACGATGAAAAGACCGCCATCATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGtcga                                                                                                           >  2:1499972/1‑150 (MQ=255)
                                            cGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGtcgatc                                                                                                         <  1:919194/150‑1 (MQ=255)
                                             gAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGCCGCCGTTGCCGGGATTTTGCTAATGCGTGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGtcgatcg                                                                                                        >  1:1177725/1‑150 (MQ=255)
                                                  cGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCACGCTTACAGGATGCTCAGTGGGTACGCCACGATGGGGCGCCGTTCGGCGATCGATggg                                                                                                   <  2:1455701/150‑1 (MQ=255)
                                                   gTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGTGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCGGTTCGTCGAGCGATgggg                                                                                                  >  1:175345/1‑150 (MQ=255)
                                                     aCAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACa                                                                                                                                                     >  2:1513505/1‑97 (MQ=255)
                                                     aCAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACa                                                                                                                                                     <  1:1513505/97‑1 (MQ=255)
                                                              cccAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGtagtag                                                                                       >  1:1222251/1‑150 (MQ=255)
                                                               ccAGGGCGTGACACTTGTCAGCCAAGAAAACGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGa                                                                                      >  2:1116174/1‑150 (MQ=255)
                                                                cAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGac                                                                                     <  1:777414/150‑1 (MQ=255)
                                                                    gCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCa                                                                                 <  1:941859/150‑1 (MQ=255)
                                                                      gTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATc                                                                               <  2:1222251/150‑1 (MQ=255)
                                                                         acacTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCgtag                                                                                          >  1:1729855/1‑136 (MQ=255)
                                                                         acacTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCgtag                                                                                          <  2:1729855/136‑1 (MQ=255)
                                                                           acTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGgaag                                                                          >  1:1305733/1‑150 (MQ=255)
                                                                            cTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGgaaga                                                                         >  2:1697240/1‑150 (MQ=255)
                                                                              tGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAAc                                                                       <  1:454385/150‑1 (MQ=255)
                                                                                tCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACgg                                                                     <  1:1697240/150‑1 (MQ=255)
                                                                                   gCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAg                                                                  >  2:117826/1‑150 (MQ=255)
                                                                                    ccAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCa                                                                                                                          >  2:133909/1‑93 (MQ=255)
                                                                                    ccAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCa                                                                                                                          <  1:133909/93‑1 (MQ=255)
                                                                                         aaaaaGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGGTGCTCAGAGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGATGACCGGTAGGTGGcc                                                            >  2:1054217/1‑150 (MQ=255)
                                                                                         aaaaaGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGcc                                                            <  2:1065004/150‑1 (MQ=255)
                                                                                          aaaaGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCt                                                           <  2:1060567/150‑1 (MQ=255)
                                                                                          aaaaGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCt                                                           <  1:530643/150‑1 (MQ=255)
                                                                                           aaaGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTg                                                          >  1:1710846/1‑150 (MQ=255)
                                                                                              gCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGacg                                                       >  2:933234/1‑150 (MQ=255)
                                                                                               cccGGCGTGATTGCCGGGCTTTTTCAGGTCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGacgca                                                     >  2:202219/1‑150 (MQ=255)
                                                                                                             cGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGCTCAg                                        <  1:1379792/150‑1 (MQ=255)
                                                                                                             cGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGCTCAg                                        >  2:349825/1‑150 (MQ=255)
                                                                                                                         ggTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGTCTGACGCACGCGGAAGCTCAGGTCTTTCGCCgg                            <  2:146569/150‑1 (MQ=255)
                                                                                                                            gCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGCTCAGGTCTTTCGCCGGGcc                         <  1:1116174/150‑1 (MQ=255)
                                                                                                                                 atcctgtaaGCCTGGCTTGCAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGCTCAGGTCTTTCGCCTGGCCTTCCt                    >  2:586727/8‑150 (MQ=255)
                                                                                                                                  cACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGCTCAGGTCTTTCGCCGGGCCTTCCTg                   <  2:1612617/150‑1 (MQ=255)
                                                                                                                                       cAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGCTCAGGTCTTTCGCCGGGCCTTCCTgcagca              >  2:590169/1‑150 (MQ=255)
                                                                                                                                            cAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGa                                                                                      <  1:152818/73‑1 (MQ=255)
                                                                                                                                            cAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGa                                                                                      >  2:152818/1‑73 (MQ=255)
                                                                                                                                            cAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGCTCAGGTCTTTCGCCgg                            >  1:1905575/1‑131 (MQ=255)
                                                                                                                                            cAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGCTCAGGTCTTTCGCCgg                            <  2:1905575/131‑1 (MQ=255)
                                                                                                                                             aGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGaca                                                                                    <  2:1060834/74‑1 (MQ=255)
                                                                                                                                             aGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGaca                                                                                    >  1:1060834/1‑74 (MQ=255)
                                                                                                                                               gCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGCTCAGGTCTTTCGCCGGGCCTTCCTGCAGCACGTACTTg      >  2:1897808/1‑150 (MQ=255)
                                                                                                                                                 ttACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGCTCAGGTCTTTCGCCGGGCCTTCCTGCAGCACGTACTTGAt    >  2:1035304/1‑150 (MQ=255)
                                                                                                                                                 ttACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTACGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGCTCAGGTCTTTCTCCGGGCCTTCCTGCAGCACGTACTTGAt    <  2:1177725/150‑1 (MQ=255)
                                                                                                                                                   aCAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGCTCAGGTCTTTCGCCGGGCCTTCCTGCAGCACGTACTTGATGt  <  2:1710846/150‑1 (MQ=255)
                                                                                                                                                   aCAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGCTCAGGTCTTTCGCCGGGCCTTCCTGCAGCACGTACTTGATGt  >  2:1311800/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
TTGCCTAAGACGTTGAAATAGCCCCGCATTGTACGGGGGGGAAACGAAGGCGTACAGGATGCCCCAGGGCGTGACACTTGTCAGCCAAGAAAAAGCCCGGCGTGATTGCCGGGCTTTTTCAGGTGCGGGGCACTGCAAGCCAGGCTTACAGGATGCTCAGCGGGTACTCCACGATCAGGCGCAGTTCGTCGATCGATGGGGAATCGTAGTAGACACCATCGGAAGAACGGTAGGTGGCCTGACGCACGCGGAAGCTCAGGTCTTTCGCCGGGCCTTCCTGCAGCACGTACTTGATGT  >  NC_002947/1385052‑1385348

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