BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,542,218	A→C	6.5%	N61T (AAT→ACT)	fimE →	regulator for fimA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,542,218	0	A	C	6.5%	423.0 / 21.5	108	N61T (AAT→ACT)	fimE	regulator for fimA
Reads supporting (aligned to +/- strand): ref base A (53/48); new base C (3/4); total (56/52)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 120 total aligned reads displayed.

TTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGG  >  NC_000913/4542120‑4542315

ttattGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGaataaa                                                                                                   >  1:1055725/1‑99 (MQ=255)

ttattGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGaataaa                                                                                                   >  1:610873/1‑99 (MQ=255)

gacaaacatgatctcatatcagggacttgttcgcaccttccTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGaataaa                                                                                                   <  1:1775112‑M2/58‑1 (MQ=255)

  caaacatgatctcatatcagggacttgttcgcaccttccTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGaataaata                                                                                                 >  2:1782333‑M2/40‑99 (MQ=255)

  caaacatgatctcatatcagggacttgttcgcaccttccTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGaataaata                                                                                                 >  2:151312‑M2/40‑99 (MQ=255)

   ttGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAatat                                                                                                >  2:1425157/1‑99 (MQ=255)

     gTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTc                                                                                              >  1:1797992/1‑99 (MQ=255)

      tCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCg                                                                                             >  2:1124458/1‑99 (MQ=255)

       cTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGc                                                                                            <  1:775830/99‑1 (MQ=255)

       cTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGc                                                                                            <  1:1958222/99‑1 (MQ=255)

       cTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGc                                                                                            <  2:1563466/99‑1 (MQ=255)

       cTTATTCTGTTGGCACATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGc                                                                                            <  1:859793/99‑1 (MQ=255)

       atgatctcatatcagggacttgttcgcaccttccTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGc                                                                                            <  2:2455665‑M2/65‑1 (MQ=255)

        tgatctcatatcagggacttgttcgcaccttccTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGcc                                                                                           >  1:26598‑M2/34‑99 (MQ=255)

          aTTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGa                                                                                         <  1:1010517/99‑1 (MQ=255)

           ttCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGAc                                                                                        <  1:99948/99‑1 (MQ=255)

           tctcatatcagggacttgttcgcaccttccTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGAc                                                                                        >  2:2152993‑M2/31‑99 (MQ=255)

            tCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACt                                                                                       <  1:1056521/99‑1 (MQ=255)

            tCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACt                                                                                       >  1:2121343/1‑99 (MQ=255)

               gTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTgaa                                                                                    <  2:2328834/99‑1 (MQ=255)

               gTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGgg                                                                                    >  1:2117848‑M2/1‑97 (MQ=255)

                ttGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggt                                                                                   >  2:698380‑M2/1‑96 (MQ=255)

                  tcagggacttgttcgcaccttccTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTgaagaa                                                                                 <  1:2246107‑M2/76‑1 (MQ=255)

                   gCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAACTATTCGCCGACTGAAGAAc                                                                                <  2:313892/99‑1 (MQ=255)

                   gCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaat                                                                                >  1:162121‑M2/1‑93 (MQ=255)

                           tgttcgcaccttccTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTc                                                                        >  2:1329355‑M2/15‑99 (MQ=255)

                           tgttcgcaccttccTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTc                                                                        >  1:1875899‑M2/15‑99 (MQ=255)

                           gtaTGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGCACGGATTTTc                                                                        >  2:843308/3‑99 (MQ=255)

                            cATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAACTATTCGCCGACTGAAGAACGGATTTTCt                                                                       <  2:725856/99‑1 (MQ=255)

                            cATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaac                                                                       >  2:1693366‑M2/1‑84 (MQ=255)

                            cATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCt                                                                       >  1:1944505/1‑99 (MQ=255)

                              tGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAACTATTCGCCGACTGAAGAACGGATTTTCTAc                                                                     >  2:2261850/1‑99 (MQ=255)

                              tGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaactt                                                                     <  1:1158260‑M2/99‑18 (MQ=255)

                              tGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTAc                                                                     <  2:1869706/99‑1 (MQ=255)

                               gggATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaactta                                                                    >  2:405300‑M2/1‑81 (MQ=255)

                                  accttccTTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGtt                                                                 <  1:2112824‑M2/92‑1 (MQ=255)

                                   tGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattga                                                                <  1:282419‑M2/99‑23 (MQ=255)

                                   tGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattga                                                                >  2:215674‑M2/1‑77 (MQ=255)

                                    gCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCa                                                               <  2:637161/99‑1 (MQ=255)

                                     cGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgata                                                              <  1:1938198‑M2/99‑25 (MQ=255)

                                      gTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatag                                                             <  2:1812186‑M2/99‑26 (MQ=255)

                                         ttAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgt                                                          >  2:1935781‑M2/1‑71 (MQ=255)

                                           aGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTa                                                        <  1:1692197/99‑1 (MQ=255)

                                               aaCTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatg                                                    <  2:1057966‑M2/99‑35 (MQ=255)

                                               aaCTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGtt                                                    >  2:280338/1‑99 (MQ=255)

                                                aCTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGttt                                                   <  2:754627/99‑1 (MQ=255)

                                                   gCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTtga                                                <  1:1840672/99‑1 (MQ=255)

                                                   gCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTtga                                                >  2:1275397/1‑99 (MQ=255)

                                                    cTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcag                                               <  1:2046689‑M2/99‑40 (MQ=255)

                                                    cTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcag                                               >  2:1347485‑M2/1‑60 (MQ=255)

                                                     ttGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcaga                                              >  2:1666686‑M2/1‑59 (MQ=255)

                                                       gATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAACTATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAg                                            <  2:1214085/99‑1 (MQ=255)

                                                        aTCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataa                                           >  2:1453490‑M2/1‑56 (MQ=255)

                                                        aTCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataa                                           <  1:1341052‑M2/99‑44 (MQ=255)

                                                         tCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAACTATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCg                                          >  2:73237/1‑99 (MQ=255)

                                                         tCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCg                                          <  1:1836101/99‑1 (MQ=255)

                                                         tCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCg                                          >  2:344400/1‑99 (MQ=255)

                                                         tCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCg                                          >  1:301262/1‑99 (MQ=255)

                                                          cTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgatttatgttcagataatg                                         <  2:866591‑M2/99‑46 (MQ=255)

                                                           tGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgc                                        >  1:963656‑M2/1‑53 (MQ=255)

                                                            gCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttaagataatgac                                       >  2:1756111‑M2/1‑52 (MQ=255)

                                                             cATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgccc                                      <  2:2482968‑M2/99‑49 (MQ=255)

                                                              aTTATCAGGACCTTGACCTTAATGAAGGTAGAATAACTATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAg                                     <  1:100852/99‑1 (MQ=255)

                                                               ttATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGc                                    >  2:1210210/1‑99 (MQ=255)

                                                               ttATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGc                                    <  2:356397/99‑1 (MQ=255)

                                                                     ggACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGa                              >  2:2415900/1‑99 (MQ=255)

                                                                      gaccttgacctcAATGAAGGTAGACTGAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGaa                             >  2:2272390/1‑99 (MQ=255)

                                                                       accttgaccttAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttg                            <  1:475336‑M2/99‑59 (MQ=255)

                                                                          ttgaccttAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtca                         >  1:739489‑M2/1‑38 (MQ=255)

                                                                          ttgaccttAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtca                         >  2:1835786‑M2/1‑38 (MQ=255)

                                                                          ttgaccttAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtca                         >  2:166111‑M2/1‑38 (MQ=255)

                                                                          ttgaccttAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCt                         <  1:1730935/99‑1 (MQ=255)

                                                                            gaccttAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatg                       <  1:1226812‑M2/99‑64 (MQ=255)

                                                                            gaccttAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTgg                       <  2:2313660/99‑1 (MQ=255)

                                                                             accttAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgc                      <  2:1013228‑M2/99‑65 (MQ=255)

                                                                             accttAATGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgc                      <  2:1015763‑M2/99‑65 (MQ=255)

                                                                                ttAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGAccc                   >  1:1306077/1‑99 (MQ=255)

                                                                                ttAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGAccc                   <  1:2520010/99‑1 (MQ=255)

                                                                                  aaTGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctc                 >  1:1964818‑M2/1‑30 (MQ=255)

                                                                                    tGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctcca               >  1:1595773‑M2/1‑28 (MQ=255)

                                                                                    tGAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctcca               <  2:130325‑M2/99‑72 (MQ=255)

                                                                                     gAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccac              >  2:323454‑M2/1‑27 (MQ=255)

                                                                                     gAAGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccac              <  1:1420019‑M2/99‑73 (MQ=255)

                                                                                      aaGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccacc             >  2:863691‑M2/1‑26 (MQ=255)

                                                                                      aaGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccacc             <  2:68797‑M2/99‑74 (MQ=255)

                                                                                       aGGTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccg            <  1:803‑M2/99‑75 (MQ=255)

                                                                                       aGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACg            >  2:2042856/1‑99 (MQ=255)

                                                                                        ggTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccga           >  2:1231918‑M2/1‑24 (MQ=255)

                                                                                        ggTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccga           <  2:1659445‑M2/99‑76 (MQ=255)

                                                                                         gTAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgat          <  1:1913322‑M2/99‑77 (MQ=255)

                                                                                          tAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgatt         >  2:2229521‑M2/1‑22 (MQ=255)

                                                                                          tAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgatt         >  1:671070‑M2/1‑22 (MQ=255)

                                                                                          tAGAATAAATATTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgatt         >  1:2510714‑M2/1‑22 (MQ=255)

                                                                                          tAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGc         >  1:1699070/1‑99 (MQ=255)

                                                                                            gAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTa       >  2:196526/1‑99 (MQ=255)

                                                                                              ataaataTTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAAc     <  2:26598/99‑1 (MQ=255)

                                                                                               taaataTTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgag    <  2:1713114‑M2/99‑83 (MQ=255)

                                                                                                aaataTTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgaga   >  1:746557‑M2/1‑16 (MQ=255)

                                                                                                 aataTTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaa  <  1:1060028‑M2/99‑85 (MQ=255)

                                                                                                 aataTTCGCCGACTGggtaatgactccaacttattgatagtgttttatgttcagataatgcccgatgactttgtcatgcagctccaccgattttgagaa  <  2:975951‑M2/99‑85 (MQ=255)

TTATTGTCTTATTCTGTTGGCATATCGGCATGGGATGCGTATTAGTGAACTGCTTGATCTGCATTATCAGGACCTTGACCTTAATGAAGGTAGAATAAATATTCGCCGACTGAAGAACGGATTTTCTACCGTTCACCCGTTACGTTTTGATGAGCGTGAAGCCGTGGAACGCTGGACCCAGGAACGTGCTAACTGG  >  NC_000913/4542120‑4542315

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