BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,941,734	Δ1 bp	noncoding (409/2906 nt)	rrlC →	23S ribosomal RNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,941,734	0	A	.	98.4%	312.8 / 2.9	182	noncoding (409/2906 nt)	rrlC	23S ribosomal RNA
Reads supporting (aligned to +/- strand): ref base A (0/0); major base . (109/70); minor base G (2/1); total (111/71)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.88e-01

Only 100 of 208 total aligned reads displayed.

GATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGAGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGA  >  NC_000913/3941644‑3941818

gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                                                                    >  2:357549/1‑93 (MQ=25)

gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                                                                    >  1:5783294/1‑93 (MQ=18)

gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                                                                    >  1:986148/1‑93 (MQ=25)

gATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggg                                                                                    >  2:2577134/1‑93 (MQ=25)

 aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATggggggc                                                                                   >  2:6574553/1‑92 (MQ=18)

 aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATggggggc                                                                                   >  2:7410574/1‑92 (MQ=18)

 aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggga                                                                                   <  2:5413273/93‑2 (MQ=17)

 aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggga                                                                                   >  1:1951241/1‑92 (MQ=21)

 aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggga                                                                                   >  2:8794439/1‑92 (MQ=18)

 aTACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGAGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATgggggga                                                                                   <  2:2772417/93‑2 (MQ=12)

  tACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAc                                                                                 >  2:8174517/1‑93 (MQ=12)

  tACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAc                                                                                 >  2:9753457/1‑93 (MQ=16)

  tACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGAc                                                                                 >  1:2593922/1‑93 (MQ=9)

    cAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCa                                                                               >  2:3804218/1‑93 (MQ=16)

    cAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCa                                                                               >  2:8714123/1‑93 (MQ=9)

    cAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCa                                                                               >  1:8198944/1‑93 (MQ=16)

    cAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCa                                                                               >  1:7475773/1‑93 (MQ=16)

     aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAt                                                                              >  2:9324222/1‑93 (MQ=16)

     aGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GAGACCAt                                                                              >  1:9379170/1‑93 (MQ=255)

      gggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtc                                                                             <  2:721485/93‑1 (MQ=9)

      gggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCCGCATATGGG‑GGGACCAtc                                                                             >  1:5240967/1‑93 (MQ=11)

       ggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcc                                                                            <  2:6831026/93‑1 (MQ=9)

       ggTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcc                                                                            <  1:6732391/93‑1 (MQ=9)

        gTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcct                                                                           >  2:1333088/1‑93 (MQ=16)

        gTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcct                                                                           <  1:10029079/93‑1 (MQ=9)

         tGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACTAtcctc                                                                          >  1:5016602/1‑93 (MQ=9)

         tGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcctc                                                                          <  1:9749369/93‑1 (MQ=9)

         tGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCAtcctc                                                                          <  1:8395293/93‑1 (MQ=9)

           aCAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCa                                                                        >  2:3313822/1‑93 (MQ=9)

            cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCaa                                                                       <  1:7939246/93‑1 (MQ=9)

            cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCaa                                                                       <  1:5728456/93‑1 (MQ=9)

            cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCaa                                                                       <  2:4126093/93‑1 (MQ=9)

            cAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCaa                                                                       >  1:2705428/1‑93 (MQ=16)

            cAGCCCCGTACACAAAAATGCACATATTGTGAGCCCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCaa                                                                       <  1:7236807/93‑1 (MQ=9)

            cAGCACCGTACAGAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCaa                                                                       <  1:7639696/93‑1 (MQ=11)

             aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAg                                                                      >  2:2147515/1‑93 (MQ=16)

             aGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAg                                                                      <  1:2675120/93‑1 (MQ=9)

              gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAgg                                                                     >  1:6033061/1‑93 (MQ=16)

              gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAgg                                                                     >  2:9523088/1‑93 (MQ=16)

              gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAgg                                                                     >  2:1866609/1‑93 (MQ=16)

              gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAgg                                                                     >  1:3682009/1‑93 (MQ=9)

              gCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAgc                                                                     >  1:2290366/1‑92 (MQ=9)

               ccccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGc                                                                    <  1:7554347/93‑1 (MQ=9)

               ccccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGc                                                                    <  1:4021924/93‑1 (MQ=9)

               ccccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGc                                                                    >  1:8710897/1‑93 (MQ=16)

                 ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTa                                                                  >  1:1137083/1‑93 (MQ=9)

                 ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTa                                                                  <  1:2765078/93‑1 (MQ=9)

                 ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTa                                                                  >  1:4261317/1‑93 (MQ=16)

                 ccGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTa                                                                  >  2:2822174/1‑93 (MQ=16)

                  cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaa                                                                 <  2:9185922/93‑1 (MQ=9)

                  cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaa                                                                 >  2:1070019/1‑93 (MQ=16)

                  cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaa                                                                 >  1:9912287/1‑93 (MQ=16)

                  cGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaa                                                                 <  1:2033323/93‑1 (MQ=9)

                   gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTTTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                                                                >  2:1455286/1‑93 (MQ=255)

                   gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                                                                >  2:9951917/1‑93 (MQ=16)

                   gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                                                                >  2:4453549/1‑93 (MQ=16)

                   gTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTaaa                                                                <  2:3887138/93‑1 (MQ=9)

                    tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAAt                                                               >  1:516606/1‑93 (MQ=9)

                    tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAAt                                                               >  2:6356222/1‑93 (MQ=12)

                    tACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAAt                                                               >  1:9159327/1‑93 (MQ=9)

                     acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                                                              >  2:8350931/1‑93 (MQ=16)

                     acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                                                              >  2:6502692/1‑93 (MQ=16)

                     acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                                                              <  1:4905150/93‑1 (MQ=9)

                     acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                                                              <  2:679064/93‑1 (MQ=9)

                     acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                                                              <  2:7776320/93‑1 (MQ=9)

                     acacAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATa                                                              <  1:7987062/93‑1 (MQ=9)

                      cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATAc                                                             >  1:6127570/1‑93 (MQ=16)

                      cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATAc                                                             >  2:3837796/1‑93 (MQ=16)

                      cacaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATAc                                                             <  2:2613161/93‑1 (MQ=9)

                       acaAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACt                                                            >  2:3551525/1‑93 (MQ=16)

                        caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTc                                                           <  1:6613601/93‑1 (MQ=9)

                        caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTc                                                           <  1:3858374/93‑1 (MQ=9)

                        caAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTc                                                           >  2:9743928/1‑93 (MQ=16)

                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTcc                                                          <  1:3956430/93‑1 (MQ=9)

                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTcc                                                          <  2:5894525/93‑1 (MQ=9)

                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTcc                                                          >  1:282866/1‑93 (MQ=16)

                         aaaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTcc                                                          >  2:4543727/1‑93 (MQ=16)

                          aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCct                                                         >  2:4786991/1‑93 (MQ=16)

                          aaaaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCct                                                         >  2:7062376/1‑93 (MQ=9)

                            aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGATTATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctga                                                       <  1:1166867/93‑1 (MQ=9)

                            aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctga                                                       >  1:4023853/1‑93 (MQ=16)

                            aaTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctga                                                       <  1:9812650/93‑1 (MQ=9)

                             aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgac                                                      <  1:4252903/93‑1 (MQ=9)

                             aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgac                                                      <  1:2221429/93‑1 (MQ=9)

                             aTGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgac                                                      <  2:2762317/93‑1 (MQ=9)

                              tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgact                                                     >  1:6107587/1‑93 (MQ=9)

                              tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgact                                                     >  1:6461002/1‑93 (MQ=12)

                              tGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgact                                                     >  1:3280698/1‑93 (MQ=16)

                               gCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgactg                                                    >  2:8333734/1‑93 (MQ=16)

                                cacaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCctgactga                                                   <  2:221183/93‑1 (MQ=9)

                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTcctgcctgcc                                                  >  1:9022995/1‑91 (MQ=11)

                                 acaTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCtgactgac                                                  <  2:4302176/93‑1 (MQ=9)

                                  caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc                                                 >  1:1630832/1‑93 (MQ=0)

                                  caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc                                                 >  2:1773589/1‑93 (MQ=0)

                                  caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc                                                 >  2:7377944/1‑93 (MQ=0)

                                  caTATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGG‑GGGACCATCCTCCAAGGCTAAATACTCCTGACTGAcc                                                 >  1:7124957/1‑93 (MQ=0)

                                                              gACACGTGGTCTCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATGCTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGaaa                      <  2:1968382/93‑1 (MQ=255)

                                                                  cGTGGTATCCTGTCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCg                  >  2:582387/1‑93 (MQ=17)

                                                                              tCTGAATATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCg      >  2:516165/1‑93 (MQ=17)

                                                                                  aaTATGGGGGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGa  <  2:6135675/93‑1 (MQ=18)

GATACAGGGTGACAGCCCCGTACACAAAAATGCACATATTGTGAGCTCGATGAGTAGGGCGGGACACGTGGTATCCTGTCTGAATATGGGAGGGACCATCCTCCAAGGCTAAATACTCCTGACTGACCGATAGTGAACCAGTACCGTGAGGGAAAGGCGAAAAGAACCCCGGCGA  >  NC_000913/3941644‑3941818

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