BRESEQ :: Evidence

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_000913	2,246,800	0	A	G	100.0%	350.3 / ‑6.7	105	intergenic (‑31/+263)	cirA/lysP	colicin 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 (60/44); total (60/45)

Only 100 of 116 total aligned reads displayed.

GGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCC  >  NC_000913/2246705‑2246895

ggCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaatgaaat                                                                                             >  2:1727825/1‑100 (MQ=255)

ggCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaatgaaat                                                                                             >  1:4951/1‑100 (MQ=255)

ggCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaatgaaat                                                                                             >  2:2978225/1‑100 (MQ=255)

ggCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaataaaat                                                                                             >  2:1967538/1‑100 (MQ=255)

ggCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGGAACTACGaaaatgaaat                                                                                             >  2:1194446/1‑100 (MQ=255)

 ggcATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAaaatgaaatg                                                                                            <  1:1558047/98‑1 (MQ=255)

 gCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGaaaatgaaat                                                                                             <  1:2901178/99‑1 (MQ=255)

     tGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGtt                                                                                        <  1:703210/100‑1 (MQ=255)

       cacaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTAt                                                                                      >  1:1735736/1‑100 (MQ=255)

        acaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTAt                                                                                      >  1:2920265/1‑99 (MQ=255)

        acaAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATc                                                                                     <  2:2468740/100‑1 (MQ=255)

         caAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCg                                                                                    >  2:19814/1‑100 (MQ=255)

         caAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCg                                                                                    >  2:2400567/1‑100 (MQ=255)

         caAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCg                                                                                    >  2:254032/1‑100 (MQ=255)

           aGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCt                                                                                  <  2:1035525/100‑1 (MQ=255)

            gAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTc                                                                                 <  2:757797/100‑1 (MQ=255)

             aaaTAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTca                                                                                >  2:111766/1‑100 (MQ=255)

                 aGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGATAAATACGAAAATGAAATGGATTAACGCTCAAAtc                                                                            >  1:1003774/1‑100 (MQ=255)

                    ggACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtcttc                                                                         >  1:1108357/1‑100 (MQ=255)

                    ggACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtctt                                                                          >  2:1212591/1‑99 (MQ=255)

                     gACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACTAAAATGAAATGGGTTATCGCTCACAtcttct                                                                        <  2:2709791/100‑1 (MQ=255)

                     gACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtcttct                                                                        >  2:56941/1‑100 (MQ=255)

                     gACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACAtcttct                                                                        <  1:2485917/100‑1 (MQ=255)

                       cAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCttcttc                                                                      <  2:2224049/100‑1 (MQ=255)

                         aaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCAc                                                                    >  1:16198/1‑100 (MQ=255)

                         aaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCAc                                                                    >  1:1095446/1‑100 (MQ=255)

                         aaCACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAAGTAACTACAAAAATGAAATGGGTTATCGCTCACATCTTCTTCAc                                                                    >  1:2314771/1‑100 (MQ=255)

                           cacaGCCCGACCCGTACGAAAGGGTTAAACCTAAACATTCCATATCTCAATGAGGAAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCaag                                                                   >  2:1907935/1‑97 (MQ=255)

                            acaGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTc                                                                 <  1:2773264/100‑1 (MQ=255)

                               gCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTcccc                                                              >  2:2443952/1‑100 (MQ=255)

                               gCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTcccc                                                              >  2:2508357/1‑100 (MQ=255)

                               gCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTccc                                                               >  2:2858197/1‑99 (MQ=255)

                               gCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTcc                                                                <  1:2933904/98‑1 (MQ=255)

                                cccGACCCGTACGAAAGGGTTCAACTTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCt                                                             <  1:494581/100‑1 (MQ=255)

                                 ccGACCCGTACGAAAGGGTTCAACCTAAACATTACATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCtt                                                            >  2:640324/1‑100 (MQ=255)

                                  cGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGATATCGCTCACATCTTCTTCACGTCCCCTTg                                                           >  1:2614603/1‑100 (MQ=255)

                                    aCCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCg                                                         <  1:2622623/100‑1 (MQ=255)

                                     cccGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGt                                                        >  2:2089414/1‑100 (MQ=255)

                                     cccGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGt                                                        >  1:2404196/1‑100 (MQ=255)

                                       cGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGGACACATATTGTTCACGTCACCTTCCGTAc                                                      >  2:2581777/1‑100 (MQ=255)

                                       cGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTAc                                                      >  1:1519011/1‑100 (MQ=255)

                                       cGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTAc                                                      >  2:801236/1‑100 (MQ=255)

                                        gTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACg                                                     <  1:1059593/100‑1 (MQ=255)

                                          aCGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgc                                                   <  2:1572993/100‑1 (MQ=255)

                                          aCGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgc                                                   >  1:1275570/1‑100 (MQ=255)

                                          aCGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgc                                                   <  1:1091384/100‑1 (MQ=255)

                                           cGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgca                                                  >  1:64537/1‑100 (MQ=255)

                                           cGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgca                                                  >  2:332538/1‑100 (MQ=255)

                                            gAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGgcag                                                 <  1:1407978/100‑1 (MQ=255)

                                            gAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATAGAAGATCACACACAGCTTCTTCACGTCCCCTTGCGTACGgc                                                   >  1:135003/1‑98 (MQ=255)

                                                gggTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATc                                             >  1:758021/1‑100 (MQ=255)

                                                gggTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATc                                             >  1:1007222/1‑100 (MQ=255)

                                                 gtttCAACCTAATCTTTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATc                                             <  1:1130127/97‑1 (MQ=255)

                                                 ggTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATcg                                            >  2:143346/1‑100 (MQ=255)

                                                   ttCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATcgcg                                          >  2:1202639/1‑100 (MQ=255)

                                                   ttCAACCTAAACATTCCATAGCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACTTCCCCTTGCGTACGGCAGCATcgcg                                          <  2:448794/100‑1 (MQ=255)

                                                     cAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCgg                                         >  2:1985527/1‑99 (MQ=255)

                                                     cAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGt                                        <  1:1727845/100‑1 (MQ=255)

                                                      aaCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCAAATCTTATTCACATCACCTTGCGTACGGCAGCAtagagg                                         >  1:2374874/1‑93 (MQ=255)

                                                       aCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTgg                                      <  2:2322749/100‑1 (MQ=255)

                                                       aCCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTg                                       >  1:1954874/1‑99 (MQ=255)

                                                         cTAAACATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGAc                                    >  1:1390684/1‑100 (MQ=255)

                                                             aCATTACATACCTCCGTGAGCTAACTACGAAAAGGCAGTCGCGTAGTGCTCACATCTTCTTCTAGTCGCGTTTCGTACTgctgtgtcgt                                           <  2:2552852/89‑9 (MQ=255)

                                                              cATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAtt                               >  2:1783170/1‑100 (MQ=255)

                                                              cATTCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAt                                <  2:178219/99‑1 (MQ=255)

                                                                ttCCATATCTCCATGAGGTAACTACGAAAATGAAATGTGTTATCGCTCACATCTTCTTCACGTCCCCTTACGTACGGCAGCATCGCGGTGGACTTAtttt                             >  2:2085606/1‑100 (MQ=255)

                                                                ttCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAtttt                             <  2:2050946/100‑1 (MQ=255)

                                                                 tCCATATCTCCATTAAGTAACTACTAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGAGTTAtttt                             <  2:1580608/99‑1 (MQ=255)

                                                                 tCCATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAttttt                            >  1:1779499/1‑100 (MQ=255)

                                                                  ccATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAtttttt                           <  1:1811167/100‑1 (MQ=255)

                                                                  ccATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTAtttttt                           >  2:2647806/1‑100 (MQ=255)

                                                                   cATATCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTa                          >  2:258104/1‑100 (MQ=255)

                                                                    atatCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTAt                         <  2:1564902/100‑1 (MQ=255)

                                                                    atatCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTAt                         <  1:1946914/100‑1 (MQ=255)

                                                                     tatCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTAt                         >  1:1463577/1‑99 (MQ=255)

                                                                     tatCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATg                        <  2:758017/100‑1 (MQ=255)

                                                                     tatCTCCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATg                        >  1:1566875/1‑100 (MQ=255)

                                                                       tctcCATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCa                      <  1:2853393/100‑1 (MQ=255)

                                                                          ccATGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGt                   <  1:565286/100‑1 (MQ=255)

                                                                             tGAGGTAACTACGAAAATGAAATGTGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGTACTTATTTTTTATGCAGGTGat                <  2:2002862/100‑1 (MQ=255)

                                                                             tGAGGTAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGat                <  2:1116681/100‑1 (MQ=255)

                                                                              gAGGTAACTACGAAAATGAAATGGGTTATGGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTTTGCAGGTGat                >  2:1972490/1‑99 (MQ=255)

                                                                               aGGTAACTACGCAAATGAAATGGTTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGa                 <  2:1500655/97‑1 (MQ=255)

                                                                                  tAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATcc           <  2:541839/100‑1 (MQ=255)

                                                                                  tAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATcc           <  2:4951/100‑1 (MQ=255)

                                                                                  tAACTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATcc           >  2:473008/1‑100 (MQ=255)

                                                                                     cTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCTCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAg        <  1:1243170/100‑1 (MQ=255)

                                                                                     cTACGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGCTCATccc          >  2:1730090/1‑98 (MQ=255)

                                                                                       aCGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTTCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaa      >  1:1275162/1‑100 (MQ=255)

                                                                                       aCGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaa      >  2:1170886/1‑100 (MQ=255)

                                                                                       aCGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGATGGACTTATTTTTTATGCAGGTGATCATCCCAGaa      >  1:1203379/1‑100 (MQ=255)

                                                                                        cGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACTGCAGCATCCCGGTGGACTTATTTTTTATGCAGGTGATCATACCAGaa      >  1:1302244/1‑99 (MQ=255)

                                                                                        cGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaaa     >  2:1728895/1‑100 (MQ=255)

                                                                                        cGAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGaaa     >  2:1069574/1‑100 (MQ=255)

                                                                                         gAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAg    >  2:878245/1‑100 (MQ=255)

                                                                                         gAAAATGAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAg    <  1:2858229/100‑1 (MQ=255)

                                                                                          aaagtaaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGTAGGTGATCATCCCAGAAAg    <  2:708315/99‑1 (MQ=255)

                                                                                          aaaatgaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGc   >  2:2342209/1‑100 (MQ=255)

                                                                                          aaaatgaaatGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGc   >  2:551826/1‑100 (MQ=255)

                                                                                           aaatgaaatgGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGcc  >  2:950055/1‑100 (MQ=255)

GGCCATGCACAAGAAATAGCGGACAAACACAGCCCGACCCGTACGAAAGGGTTCAACCTAAACATTCCATATCTCCATGAGGTAACTACGAAAATAAAATGGGTTATCGCTCACATCTTCTTCACGTCCCCTTGCGTACGGCAGCATCGCGGTGGACTTATTTTTTATGCAGGTGATCATCCCAGAAAGCC  >  NC_000913/2246705‑2246895

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