Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,246,800 A→G 100% intergenic (‑31/+263) cirA ← / ← lysP colicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,246,8000AG98.8% 291.0 / ‑4.5 84intergenic (‑31/+263)cirA/lysPcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (47/36);  total (47/37)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

AGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGC  >  NC_000913/2246704‑2246894
                                                                                                |                                                                                              
aGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaatgaaa                                                                                             >  1:1781896/1‑100 (MQ=255)
 ggCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCCCCATGAGGTAACTACGaaaatgaaat                                                                                            <  2:822692/100‑1 (MQ=255)
 ggCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAAAATTCCATATCTCCCTGAGGTAACTACGaaaatgaaa                                                                                             >  2:2464928/1‑99 (MQ=255)
    cATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATggg                                                                                         <  1:441164/100‑1 (MQ=255)
    cATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATggg                                                                                         <  1:443443/100‑1 (MQ=255)
    cATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATggg                                                                                         <  2:689482/100‑1 (MQ=255)
    cATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATggg                                                                                         <  1:8096/100‑1 (MQ=255)
     aagCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGt                                                                                        >  1:1729536/3‑100 (MQ=255)
     aTGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGt                                                                                        <  1:861358/100‑1 (MQ=255)
     aTGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGt                                                                                        <  1:2025824/100‑1 (MQ=255)
     aTGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGt                                                                                        <  1:2157277/100‑1 (MQ=255)
       gCACAAGAAATAGCGGACAAACACCGCCCGTCCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGtt                                                                                       <  2:1561001/99‑1 (MQ=255)
        cacaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTAt                                                                                     >  1:1312595/1‑100 (MQ=255)
        cacaAGAAACAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTAt                                                                                     <  2:1657824/100‑1 (MQ=255)
         acaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATc                                                                                    <  2:2210973/100‑1 (MQ=255)
          caAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCg                                                                                   >  1:237898/1‑100 (MQ=255)
           aagaaagaGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGc                                                                                  <  1:151150/100‑1 (MQ=255)
             gAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTc                                                                                >  2:2329947/1‑100 (MQ=255)
             aaaaTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTc                                                                                >  2:1654731/2‑100 (MQ=255)
                aTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTcaca                                                                             <  2:1008047/100‑1 (MQ=255)
                  aGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtc                                                                           >  1:747252/1‑100 (MQ=255)
                     ggACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtcttc                                                                        >  2:2092305/1‑100 (MQ=255)
                     ggACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtcttc                                                                        >  1:1755591/1‑100 (MQ=255)
                        cAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCttcttc                                                                     >  2:131956/1‑100 (MQ=255)
                        cAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCttcttc                                                                     >  2:2517832/1‑100 (MQ=255)
                         aaaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCa                                                                    <  2:1729259/100‑1 (MQ=255)
                         aaaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCa                                                                    >  1:1616273/1‑100 (MQ=255)
                          gacacAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCAc                                                                   >  2:533635/2‑100 (MQ=255)
                           acacAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACg                                                                  <  1:652059/100‑1 (MQ=255)
                             acaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTc                                                                <  1:1617645/100‑1 (MQ=255)
                                gCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTcccc                                                             >  1:2272793/1‑100 (MQ=255)
                                gCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTcccc                                                             >  1:1411923/1‑100 (MQ=255)
                                gCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTcccc                                                             >  2:1978354/1‑100 (MQ=255)
                                gCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTcccc                                                             >  1:2014546/1‑100 (MQ=255)
                                 cccGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCt                                                            >  2:1855420/1‑100 (MQ=255)
                                 cccGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCt                                                            >  2:94906/1‑100 (MQ=255)
                                   cGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTg                                                          <  2:1962534/100‑1 (MQ=255)
                                    gACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTGACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGc                                                         <  2:2192821/100‑1 (MQ=255)
                                      cccGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGt                                                       >  2:738638/1‑100 (MQ=255)
                                      cccGTACGAAAGGGTTAAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACACCTCCTTCACGTCCCCTTg                                                          <  2:1476522/97‑1 (MQ=255)
                                        cGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTAc                                                     >  2:449596/1‑100 (MQ=255)
                                         gTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACg                                                    <  2:367858/100‑1 (MQ=255)
                                         gTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACg                                                    <  1:2162395/100‑1 (MQ=255)
                                          tACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACg                                                    >  2:812964/1‑99 (MQ=255)
                                               aaGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgcagca                                              <  2:1075632/100‑1 (MQ=255)
                                                  ggTTCAACCTAAACATTCCATATATCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCATATTCAAGTAACCTTGCGTACGGCAGCATcg                                           >  2:2468491/1‑100 (MQ=255)
                                                   gTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATcgc                                          >  2:1866520/1‑100 (MQ=255)
                                                   gTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATATTCTTCACGTCCCCTTGCGTACGGCAGCATcg                                           >  2:2266795/1‑99 (MQ=255)
                                                    ttCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATcgcg                                         >  2:2330009/1‑100 (MQ=255)
                                                     tCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCgg                                        <  2:2231598/100‑1 (MQ=255)
                                                     tCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCgg                                        <  2:1207646/100‑1 (MQ=255)
                                                     tCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCgg                                        <  1:131958/100‑1 (MQ=255)
                                                        aCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTgg                                     <  1:1974526/100‑1 (MQ=255)
                                                        aCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTgg                                     <  2:1188565/100‑1 (MQ=255)
                                                        aCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTgg                                     >  1:2543723/1‑100 (MQ=255)
                                                         ccTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGa                                    >  1:158534/1‑100 (MQ=255)
                                                         ccTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGa                                    >  1:1612823/1‑100 (MQ=255)
                                                         ccTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGa                                    >  1:1613470/1‑100 (MQ=255)
                                                          cTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGAc                                   >  2:229676/1‑100 (MQ=255)
                                                           tAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACt                                  >  2:2268643/1‑100 (MQ=255)
                                                           tAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACt                                  <  1:2462047/100‑1 (MQ=255)
                                                             aaCATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTa                                >  1:2166215/1‑100 (MQ=255)
                                                             aaCATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGCCCCCTTGCGTACGGCAGCATCGCGGTGGACTTa                                >  2:1244205/1‑100 (MQ=255)
                                                              aCATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAt                               >  1:2449693/1‑100 (MQ=255)
                                                               cATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAtt                              >  2:1879532/1‑100 (MQ=255)
                                                               cATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAtt                              >  2:1359022/1‑100 (MQ=255)
                                                               cATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAtt                              <  2:158532/100‑1 (MQ=255)
                                                                 ttCCATATCTCCATGAGGTAACTACGACAATGAAATGGGTTAGCGCTCACATCTTCTTCACGTACCCTTGCGTACGGCAGCATCGCGGTGGACTTAtttt                            <  2:1247691/100‑1 (MQ=255)
                                                                  tCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAttttt                           >  1:639647/1‑100 (MQ=255)
                                                                    cATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTa                         >  1:237433/1‑100 (MQ=255)
                                                                    cATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTa                         <  2:776382/100‑1 (MQ=255)
                                                                     atatCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTAt                        <  1:1425530/100‑1 (MQ=255)
                                                                         ctcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAg                    >  1:2179734/1‑100 (MQ=255)
                                                                          tcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAgg                   <  2:591231/100‑1 (MQ=255)
                                                                              tGAGGTGCCTACGAATGTGAGATGGGTTATTGCGGACTTCTTCTTCACGTCCCGTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTGAt                        <  2:2549118/91‑1 (MQ=255)
                                                                               gAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGat               <  2:1747878/99‑1 (MQ=255)
                                                                                  gTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGatcatc           >  2:856108/1‑100 (MQ=255)
                                                                                  gTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTGTATGCAGGTGATCATcc          <  1:578500/100‑1 (MQ=255)
                                                                                   tAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATcc          >  1:1784356/1‑100 (MQ=255)
                                                                                    aaCTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATccc         <  2:703969/100‑1 (MQ=255)
                                                                                    aaCTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATccc         >  1:510711/1‑100 (MQ=255)
                                                                                    aaCTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATccc         >  1:805347/1‑100 (MQ=255)
                                                                                     attaataaaatgaaatGGGTTATCGCTCACATCTTCTTCACGTCTCCATGCGTACGGAAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCa        <  2:1212864/94‑1 (MQ=255)
                                                                                        aCGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCATaa     >  1:205301/1‑100 (MQ=255)
                                                                                        aCGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaa     >  2:130812/1‑100 (MQ=255)
                                                                                         cGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaaa    >  1:232623/1‑100 (MQ=255)
                                                                                          gAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAg   <  1:1836932/100‑1 (MQ=255)
                                                                                          gAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAg   >  2:1043243/1‑100 (MQ=255)
                                                                                           aaaatgaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGc  <  2:598571/100‑1 (MQ=255)
                                                                                           aaaatgaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGc  <  1:1373809/100‑1 (MQ=255)
                                                                                                |                                                                                              
AGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGC  >  NC_000913/2246704‑2246894

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