BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,063,931	A→G	V122A (GTA→GCA)	ompL ←	putative outer membrane porin L

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,063,931	0	A	G	100.0%	1417.9 / NA	337	V122A (GTA→GCA)	ompL	putative outer membrane porin L
Reads supporting (aligned to +/- strand): ref base A (0/0); new base G (175/162); total (175/162)

Only 100 of 337 total aligned reads displayed.

TAAAATTCCAGTAGGTGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTACTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGCCACCGGGCTG  >  NC_000913/4063833‑4064028

 aaaaTTCCAGTAGGTGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTg                                                                                                >  1:6695842/1‑101 (MQ=255)

  aaaTTCCAGTAGGTGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGt                                                                                               >  1:6695846/1‑101 (MQ=255)

  aaaTTCCAGTAGGTGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGt                                                                                               >  1:4101571/1‑101 (MQ=255)

    aTTCCGGTAGGTGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCa                                                                                             >  2:6695849/1‑101 (MQ=255)

     ttCCAGTAGGTGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAg                                                                                            >  1:6695851/1‑101 (MQ=255)

     ttCCAGTAGGTGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAg                                                                                            <  1:6695505/101‑1 (MQ=255)

     ttCCAGTAGGTGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAg                                                                                            <  2:6695330/101‑1 (MQ=255)

      tCCAGTAGGTGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGa                                                                                           <  1:6695109/101‑1 (MQ=255)

      tCCAGTAGGTGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGa                                                                                           <  1:6695763/101‑1 (MQ=255)

       ccAGTAGGTGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGAt                                                                                          >  2:6695855/1‑101 (MQ=255)

        cAGTAGGTGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGAtt                                                                                         <  2:6695483/101‑1 (MQ=255)

         aGTAGGTGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTa                                                                                        >  1:6695859/1‑101 (MQ=255)

           tAGGTGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTaaa                                                                                      >  2:6695860/1‑101 (MQ=255)

             ggTGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAAcc                                                                                    >  2:6695865/1‑101 (MQ=255)

               tGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCAt                                                                                  <  1:6695317/101‑1 (MQ=255)

                gCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATg                                                                                 <  2:6695027/101‑1 (MQ=255)

                gCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATg                                                                                 >  1:6695869/1‑101 (MQ=255)

                gCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATg                                                                                 <  1:6695182/101‑1 (MQ=255)

                gCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATg                                                                                 <  2:6695304/101‑1 (MQ=255)

                 ccAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATgg                                                                                >  2:6695872/1‑101 (MQ=255)

                   aaTTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTa                                                                              >  2:6695878/1‑101 (MQ=255)

                   aaTTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTa                                                                              <  2:6695444/101‑1 (MQ=255)

                   aaTTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTa                                                                              >  1:6695876/1‑101 (MQ=255)

                    aTTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTAc                                                                             >  1:6695880/1‑101 (MQ=255)

                      ttCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACaa                                                                           >  1:6695885/1‑101 (MQ=255)

                       tCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACaaa                                                                          <  2:6695087/101‑1 (MQ=255)

                        cATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAAt                                                                         <  1:6695299/101‑1 (MQ=255)

                         aTAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAAtt                                                                        >  2:6695888/1‑101 (MQ=255)

                          tAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATTGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAAttt                                                                       >  1:6695892/1‑101 (MQ=255)

                          tAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAAttt                                                                       >  2:6695893/1‑101 (MQ=255)

                             gtgtCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTtata                                                                    >  2:6695898/1‑101 (MQ=255)

                             gtgtCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTtata                                                                    <  2:6695393/101‑1 (MQ=255)

                               gtCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGt                                                                  >  1:6695899/1‑101 (MQ=255)

                                  attattATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGCTTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGa                                                               >  2:6695905/1‑101 (MQ=255)

                                  attattATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGa                                                               <  2:6695798/101‑1 (MQ=255)

                                  attattATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGa                                                               <  2:6695253/101‑1 (MQ=255)

                                   ttattaTCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGAc                                                              >  1:6695906/1‑101 (MQ=255)

                                      ttatCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTc                                                           >  2:6695911/1‑101 (MQ=255)

                                        atCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCa                                                         <  1:6695522/101‑1 (MQ=255)

                                        atCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCa                                                         >  1:6695914/1‑101 (MQ=255)

                                           cAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGAt                                                      <  1:6695559/101‑1 (MQ=255)

                                            aGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATa                                                     >  1:6695918/1‑101 (MQ=255)

                                            aGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATa                                                     <  2:6695495/101‑1 (MQ=255)

                                             gTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATac                                                    <  1:6695286/101‑1 (MQ=255)

                                              ttCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATaca                                                   <  2:6695404/101‑1 (MQ=255)

                                              ttCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATaca                                                   >  2:6695920/1‑101 (MQ=255)

                                               tCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATacac                                                  >  1:6695923/1‑101 (MQ=255)

                                                ccccGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAAc                                                                                     <  1:6695305/65‑1 (MQ=255)

                                                ccccGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACt                                                 <  1:6695668/101‑1 (MQ=255)

                                                 cccGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTg                                                <  2:6694849/101‑1 (MQ=255)

                                                  ccGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTtata                                                                    <  1:6695369/80‑1 (MQ=255)

                                                   cGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGca                                              <  1:6695267/101‑1 (MQ=255)

                                                   cGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGca                                              >  1:6695928/1‑101 (MQ=255)

                                                   cGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGca                                              <  1:6695440/101‑1 (MQ=255)

                                                     cTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGcacc                                            <  1:6695329/101‑1 (MQ=255)

                                                     cTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGcacc                                            <  1:6695233/101‑1 (MQ=255)

                                                       cAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGcaccac                                          >  2:6695933/1‑101 (MQ=255)

                                                        aaaTCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCaccacc                                         >  2:6695937/1‑101 (MQ=255)

                                                        aaaTCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCaccacc                                         >  1:6695936/1‑101 (MQ=255)

                                                            cTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGaa                                     >  2:6695945/1‑101 (MQ=255)

                                                            cTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGaa                                     >  1:6695944/1‑101 (MQ=255)

                                                            cTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGaa                                     >  1:6695949/1‑101 (MQ=255)

                                                             tGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAAc                                    <  1:6695699/101‑1 (MQ=255)

                                                             tGTTGAACTATAGTTATTATGGTTATAGAGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAAc                                    <  1:6695251/101‑1 (MQ=255)

                                                              gTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCTCCACCCGAAcc                                   <  2:6695175/101‑1 (MQ=255)

                                                              gTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAAcc                                   <  2:6695353/101‑1 (MQ=255)

                                                              gTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAAcc                                   <  1:6695433/101‑1 (MQ=255)

                                                              gTTGAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAAcc                                   >  1:6695952/1‑101 (MQ=255)

                                                                 gAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGAt                                <  1:6695478/101‑1 (MQ=255)

                                                                 gAACTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGAt                                >  2:6695955/1‑101 (MQ=255)

                                                                   aCTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATAc                              >  2:6695959/1‑101 (MQ=255)

                                                                    cTATAGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATAct                             >  1:6695961/1‑101 (MQ=255)

                                                                      ataGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATActct                           >  1:6695967/1‑101 (MQ=255)

                                                                        aGTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTa                         <  2:6695612/101‑1 (MQ=255)

                                                                         gTTATTATGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTAt                        >  2:6695969/1‑101 (MQ=255)

                                                                            attatGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCat                     >  1:6695974/1‑101 (MQ=255)

                                                                            attatGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCat                     >  2:6695975/1‑101 (MQ=255)

                                                                             ttatGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCatt                    >  1:6695978/1‑101 (MQ=255)

                                                                              tatGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCatta                   <  2:6695368/101‑1 (MQ=255)

                                                                              tatGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCatta                   <  2:6695281/101‑1 (MQ=255)

                                                                                tGGTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCattaat                 <  1:6695597/101‑1 (MQ=255)

                                                                                 ggTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGAt                                <  2:6695379/85‑1 (MQ=255)

                                                                                  gTTATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCattaatta               >  1:6695983/1‑101 (MQ=255)

                                                                                   ttATAGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCAttaattaa              >  2:6695985/1‑101 (MQ=255)

                                                                                     ataGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGc            >  1:6695990/1‑101 (MQ=255)

                                                                                     ataGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGc            <  1:6695392/101‑1 (MQ=255)

                                                                                      taGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGcc           >  2:6695993/1‑101 (MQ=255)

                                                                                       aGCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGCCa          >  1:6695995/1‑101 (MQ=255)

                                                                                        gCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGCCAc         <  2:6695466/101‑1 (MQ=255)

                                                                                        gCGATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGCCAc         <  1:6695465/101‑1 (MQ=255)

                                                                                          gATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGCCACCg       >  1:6696000/1‑101 (MQ=255)

                                                                                          gATTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGCCACCg       >  2:6696002/1‑101 (MQ=255)

                                                                                           aTTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGCCACCgg      <  2:6695568/101‑1 (MQ=255)

                                                                                           aTTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGCCACCgg      <  2:6695425/101‑1 (MQ=255)

                                                                                           aTTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGCCACCgg      >  1:6696006/1‑101 (MQ=255)

                                                                                           aTTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGCCACCgg      <  1:6695401/101‑1 (MQ=255)

                                                                                           aTTTCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGAAACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGCCACCgg      <  2:6695658/101‑1 (MQ=255)

                                                                                             ttCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGCCACCGGGc    <  2:6695544/101‑1 (MQ=255)

                                                                                              tCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGCCACCGGGCt   <  2:6695344/101‑1 (MQ=255)

                                                                                              tCGTGCTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGCCACCGGGCt   >  1:6696010/1‑101 (MQ=255)

TAAAATTCCAGTAGGTGCCAATTTCATAGGTGTCATTATTATCCAGTTCCCCGCTCAAATCTGTTGAACTATAGTTATTATGGTTATAGCGATTTCGTACTGTCAGATTAAACCATGGTACAAATTTATAGTTGACGTCCAGATACACTGCACCACCCGAACCGATACTCTTATCATTAATTAAGCCACCGGGCTG  >  NC_000913/4063833‑4064028

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 33 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`