BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,807,191	G→A	P23S (CCT→TCT)	b1726 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,807,191	0	G	A	90.0%	639.3 / 41.9	201	P23S (CCT→TCT)	b1726	hypothetical protein
Reads supporting (aligned to +/- strand): ref base G (13/7); new base A (97/84); total (110/91)
Fisher's exact test for biased strand distribution p-value = 3.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 203 total aligned reads displayed.

AGTCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGACGCGTCG  >  NC_000913/1807101‑1807283

aGTCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGAta                                                                                            >  2:164057/1‑93 (MQ=255)

aGTCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGAta                                                                                            >  2:1673780/1‑93 (MQ=255)

aGTCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGAta                                                                                            <  2:3340980/93‑1 (MQ=255)

 gTCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGtat                                                                                           >  2:4065962/1‑93 (MQ=255)

  tCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGAtata                                                                                          <  1:730631/93‑1 (MQ=255)

   cGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATaa                                                                                         <  2:6333615/93‑1 (MQ=255)

   cGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATaa                                                                                         <  2:5860831/93‑1 (MQ=255)

   cGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATaa                                                                                         >  1:49793/1‑93 (MQ=255)

    gAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATaaa                                                                                        <  2:4244126/93‑1 (MQ=255)

    gAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATaaa                                                                                        >  1:2922825/1‑93 (MQ=255)

     aGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATaaaa                                                                                       >  2:2198185/1‑93 (MQ=255)

     aGCATGACCGCAATAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATaaaa                                                                                       <  2:4743737/93‑1 (MQ=255)

      gCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATaaaaa                                                                                      >  1:6518913/1‑93 (MQ=255)

       cATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAAt                                                                                     <  2:4103065/93‑1 (MQ=255)

         tGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCa                                                                                   >  1:2251476/1‑93 (MQ=255)

           aCCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCacc                                                                                 >  2:1793653/1‑93 (MQ=255)

            ccGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCaccc                                                                                <  2:3115578/93‑1 (MQ=255)

              gCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCacccaa                                                                              >  1:5839373/1‑93 (MQ=255)

              gCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCacccaa                                                                              <  1:3111156/93‑1 (MQ=255)

               cATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCacccaac                                                                             >  1:2501861/1‑93 (MQ=255)

                 ttAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCacccaaccc                                                                           >  1:1480006/1‑93 (MQ=255)

                     tGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCaaaa                                                                       >  1:2017505/1‑93 (MQ=255)

                        cttccttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATcc                                                                    <  1:1427802/93‑1 (MQ=255)

                        cttccttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATct                                                                    >  2:1877622/1‑92 (MQ=255)

                        cttccttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATcc                                                                    <  1:542609/93‑1 (MQ=255)

                         ttccttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCt                                                                   >  1:1160626/1‑93 (MQ=255)

                         ttccttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCt                                                                   >  2:268120/1‑93 (MQ=255)

                          tccttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCtt                                                                  >  1:3129303/1‑93 (MQ=255)

                          tccttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCtt                                                                  <  2:4324042/93‑1 (MQ=255)

                          tccttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCtt                                                                  >  2:2877056/1‑93 (MQ=255)

                            cttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCtttt                                                                <  2:3875999/93‑1 (MQ=255)

                             ttGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTc                                                               <  2:2072987/93‑1 (MQ=255)

                                 tttCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCaaaa                                                           <  2:4224908/93‑1 (MQ=255)

                                 tttCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCaaaa                                                           <  1:5465005/93‑1 (MQ=255)

                                  ttCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAAc                                                          >  2:5398640/1‑93 (MQ=255)

                                  ttCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAAc                                                          <  2:5436880/93‑1 (MQ=255)

                                     cTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAACGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGc                                                       >  2:2500843/1‑93 (MQ=255)

                                       gCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCaa                                                     >  1:5551192/1‑93 (MQ=255)

                                        ccGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAAt                                                    >  1:5768010/1‑93 (MQ=255)

                                        ccGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAAt                                                    >  1:2602715/1‑93 (MQ=255)

                                         cGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATg                                                   <  1:4877177/93‑1 (MQ=255)

                                         cGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATg                                                   <  2:5320974/93‑1 (MQ=255)

                                         cGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATg                                                   <  1:2663463/93‑1 (MQ=255)

                                          gggTATTCATACATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGc                                                  >  1:2115287/1‑93 (MQ=255)

                                          gggTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGc                                                  <  2:708007/93‑1 (MQ=255)

                                          gggTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGc                                                  >  1:5496771/1‑93 (MQ=255)

                                              aTTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTcc                                              <  1:4077465/93‑1 (MQ=255)

                                              aTTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTcc                                              <  1:9206/93‑1 (MQ=255)

                                               ttCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTGCa                                             >  1:3845627/1‑93 (MQ=255)

                                                  aTAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCt                                          <  2:2468969/93‑1 (MQ=255)

                                                  aTAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCt                                          >  2:5225621/1‑93 (MQ=255)

                                                    aaaTTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCttg                                        >  2:2104503/1‑93 (MQ=255)

                                                      aTTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCttgtt                                      >  1:4259269/1‑93 (MQ=255)

                                                      aTTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCttgtt                                      <  2:3763764/93‑1 (MQ=255)

                                                      aTTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCttgtt                                      >  1:5251072/1‑93 (MQ=255)

                                                      aTTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCttgtt                                      >  1:5372869/1‑93 (MQ=255)

                                                      aTTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCttgtt                                      >  2:4476937/1‑93 (MQ=255)

                                                       tttCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCttgttg                                     >  2:1319609/1‑93 (MQ=255)

                                                         tCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGAt                                   >  2:4900791/1‑93 (MQ=255)

                                                         tCAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGAt                                   <  1:3236236/93‑1 (MQ=255)

                                                          cAGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATa                                  >  1:1897055/1‑93 (MQ=255)

                                                          cAGCAAAGAAATCAATGTAGAGATAAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATa                                  >  1:2474102/1‑93 (MQ=255)

                                                           aGCAAAGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATAc                                 <  2:4577132/93‑1 (MQ=255)

                                                            gCAAAGAAGTCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACg                                >  1:1894019/1‑93 (MQ=255)

                                                              aaaGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTc                              >  1:4744578/1‑93 (MQ=255)

                                                               aaGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAGTCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCa                             <  2:4666285/93‑1 (MQ=255)

                                                               aaGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCa                             >  1:4120840/1‑93 (MQ=255)

                                                               aaGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCa                             <  1:4026708/93‑1 (MQ=255)

                                                               aaGAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCa                             <  1:5093756/93‑1 (MQ=255)

                                                                aGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCAt                            <  2:4088311/93‑1 (MQ=255)

                                                                 gAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATg                           >  1:5418381/1‑93 (MQ=255)

                                                                 gAAATCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATg                           >  2:1213233/1‑93 (MQ=255)

                                                                  aaaTCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGa                          <  2:3514787/93‑1 (MQ=255)

                                                                  aaaTCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGa                          >  1:5698544/1‑93 (MQ=255)

                                                                    aTCAATGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAACTCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGt                        <  2:3421650/93‑1 (MQ=255)

                                                                       aaTGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTtct                     >  1:5889404/1‑93 (MQ=255)

                                                                       aaTGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTtct                     >  1:1833189/1‑93 (MQ=255)

                                                                       aaTGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTtct                     >  1:1660304/1‑93 (MQ=255)

                                                                         tGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTcc                   >  1:3668021/1‑93 (MQ=255)

                                                                         tGTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTcc                   <  2:5761855/93‑1 (MQ=255)

                                                                          gTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCt                  <  1:100436/93‑1 (MQ=255)

                                                                          gTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCt                  <  1:5272157/93‑1 (MQ=255)

                                                                          gTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCt                  <  2:6083448/93‑1 (MQ=255)

                                                                          gTAGAGATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCt                  <  1:3332637/93‑1 (MQ=255)

                                                                            agagATCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTg                >  2:6056845/1‑93 (MQ=255)

                                                                             gagaTCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGa               >  1:1116922/1‑93 (MQ=255)

                                                                              agaTCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGAt              >  2:6282135/1‑93 (MQ=255)

                                                                               gaTCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGAtt             >  2:3136088/1‑93 (MQ=255)

                                                                                 tCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCttgttg                                     <  1:124567/67‑1 (MQ=255)

                                                                                 tCAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTaa           >  1:642441/1‑93 (MQ=255)

                                                                                  cAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAg          >  1:5717145/1‑93 (MQ=255)

                                                                                  cAATGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAg          <  2:1590408/93‑1 (MQ=255)

                                                                                   aaTGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAATTTGCAATTCGTTCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGa         >  2:6640544/1‑93 (MQ=255)

                                                                                   aaTGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGa         <  1:356013/93‑1 (MQ=255)

                                                                                    aTGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGAc        <  1:5995118/93‑1 (MQ=255)

                                                                                    aTGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGAc        <  2:6133244/93‑1 (MQ=255)

                                                                                    aTGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGAc        >  1:4354396/1‑93 (MQ=255)

                                                                                    aTGCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGAc        <  2:642401/93‑1 (MQ=255)

                                                                                      gCAGATATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGAcgc      <  2:5291083/93‑1 (MQ=255)

                                                                                          gTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGACGcgtcg  <  2:3973358/93‑1 (MQ=255)

AGTCGAGCATGACCGCATTAATGCCTTCCTTGTTTTCCTGCCGGGTATTCATAAATTTCAGCAAAGAAATCAATGTAGAGATCAATGCAGGTATAAAAATCACCCAACCCAAAATCCTTTTCAAAACAGCAATGCGTCCAGCTTGTTGATACGTCATGAGTTCTCCTTGATTAAGACGCGTCG  >  NC_000913/1807101‑1807283

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