Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,246,800 A→G 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,8000AG94.4% 288.8 / 3.3 90intergenic (‑31/+263)cirA/lysPcolicin IA outer membrane receptor and translocator; ferric iron‑catecholate transporter/lysine transporter
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (47/37);  total (49/41)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

CAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGC  >  NC_000913/2246703‑2246894
                                                                                                 |                                                                                              
caGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaatgaa                                                                                              >  1:20144/1‑100 (MQ=255)
  ggCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaatgaaat                                                                                            <  1:1305440/100‑1 (MQ=255)
  cccatctcacaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaatgaaat                                                                                            <  1:568109/93‑1 (MQ=255)
    ccATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATgg                                                                                          <  1:976476/100‑1 (MQ=255)
    ccATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATgg                                                                                          >  2:1013190/1‑100 (MQ=255)
      aTGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGt                                                                                        <  1:926879/100‑1 (MQ=255)
      aTGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGt                                                                                        <  1:780547/100‑1 (MQ=255)
      aTGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGt                                                                                        <  1:1454441/100‑1 (MQ=255)
      aTGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGt                                                                                        <  1:1648229/100‑1 (MQ=255)
        gctcaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAGGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTGCGAAAATTAAATTGGt                                                                                        <  1:57775/95‑1 (MQ=255)
         cacaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaatgaaatagatcgg                                                                                     >  2:568111/1‑93 (MQ=255)
         cacaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTAt                                                                                     >  1:1377843/1‑100 (MQ=255)
         cacaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTAt                                                                                     >  2:328806/1‑100 (MQ=255)
          acaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATc                                                                                    >  1:195639/1‑100 (MQ=255)
           caAGAAATAGCGGACAAACACACCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATAAGGAAACTACGAAAATAAAATGAGTTAtag                                                                                   >  2:719662/1‑98 (MQ=255)
             aGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCt                                                                                 >  1:1289234/1‑100 (MQ=255)
               aaaTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTca                                                                               <  1:1375930/100‑1 (MQ=255)
               aaaTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTca                                                                               <  2:1777292/100‑1 (MQ=255)
               aaaTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTca                                                                               <  2:1320123/100‑1 (MQ=255)
               aaaGAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTca                                                                               <  1:513658/100‑1 (MQ=255)
                aaTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTcac                                                                              <  1:1418148/100‑1 (MQ=255)
                aaTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTcac                                                                              >  1:783131/1‑100 (MQ=255)
                 aTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTcaca                                                                             <  2:879130/100‑1 (MQ=255)
                 aTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTcaca                                                                             >  1:1900564/1‑100 (MQ=255)
                  tAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAt                                                                            >  2:1616188/1‑100 (MQ=255)
                  tAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCCCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAt                                                                            >  1:610455/1‑100 (MQ=255)
                   aGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtc                                                                           >  1:1065701/1‑100 (MQ=255)
                      ggACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtcttc                                                                        >  2:657188/1‑100 (MQ=255)
                       gACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtcttct                                                                       >  1:1627799/1‑100 (MQ=255)
                       gACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtcttct                                                                       >  2:1131545/1‑100 (MQ=255)
                        aCAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATcttctt                                                                      <  1:1869936/100‑1 (MQ=255)
                        aCAAACACACCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATcttctt                                                                      <  1:1613857/100‑1 (MQ=255)
                         cAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCttcttc                                                                     <  2:1289241/100‑1 (MQ=255)
                         cAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCttcttc                                                                     >  2:1700916/1‑100 (MQ=255)
                         cAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCttcttc                                                                     >  1:975752/1‑100 (MQ=255)
                           aaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCa                                                                    <  1:719661/99‑1 (MQ=255)
                           aaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCAc                                                                   <  2:324946/100‑1 (MQ=255)
                           aaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCAc                                                                   >  1:1536416/1‑100 (MQ=255)
                            acacAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACg                                                                  <  2:1876140/100‑1 (MQ=255)
                            acacAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACg                                                                  >  1:1711268/1‑100 (MQ=255)
                             cacaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGt                                                                 >  1:392911/1‑100 (MQ=255)
                             cacaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGt                                                                 >  2:58291/1‑100 (MQ=255)
                             cacaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGt                                                                 >  1:847602/1‑100 (MQ=255)
                             cacaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAACTTCTTCACGt                                                                 >  1:1615774/1‑100 (MQ=255)
                              acaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTc                                                                >  2:1266436/1‑100 (MQ=255)
                                 gCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTcccc                                                             <  2:978943/100‑1 (MQ=255)
                                 gCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTcccc                                                             <  1:1405584/100‑1 (MQ=255)
                                       cccGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGt                                                       >  1:91411/1‑100 (MQ=255)
                                       cccGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGt                                                       <  2:343842/100‑1 (MQ=255)
                                         cGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTAc                                                     >  1:595000/1‑100 (MQ=255)
                                         cGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTAc                                                     >  1:595070/1‑100 (MQ=255)
                                            aCGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgc                                                  <  1:863662/100‑1 (MQ=255)
                                            aCGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgc                                                  <  2:648023/100‑1 (MQ=255)
                                                  gggTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATc                                            >  2:795877/1‑100 (MQ=255)
                                                   ggTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATcg                                           >  2:1530478/1‑100 (MQ=255)
                                                    gTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATcgc                                          <  2:1015486/100‑1 (MQ=255)
                                                      tCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCgg                                        <  2:284038/100‑1 (MQ=255)
                                                       cAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGt                                       >  2:1789391/1‑100 (MQ=255)
                                                        aaCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTg                                      <  2:876574/100‑1 (MQ=255)
                                                         aCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTgg                                     <  1:314557/100‑1 (MQ=255)
                                                          ccTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGa                                    >  1:306727/1‑100 (MQ=255)
                                                          ccTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGa                                    >  1:287738/1‑100 (MQ=255)
                                                            tAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACt                                  <  1:498861/100‑1 (MQ=255)
                                                              aaCATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTa                                <  2:1536424/100‑1 (MQ=255)
                                                                 aTTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAttt                             <  2:364125/100‑1 (MQ=255)
                                                                 aTTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAttt                             >  2:1681585/1‑100 (MQ=255)
                                                                   tCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAttttt                           >  2:1799738/1‑100 (MQ=255)
                                                                     cATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTa                         >  2:1476097/1‑100 (MQ=255)
                                                                         tctcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCa                     >  1:217216/1‑100 (MQ=255)
                                                                         tctcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCa                     <  1:1931486/100‑1 (MQ=255)
                                                                         tctcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCa                     >  2:705035/1‑100 (MQ=255)
                                                                          ctcCATGAGGTAACTACGAAAATGAAATGGGTTATCACTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAg                    >  1:1605526/1‑100 (MQ=255)
                                                                           tcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAgg                   <  1:750159/100‑1 (MQ=255)
                                                                           tcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAgg                   <  2:617688/100‑1 (MQ=255)
                                                                           tcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAgg                   >  2:1485901/1‑100 (MQ=255)
                                                                             caggAGGTCACTACGAAAATGAAATGGGTTATCGCGAAAATCTTCTTAACGTCCCCGTCCCTACGGCACCATCGCGGTTGACTTATTTTCTATCCAGGta                 >  1:1771664/4‑99 (MQ=255)
                                                                             cATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTg                 <  2:607274/100‑1 (MQ=255)
                                                                               tGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGat               >  2:344740/1‑100 (MQ=255)
                                                                                ccgatctcTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGatc              <  1:115799/93‑1 (MQ=255)
                                                                                  ggTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGatcat            <  2:783132/100‑1 (MQ=255)
                                                                                   gTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGatcatc           >  2:576124/1‑100 (MQ=255)
                                                                                    tAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATcc          <  1:1396966/100‑1 (MQ=255)
                                                                                      aCTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCa        >  1:1874035/1‑100 (MQ=255)
                                                                                      aCTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCa        <  1:795875/100‑1 (MQ=255)
                                                                                       cTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGatcagatcgg       >  2:115799/1‑94 (MQ=255)
                                                                                       cTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAg       >  1:1383438/1‑100 (MQ=255)
                                                                                       cTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGAAGCATCGCGGAGGACTTATTTTTTATGCAGGTGATCATCCCa        >  2:1160033/1‑99 (MQ=255)
                                                                                       cTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTACCCTTGCGTACCGCAGCATTGTGGTGGACTTATTTTTTATGCAGGTGATCATCCCAg       >  1:313029/1‑100 (MQ=255)
                                                                                        tACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGa      >  2:289000/1‑100 (MQ=255)
                                                                                         aCGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaa     <  1:934623/100‑1 (MQ=255)
                                                                                          cGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaaa    <  2:595072/100‑1 (MQ=255)
                                                                                          cGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaaa    >  1:949675/1‑100 (MQ=255)
                                                                                          cGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaaa    <  2:595002/100‑1 (MQ=255)
                                                                                          cGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaaa    >  2:288859/1‑100 (MQ=255)
                                                                                            aaaatgaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGc  >  1:286921/1‑100 (MQ=255)
                                                                                                 |                                                                                              
CAGGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGC  >  NC_000913/2246703‑2246894

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: