BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,071,491	T→A	9.2%	noncoding (175/195 nt)	isrC →	small RNA IsrC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,071,491	0	T	A	9.2%	284.7 / 15.7	100	noncoding (175/195 nt)	isrC	small RNA IsrC
Reads supporting (aligned to +/- strand): ref base T (42/48); new base A (4/5); total (46/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.10e-01

Only 100 of 105 total aligned reads displayed.

TTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGC  >  NC_000913/2071342‑2071640

ttATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGt                                                                                                                                                       <  2:163583/150‑1 (MQ=255)

ttATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGt                                                                                                                                                       <  2:157806/150‑1 (MQ=255)

 tataTCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGtt                                                                                                                                                      <  1:202530/150‑1 (MQ=255)

 tataTCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGtt                                                                                                                                                      >  1:1108081/1‑150 (MQ=255)

      cgatcgatAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTgg                                                                                                                                                 <  1:37806/150‑1 (MQ=255)

      cgatcgatAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTgg                                                                                                                                                 <  1:757033/150‑1 (MQ=255)

      cgatcgatAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTgg                                                                                                                                                 >  2:1346070/1‑150 (MQ=255)

           gatAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTAcc                                                                                                                                            <  1:125368/150‑1 (MQ=255)

             tAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTACCCACAAACGTGGATACCgg                                                                                                                                          <  2:398940/150‑1 (MQ=255)

             tAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTACCCACAAACGTGGATACCgg                                                                                                                                          <  1:1236015/150‑1 (MQ=255)

               aGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTggg                                                                                                                                                >  2:1390969/1‑142 (MQ=255)

               aGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTggg                                                                                                                                                <  1:1390969/142‑1 (MQ=255)

               aGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTACCCACAAACGTGGATACCGGca                                                                                                                                        <  1:585683/150‑2 (MQ=255)

                gCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTACCCACAAACGTGGATACCGGcag                                                                                                                                       <  1:1297104/150‑3 (MQ=255)

                gCTAATAATAAACTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTAAGGGAAAAACAGAAACACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCtt                                                                                                                                       >  2:207054/1‑150 (MQ=255)

                   aataataaCCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCttttt                                                                                                                                    <  1:1073210/150‑1 (MQ=255)

                     taataaCCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTa                                                                                                                                  <  1:1580735/150‑1 (MQ=255)

                     taataaCCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTa                                                                                                                                  <  1:483738/150‑1 (MQ=255)

                       ataaCCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTAtt                                                                                                                                <  1:228490/150‑1 (MQ=255)

                       ataaCCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTAtt                                                                                                                                <  1:529451/150‑1 (MQ=255)

                        taaCCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTc                                                                                                                               <  2:362623/150‑1 (MQ=255)

                          aCCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCAc                                                                                                                             <  2:522150/150‑1 (MQ=255)

                           ccTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCAcc                                                                                                                            <  2:531171/150‑1 (MQ=255)

                             tttGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCt                                                                                                                          <  2:240044/150‑1 (MQ=255)

                              ttGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTc                                                                                                                         <  1:461789/150‑1 (MQ=255)

                                  cAGTAACATGCACAGATACGTACAGAAAGACATTAAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATc                                                                                                                     >  2:1494087/1‑150 (MQ=255)

                                   aGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCt                                                                                                                    >  1:1496999/1‑150 (MQ=255)

                                    gTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTa                                                                                                                   <  2:923071/150‑1 (MQ=255)

                                      aaCATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACACTCAATCTAAg                                                                                                                 >  2:318293/1‑150 (MQ=255)

                                          tGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGaaa                                                                                                             >  2:521185/1‑150 (MQ=255)

                                              caGATACGTACAGAAAGACATTCAGCGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCt                                                                                                         >  2:1101263/1‑150 (MQ=255)

                                               aGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCtg                                                                                                        <  2:1008244/150‑1 (MQ=255)

                                                 aTACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCtgat                                                                                                      <  1:786010/150‑1 (MQ=255)

                                                 aTACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCtgat                                                                                                      <  2:1527984/150‑1 (MQ=255)

                                                  tACGTACAGAAAGACGTTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCtgatg                                                                                                     <  2:1360219/150‑1 (MQ=255)

                                                  tACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCtgatg                                                                                                     <  2:865466/150‑1 (MQ=255)

                                                    cGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGCACCTCCGGCACAGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGaa                                                                                                   >  1:838474/1‑150 (MQ=255)

                                                       aCAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACg                                                                                                >  1:427065/1‑150 (MQ=255)

                                                        cAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGa                                                                                               >  1:701607/1‑150 (MQ=255)

                                                          gAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACa                                                                                             >  2:1285816/1‑150 (MQ=255)

                                                              gACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTg                                                                                         <  2:1523923/150‑1 (MQ=255)

                                                                cATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGaa                                                                                       >  1:892509/1‑150 (MQ=255)

                                                                 aTTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAAt                                                                                      >  1:321361/1‑150 (MQ=255)

                                                                 aTTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAAt                                                                                      >  1:1040917/1‑150 (MQ=255)

                                                                     aGGGAACAACAGAACGACAATTCATAAACTCCCACAGCCGGACCTCCGGCACTGTAACGCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGATACGACATCTGAATAcca                                                                                  >  2:1326904/1‑149 (MQ=255)

                                                                     aGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCt                                                                                  <  1:895331/150‑1 (MQ=255)

                                                                     aGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCt                                                                                  <  2:410309/150‑1 (MQ=255)

                                                                        gAACAACATAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCt                                                                               <  1:760013/150‑1 (MQ=255)

                                                                            aacaGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAg                                                                           >  2:1391833/1‑150 (MQ=255)

                                                                            aacaGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAg                                                                           >  2:243859/1‑150 (MQ=255)

                                                                              caGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGc                                                                         >  2:1418396/1‑150 (MQ=255)

                                                                               aGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCt                                                                        >  1:172843/1‑150 (MQ=255)

                                                                                 aaCCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTgg                                                                      <  1:216105/150‑1 (MQ=255)

                                                                                 aaCCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTgg                                                                      <  2:1533261/150‑1 (MQ=255)

                                                                                 aaCCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTgg                                                                      >  1:463964/1‑150 (MQ=255)

                                                                                  aCCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGt                                                                     <  2:126673/150‑1 (MQ=255)

                                                                                    cacaATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTAt                                                                   >  2:36654/1‑150 (MQ=255)

                                                                                     acaATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATg                                                                  >  2:1038082/1‑150 (MQ=255)

                                                                                     acaATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATg                                                                  <  1:1101263/150‑1 (MQ=255)

                                                                                      caATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATgg                                                                 >  2:793619/1‑150 (MQ=255)

                                                                                           cAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATca                                                            >  2:453219/1‑150 (MQ=255)

                                                                                           cAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATca                                                            <  1:901577/150‑1 (MQ=255)

                                                                                           cAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATca                                                            <  1:1346070/150‑1 (MQ=255)

                                                                                             gAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATcaca                                                          <  1:1418396/150‑1 (MQ=255)

                                                                                             gAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATcaca                                                          >  1:380007/1‑150 (MQ=255)

                                                                                               aaCTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTTTGGAATCACATg                                                        >  1:1162739/1‑150 (MQ=255)

                                                                                                 cTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGAc                                                      >  1:705759/1‑150 (MQ=255)

                                                                                                  tCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACg                                                     >  1:488706/1‑150 (MQ=255)

                                                                                                     cagtgCCGGAGGTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGgc                                                  >  2:1535622/5‑150 (MQ=255)

                                                                                                       caGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGgcgc                                                >  2:341861/1‑150 (MQ=255)

                                                                                                         gCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCtt                                              >  2:109285/1‑150 (MQ=255)

                                                                                                             gACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGt                                          <  1:1495933/150‑1 (MQ=255)

                                                                                                               ccTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTgg                                        >  1:347069/1‑150 (MQ=255)

                                                                                                               ccTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTgg                                        <  1:168159/150‑1 (MQ=255)

                                                                                                                cTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGt                                       <  1:1612941/150‑1 (MQ=255)

                                                                                                                 tCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGtt                                      <  1:207054/150‑1 (MQ=255)

                                                                                                                  ccGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTg                                     >  1:801275/1‑150 (MQ=255)

                                                                                                                   cGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGc                                    >  1:1502991/1‑150 (MQ=255)

                                                                                                                       aCTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTcc                                <  2:1108081/150‑1 (MQ=255)

                                                                                                                           tAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAAc                            <  1:1317809/150‑1 (MQ=255)

                                                                                                                               ccTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGc                        >  1:1225985/1‑150 (MQ=255)

                                                                                                                                  ttACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCcg                     <  2:347069/150‑1 (MQ=255)

                                                                                                                                   tACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCcgc                    >  1:973175/1‑150 (MQ=255)

                                                                                                                                     aaagCCGGTACCCACAAACGTGGATACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCgcgc                  >  2:1231363/4‑150 (MQ=255)

                                                                                                                                      aagCCGGTACCCACAAACGTGGATACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCgcgcc                 >  1:1353220/3‑149 (MQ=255)

                                                                                                                                       agCCGGTACCCACAAACGTGGATACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAACCGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGGTGCCTCCGAACTGGCCCgcgcgc                >  1:719607/2‑148 (MQ=255)

                                                                                                                                        gCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACg               >  1:1022936/1‑150 (MQ=255)

                                                                                                                                        gCCGGTACCCACAAACGTGGATACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTg                                     <  1:586864/128‑1 (MQ=255)

                                                                                                                                        gCCGGTACCCACAAACGTGGATACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTg                                     >  2:586864/1‑128 (MQ=255)

                                                                                                                                         ccGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACgg              >  1:1458921/1‑150 (MQ=255)

                                                                                                                                          cGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACggg             <  2:890488/150‑1 (MQ=255)

                                                                                                                                          cGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACggg             <  1:109285/150‑1 (MQ=255)

                                                                                                                                           ggTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACgggg            <  1:867788/150‑1 (MQ=255)

                                                                                                                                           ggTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACgggg            <  1:1285816/150‑1 (MQ=255)

                                                                                                                                            gTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGt           >  1:1406649/1‑150 (MQ=255)

                                                                                                                                            gTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGt           <  1:1038082/150‑1 (MQ=255)

                                                                                                                                            gTATCCACG‑TTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTa          <  1:639188/150‑1 (MQ=255)

                                                                                                                                               tCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGTTTGCCTCCGAACTGGCCCGCGCACGGGGTaaa        <  2:939288/150‑1 (MQ=255)

                                                                                                                                                   cGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACgtg    <  1:793619/150‑1 (MQ=255)

                                                                                                                                                     tttGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACgtggc  <  2:120894/150‑1 (MQ=255)

TTATATCGATCGATAAGCTAATAATAACCTTTGTCAGTAACATGCACAGATACGTACAGAAAGACATTCAGGGAACAACAGAACCACAATTCAGAAACTCCCACAGCCGGACCTCCGGCACTGTAACCCTTTACCTGCCGGTATCCACGTTTGTGGGTACCGGCTTTTTTATTCACCCTCAATCTAAGGAAAAGCTGATGAAACGACATCTGAATACCTGCTACAGGCTGGTATGGAATCACATGACGGGCGCTTTCGTGGTTGCCTCCGAACTGGCCCGCGCACGGGGTAAACGTGGC  >  NC_000913/2071342‑2071640

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