BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,624,135	A→G	N100N (AAT→AAC)	b3485 ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,624,135	0	A	G	99.8%	1371.3 / ‑5.2	426	N100N (AAT→AAC)	b3485	putative transporter
Reads supporting (aligned to +/- strand): ref base A (0/1); new base G (198/227); total (198/228)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.96e-01

Only 100 of 438 total aligned reads displayed.

AAATGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATATTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGACGATA  >  NC_000913/3624045‑3624224

aaaTGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGtt                                                                                         <  2:197061/93‑1 (MQ=255)

 aaTGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGttt                                                                                        <  2:7366881/93‑1 (MQ=255)

  aTGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTa                                                                                       <  1:4586947/93‑1 (MQ=255)

  aTGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTa                                                                                       >  2:3977194/1‑93 (MQ=255)

  aTGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTa                                                                                       <  1:5767430/93‑1 (MQ=255)

   tGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTAt                                                                                      >  1:2525639/1‑93 (MQ=255)

   acccTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTAt                                                                                      <  2:6250955/91‑1 (MQ=255)

      cTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATcgc                                                                                   >  2:7671988/1‑93 (MQ=255)

       tGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATcgcg                                                                                  >  1:6063114/1‑93 (MQ=255)

         gCTCATGCGCGTGGCATCGACGTTCGCCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGaa                                                                                <  2:6420049/93‑1 (MQ=255)

         gCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGaa                                                                                <  1:7534080/93‑1 (MQ=255)

         gCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGaa                                                                                <  1:3502550/93‑1 (MQ=255)

         gCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGaa                                                                                <  1:3010963/93‑1 (MQ=255)

          cTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAg                                                                               >  1:5072722/1‑93 (MQ=255)

          cTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAg                                                                               <  2:4412839/93‑1 (MQ=255)

          cTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAg                                                                               <  1:7390876/93‑1 (MQ=255)

           tCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAgg                                                                              >  1:3902518/1‑93 (MQ=255)

           tCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAgg                                                                              >  2:3789418/1‑93 (MQ=255)

             aTGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGta                                                                            <  2:4654837/93‑1 (MQ=255)

             aTGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGta                                                                            <  1:5531385/93‑1 (MQ=255)

             aTGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGta                                                                            <  2:3342369/93‑1 (MQ=255)

               gcgcGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGtata                                                                          >  1:5004804/1‑93 (MQ=255)

                cgcgTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAg                                                                         <  1:6835063/93‑1 (MQ=255)

                cgcgTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAg                                                                         <  1:7122183/93‑1 (MQ=255)

                 gcgTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGc                                                                        <  1:3672936/93‑1 (MQ=255)

                 gcgTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGc                                                                        <  2:5874721/93‑1 (MQ=255)

                   gtagCATCGTCGATCACCTGAATATCCGGCTGGCGTCCGGCGAAGACATCACGCTGAAAATTAGGCAAATTGTTTATCGCGAAGGTATAGCga                                                                      <  2:7043049/90‑2 (MQ=255)

                   gTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGt                                                                      >  1:7900208/1‑93 (MQ=255)

                    tGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTc                                                                     <  1:3252692/93‑1 (MQ=255)

                    tGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTc                                                                     >  2:2434528/1‑93 (MQ=255)

                     ggCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTcc                                                                    <  1:3135789/93‑1 (MQ=255)

                      gCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCg                                                                   <  1:4679647/93‑1 (MQ=255)

                         tCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCg                                                                >  2:123216/1‑93 (MQ=255)

                         tCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCg                                                                >  2:7622477/1‑93 (MQ=255)

                         tCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCg                                                                >  2:2265578/1‑93 (MQ=255)

                         tCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCg                                                                >  1:2813549/1‑93 (MQ=255)

                         tCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCAGGCg                                                                >  1:1029819/1‑93 (MQ=255)

                         tCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTACGGct                                                                >  2:919506/1‑92 (MQ=255)

                           gACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTc                                                              >  1:2252484/1‑93 (MQ=255)

                             cGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCa                                                            <  2:3178183/93‑1 (MQ=255)

                              gTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAg                                                           >  2:1540964/1‑93 (MQ=255)

                              gTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAg                                                           <  2:5081541/93‑1 (MQ=255)

                                tCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTc                                                         <  2:5782694/93‑1 (MQ=255)

                                tCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGGATAGCGTCCGGCGTCCAGTc                                                         >  2:1288383/1‑93 (MQ=255)

                                 cACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTcc                                                        <  2:4840703/93‑1 (MQ=255)

                                 cACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTcc                                                        <  1:3283743/93‑1 (MQ=255)

                                 cACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTcc                                                        <  1:2182823/93‑1 (MQ=255)

                                       aaTATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGcatc                                                  <  2:6873831/93‑1 (MQ=255)

                                         tatCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGcatcca                                                <  2:2116461/93‑1 (MQ=255)

                                         tatCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGcatcca                                                >  1:544634/1‑93 (MQ=255)

                                           tCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGcatccatc                                              <  1:6717360/93‑1 (MQ=255)

                                           accGGCTGGCGTCCGGCGAGGACAGCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGcatccatc                                              <  2:5587336/92‑1 (MQ=255)

                                              ggCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCa                                           <  1:4243553/93‑1 (MQ=255)

                                               gCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCTGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCAt                                          <  2:3264918/93‑1 (MQ=255)

                                               gCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCAt                                          >  1:872216/1‑93 (MQ=255)

                                                  ggCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATcgg                                       >  1:4048009/1‑93 (MQ=255)

                                                    cGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATcggcg                                     <  1:3824210/93‑1 (MQ=255)

                                                    cGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATcggcg                                     >  1:2055075/1‑93 (MQ=255)

                                                     gTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATcggcgg                                    >  1:1668768/1‑93 (MQ=255)

                                                     gTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATcggcgg                                    >  1:1454214/1‑93 (MQ=255)

                                                      tCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGt                                   >  1:3328829/1‑93 (MQ=255)

                                                       ccGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTg                                  <  1:5239251/93‑1 (MQ=255)

                                                        cGGCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGa                                 <  2:4966546/93‑1 (MQ=255)

                                                          gCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGatc                               <  1:2586412/93‑1 (MQ=255)

                                                          gCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGatc                               <  2:452305/93‑1 (MQ=255)

                                                          gCGAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGatc                               <  1:6575895/93‑1 (MQ=255)

                                                           cGAGGACATCACGCTGAAAATTAGGCGGAATGGTTATCGCGAAGGTATAGCGTCCGGCGGCCAGACCGGCATCCATCGCATCGGCGGCGagca                              >  1:7029039/1‑93 (MQ=255)

                                                            gacgACGTCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGatcat                             <  2:2532561/90‑1 (MQ=255)

                                                            gAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGatcat                             >  1:5201795/1‑93 (MQ=255)

                                                            gAGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGCATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGagcat                             >  1:4877193/1‑93 (MQ=255)

                                                             aGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGatcatc                            >  1:2888820/1‑93 (MQ=255)

                                                             aGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGatcatc                            >  1:992075/1‑93 (MQ=255)

                                                             aGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGatcatc                            >  2:6628656/1‑93 (MQ=255)

                                                             aGGACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGatcatc                            <  2:1143880/93‑1 (MQ=255)

                                                              ggACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCctct                           >  1:5867364/1‑93 (MQ=255)

                                                               gACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCAtctc                          <  2:6941188/93‑1 (MQ=255)

                                                               gACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCAtctc                          >  2:2862572/1‑93 (MQ=255)

                                                               gACATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCAtctc                          <  1:6323682/93‑1 (MQ=255)

                                                                aCATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTcc                         <  2:605247/93‑1 (MQ=255)

                                                                aCATCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTcc                         <  2:6732879/93‑1 (MQ=255)

                                                                  aTCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCgg                       >  2:1618207/1‑93 (MQ=255)

                                                                   tCACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGt                      <  1:7712733/93‑1 (MQ=255)

                                                                    cACGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTg                     <  2:4349520/93‑1 (MQ=255)

                                                                     aCGCTGAAAAGTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTgg                    <  1:6177395/93‑1 (MQ=255)

                                                                      cGCTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGc                   <  1:3625795/93‑1 (MQ=255)

                                                                        cTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCaa                 >  2:857221/1‑93 (MQ=255)

                                                                        cTGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCaa                 >  1:2462126/1‑93 (MQ=255)

                                                                         tGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCaaa                >  1:6473018/1‑93 (MQ=255)

                                                                         tGAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGGGGCaaa                >  1:1204939/1‑93 (MQ=255)

                                                                          gAAAATTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCaaaa               <  2:1463407/93‑1 (MQ=255)

                                                                             aaTTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAAcc            >  2:518965/1‑93 (MQ=255)

                                                                             aaTTAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAAcc            >  1:4831320/1‑93 (MQ=255)

                                                                               ttAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATTTCCGGTGGCAAAAACCAc          >  2:7268793/1‑93 (MQ=255)

                                                                               ttAGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCAc          <  2:5484328/93‑1 (MQ=255)

                                                                                 aGGCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACgg        >  1:4337182/1‑93 (MQ=255)

                                                                                  ggCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGa       >  1:4452497/1‑93 (MQ=255)

                                                                                   gCGGAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGAc      >  1:1864770/1‑93 (MQ=255)

                                                                                      gAATGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGACGAt   <  2:6527509/93‑1 (MQ=255)

                                                                                       aaTGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGTCATAAACCACGGACGATa  <  2:7744121/93‑1 (MQ=255)

                                                                                       aaTGTTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGACGATa  <  2:4072147/93‑1 (MQ=255)

AAATGCCTGGCTCATGCGCGTGGCATCGACGTTCACCTGAATATCCGGCTGGCGTCCGGCGAGGACATCACGCTGAAAATTAGGCGGAATATTTATCGCGAAGGTATAGCGTCCGGCGTCCAGTCCGGCATCCATCTCATCGGCGGTGATCATCTCCGGTGGCAAAAACCACGGACGATA  >  NC_000913/3624045‑3624224

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