BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,566,541	T→C	I103V (ATA→GTA)	b1491 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,566,541	0	T	C	99.1%	369.6 / ‑4.7	116	I103V (ATA→GTA)	b1491	hypothetical protein
Reads supporting (aligned to +/- strand): ref base T (1/0); new base C (49/65); total (51/65)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 121 total aligned reads displayed.

TTACCGGTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAG  >  NC_000913/1566453‑1566629

ttACCGGTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCg                                                                                      >  1:3381542/1‑93 (MQ=255)

ttACCGGTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCg                                                                                      >  2:232546/1‑93 (MQ=255)

ttACCGGACATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCg                                                                                      >  2:1556374/1‑93 (MQ=255)

 tACCGGTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCga                                                                                     <  1:799402/93‑1 (MQ=255)

   ccGGTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGTAAAAGACCGTGTTTACGCCgaga                                                                                   >  2:2588563/1‑93 (MQ=255)

   ccGGTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCgaga                                                                                   <  2:1363119/93‑1 (MQ=255)

     ggTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACg                                                                                 >  1:3695172/1‑93 (MQ=255)

     ggTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACg                                                                                 >  1:4184244/1‑93 (MQ=255)

     ggTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACg                                                                                 <  2:341607/93‑1 (MQ=255)

      gTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGt                                                                                <  1:3004209/93‑1 (MQ=255)

      gTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGt                                                                                <  2:1116456/93‑1 (MQ=255)

       tCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGtt                                                                               <  1:1626334/93‑1 (MQ=255)

         aTAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGc                                                                             <  2:833240/93‑1 (MQ=255)

         aTAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGc                                                                             <  1:1128259/93‑1 (MQ=255)

         aTAAGATCGGACCACAGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGc                                                                             <  2:4212795/93‑1 (MQ=255)

          tAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCa                                                                            >  2:2675499/1‑93 (MQ=255)

          tAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCa                                                                            >  2:2115239/1‑93 (MQ=255)

          tAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCa                                                                            >  1:3438390/1‑93 (MQ=255)

          tAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCa                                                                            <  1:214559/93‑1 (MQ=255)

            aGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAga                                                                          <  2:2891419/93‑1 (MQ=255)

            aGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAga                                                                          <  1:1329804/93‑1 (MQ=255)

             gATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAgat                                                                         <  1:2461781/93‑1 (MQ=255)

               tCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAgatga                                                                       >  1:1652963/1‑93 (MQ=255)

                 ggACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATc                                                                     >  1:2178727/1‑93 (MQ=255)

                  gACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATcc                                                                    >  1:159610/1‑93 (MQ=255)

                   aCCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGGTGATCCa                                                                   >  1:1710406/1‑93 (MQ=255)

                   aCCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCa                                                                   <  2:3458407/93‑1 (MQ=255)

                     cACGGCAAATTTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGAGTTTAAGCCGAGACGGTGCAGATGTTTCAGt                                                                 >  2:2007376/1‑93 (MQ=255)

                     cACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGt                                                                 >  1:956429/1‑93 (MQ=255)

                         gCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGt                                                             >  2:740892/1‑93 (MQ=255)

                          cAAAATTTTCGATGGCCCCAGGGCGGTTCCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTc                                                            <  2:3823387/93‑1 (MQ=255)

                          cAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGTCGAGACGTTGCAGATGATCCAGTTTGTc                                                            >  1:761737/1‑93 (MQ=255)

                             aaTTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGTAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGat                                                         <  2:265888/93‑1 (MQ=255)

                             aaTTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGat                                                         <  2:2081135/93‑1 (MQ=255)

                             aaTTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGat                                                         <  2:265889/93‑1 (MQ=255)

                             aaTTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGat                                                         <  1:672621/93‑1 (MQ=255)

                               ttttCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGatca                                                       <  1:4363269/93‑1 (MQ=255)

                               ttttCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGatca                                                       >  2:2397540/1‑93 (MQ=255)

                               ttttCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGatca                                                       >  1:1478398/1‑93 (MQ=255)

                               ttttCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGatca                                                       <  2:3712685/93‑1 (MQ=255)

                               ttttCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGatca                                                       <  2:3712686/93‑1 (MQ=255)

                                tttCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGatcat                                                      <  2:3943821/93‑1 (MQ=255)

                                 ttCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGatcatc                                                     >  1:3953657/1‑93 (MQ=255)

                                    gATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATAATCGcc                                                  <  2:4486483/93‑1 (MQ=255)

                                    gATGGCCACAGGGCGGTACCGTCCGGCGTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGcc                                                  <  1:1417631/93‑1 (MQ=255)

                                     aTGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCt                                                 >  1:2238505/1‑93 (MQ=255)

                                     aTGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCt                                                 <  2:4599143/93‑1 (MQ=255)

                                      tGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCtg                                                <  1:3865774/93‑1 (MQ=255)

                                       ggCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCtgt                                               >  1:806947/1‑93 (MQ=255)

                                       ggCAACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCtgt                                               <  1:2856164/93‑1 (MQ=255)

                                        gCCACAGGGCGTTACCGTCCGGCATGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCtgtt                                              <  2:3200996/93‑1 (MQ=255)

                                          cacaGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCtgttgt                                            <  1:1972191/93‑1 (MQ=255)

                                          cacaGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCtgttgt                                            >  2:1455138/1‑93 (MQ=255)

                                             aGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCtgttgttgt                                         <  2:2333998/93‑1 (MQ=255)

                                             aGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCtgttgttgt                                         <  1:4537167/93‑1 (MQ=255)

                                             aGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCtgttgttgt                                         <  2:2078296/93‑1 (MQ=255)

                                              gggCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTa                                        <  2:1346494/93‑1 (MQ=255)

                                                gCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTaca                                      <  2:4073088/93‑1 (MQ=255)

                                                  ggTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACg                                    <  1:1701470/93‑1 (MQ=255)

                                                  ggTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACg                                    >  2:1012530/1‑93 (MQ=255)

                                                  ggTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACg                                    <  1:1061889/93‑1 (MQ=255)

                                                   gTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACgg                                   <  2:2741006/93‑1 (MQ=255)

                                                   gTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACgg                                   >  1:3492242/1‑93 (MQ=255)

                                                   gTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACgg                                   <  2:4586537/93‑1 (MQ=255)

                                                   gTACCGTCCGGCTTGACCTGGAAAAAGACCGGGGTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACgg                                   >  2:694582/1‑93 (MQ=255)

                                                    tACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACggg                                  <  2:1105447/93‑1 (MQ=255)

                                                     aCCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGc                                 >  1:102308/1‑93 (MQ=255)

                                                      ccGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGcc                                <  1:3252642/93‑1 (MQ=255)

                                                      ccGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGcc                                >  2:4423506/1‑93 (MQ=255)

                                                      ccGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGcc                                <  2:39889/93‑1 (MQ=255)

                                                       cGTCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGccc                               <  1:2452433/93‑1 (MQ=255)

                                                         tCCGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCgg                             >  1:4004262/1‑93 (MQ=255)

                                                          ccGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGc                            <  1:3180731/93‑1 (MQ=255)

                                                          ccGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGc                            <  2:102077/93‑1 (MQ=255)

                                                          ccGGCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGc                            <  1:1017820/93‑1 (MQ=255)

                                                          ccGGCTTGACCTGGAAAAACACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTTTACACGGGCCCGGc                            <  1:3180730/93‑1 (MQ=255)

                                                             gCTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtgg                         <  1:3829474/93‑1 (MQ=255)

                                                              cTTGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtggt                        <  2:3432064/93‑1 (MQ=255)

                                                                tGACCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCtggtgg                      >  2:744677/1‑93 (MQ=255)

                                                                 gACCTGGAAAAAGACCGTGTTGACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTggg                     <  2:2302586/93‑1 (MQ=255)

                                                                  aCCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTgggg                    <  2:4505919/93‑1 (MQ=255)

                                                                  aCCTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGCCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTgggg                    <  2:1920607/93‑1 (MQ=255)

                                                                   ccTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGt                   <  1:263534/93‑1 (MQ=255)

                                                                    cTGGAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGtt                  >  1:2151978/1‑93 (MQ=255)

                                                                      ggAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTAc                >  1:1625677/1‑93 (MQ=255)

                                                                       gAAAAAGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACt               >  1:104210/1‑93 (MQ=255)

                                                                         aaaaGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGGGGGGTTACTaa             >  2:3767254/1‑93 (MQ=255)

                                                                         aaaaGACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCACGGCTGGTGGGGTTACTat             >  2:3589930/1‑92 (MQ=255)

                                                                             gACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGtt         >  1:2543795/1‑93 (MQ=255)

                                                                             gACCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGtt         >  1:3747651/1‑93 (MQ=255)

                                                                              aCCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTa        >  1:1630170/1‑93 (MQ=255)

                                                                              aCCGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGCTGTTGTAGACGTCACCTCTCGGTGGGGTTACTAATGTTa        >  1:1630169/1‑93 (MQ=255)

                                                                               ccGTGTTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTaa       >  1:1194886/1‑93 (MQ=255)

                                                                                 gtgtTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAAcc     <  2:4229415/93‑1 (MQ=255)

                                                                                 gtgtTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAAcc     >  1:446437/1‑93 (MQ=255)

                                                                                   gtTTACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGa   >  1:1450580/1‑93 (MQ=255)

                                                                                    tttACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAg  >  1:4643161/1‑93 (MQ=255)

                                                                                    tttACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAg  <  1:286004/93‑1 (MQ=255)

                                                                                    tttACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAg  >  1:4099326/1‑93 (MQ=255)

                                                                                    tttACGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAg  >  1:2014081/1‑93 (MQ=255)

TTACCGGTCATAAGATCGGACCACGGCAAAATTTTCGATGGCCACAGGGCGGTACCGTCCGGCTTGACCTGGAAAAAGACCGTGTTTATGCCGAGACGTTGCAGATGATCCAGTTTGTCGATCATCGCCTGTTGTTGTACACGGGCCCGGCTGGTGGGGTTACTAATGTTAACCGAG  >  NC_000913/1566453‑1566629

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