Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,966,720	C→T	100%	R102C (CGC→TGC)	rho →	transcription termination factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,966,720	0	C	T	100.0%	275.3 / NA	83	R102C (CGC→TGC)	rho	transcription termination factor
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (54/29);  total (54/29)

CAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTT  >  NC_000913/3966623‑3966815                                                                                                  |                                                                                                cAGCTCCTACCTCGCCTGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgc                                                                                               <  2:1550240/100‑1 (MQ=255)   gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgccc                                                                                             >  1:2375945/1‑100 (MQ=255)   gCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgcc                                                                                              <  2:554759/99‑1 (MQ=255)    cTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgcccg                                                                                            <  1:1527102/100‑1 (MQ=255)      cgatcttCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgcccgcc                                                                                          <  2:1259882/94‑1 (MQ=255)      ccTACCTCGCCGGTCCTGTTGACACCTACGTTTCCCTTAGCCAAATCCGCCGTTTAATGCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgcccgcc                                                                                          <  2:578549/100‑1 (MQ=255)      ccTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgcccgcc                                                                                          >  2:2229220/1‑100 (MQ=255)        tACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGa                                                                                        >  1:1919825/1‑100 (MQ=255)           cTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaag                                                                                     <  2:1835986/100‑1 (MQ=255)           cTCGCCGGTCCTGATGACATCTACGTGTCACCTAGCCCAATCCCCAGTTTCAACCCCAGCACTGGGGATACTCTTCCTGGTCAGCTTTGCCGGCCGaac                                                                                      >  2:2089975/1‑98 (MQ=255)            tCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTtgcccgccagatcg                                                                                    >  1:1259889/1‑94 (MQ=255)             cGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAaagaa                                                                                   >  2:1094333/1‑100 (MQ=255)               ccGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAgg                                                                                 >  1:1078734/1‑100 (MQ=255)               ccGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAgg                                                                                 >  1:1479633/1‑100 (MQ=255)                    cttgatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAAc                                                                            <  1:778384/98‑1 (MQ=255)                     cTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACg                                                                           >  2:2387196/1‑100 (MQ=255)                      tgatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAAcag                                                                          >  2:778382/1‑98 (MQ=255)                      tgatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGc                                                                          >  2:1462580/1‑100 (MQ=255)                       gatgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGACAGAAGGTGAACGCt                                                                         >  2:1088830/1‑100 (MQ=255)                         tgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTAt                                                                       >  1:1590094/1‑100 (MQ=255)                         tgaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTAt                                                                       >  1:1243311/1‑100 (MQ=255)                          gaCATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTAtt                                                                      <  1:1155883/100‑1 (MQ=255)                            cATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTAtttt                                                                    >  1:2017258/1‑100 (MQ=255)                             aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTAtttt                                                                    >  2:954042/1‑99 (MQ=255)                             aTCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTg                                                                   <  1:2187847/100‑1 (MQ=255)                              tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                  <  2:1243304/100‑1 (MQ=255)                              tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                  >  1:2200485/1‑100 (MQ=255)                              tCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgc                                                                  >  1:242618/1‑100 (MQ=255)                               cTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgcg                                                                 >  1:639688/1‑100 (MQ=255)                                tACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTgcgc                                                                <  1:1688701/100‑1 (MQ=255)                                  cGTTTCCCCTCGCCACATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgct                                                               <  2:883483/99‑1 (MQ=255)                                  cGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctg                                                              >  1:662470/1‑100 (MQ=255)                                   gTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgc                                                             >  1:2342716/1‑100 (MQ=255)                                    tttCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCgctgct                                                            <  2:202881/100‑1 (MQ=255)                                     ttCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAAATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGctgctg                                                           >  2:275203/1‑100 (MQ=255)                                        cccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                        >  1:2336245/1‑100 (MQ=255)                                        cccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGaaa                                                        >  2:108434/1‑100 (MQ=255)                                         gcTAGCAAAATCCGCCGTTTCGACCTCAGCACTGGTGCTACCATCTCCGGTAAGATTTGTCCGCCGAAAGAAGGTGACCGCTATTTTCCCctgctgtca                                                        >  2:1418875/2‑96 (MQ=255)                                         ccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAg                                                       >  2:186327/1‑100 (MQ=255)                                         ccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAg                                                       <  1:103546/100‑1 (MQ=255)                                         ccTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAg                                                       <  2:726192/100‑1 (MQ=255)                                           tAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGtt                                                     >  1:648234/1‑100 (MQ=255)                                           tAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGtt                                                     >  2:762220/1‑100 (MQ=255)                                           tAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGtt                                                     >  2:3025/1‑100 (MQ=255)                                              ccAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTCAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTCCGCTGCTCAAAg                                                       >  2:1859013/1‑95 (MQ=255)                                               cAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                 >  1:644212/1‑100 (MQ=255)                                               cAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                 <  2:1704163/100‑1 (MQ=255)                                               cAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACg                                                 >  2:1522583/1‑100 (MQ=255)                                                aaaTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGa                                                <  2:2124065/100‑1 (MQ=255)                                                  aTCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAg                                              <  1:903554/100‑1 (MQ=255)                                                   tCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGt                                             <  1:295550/100‑1 (MQ=255)                                                    ctgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGtt                                            <  2:1114530/98‑1 (MQ=255)                                                     cgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTa                                           >  1:1176908/1‑100 (MQ=255)                                                     cgccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGACGCAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTa                                           >  2:2246449/1‑100 (MQ=255)                                                      gccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTag                                          >  1:1114535/1‑99 (MQ=255)                                                       ccgTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAAc                                         <  2:898805/100‑1 (MQ=255)                                                         gTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt                                       >  2:744495/1‑100 (MQ=255)                                                         gTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACtt                                       >  1:2151420/1‑100 (MQ=255)                                                           ttCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCg                                     <  1:2342976/100‑1 (MQ=255)                                                            tCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACTTCGTTAACtttca                                    <  1:954045/100‑4 (MQ=255)                                                             cAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                   <  2:1928549/100‑1 (MQ=255)                                                             cAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                   <  1:557383/100‑1 (MQ=255)                                                             cAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                   <  1:1299278/100‑1 (MQ=255)                                                             cAACCTCCCCACTGGTGATACCATCTCTGTTAAGAGTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGAc                                   <  2:447501/100‑1 (MQ=255)                                                              aaCCTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGCAAGTTAACGAAGTTAACTTAGACa                                  >  2:2246753/1‑100 (MQ=255)                                                                ccTCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACaaa                                >  2:493304/1‑100 (MQ=255)                                                                  tCCGCTCTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                              <  2:1837458/100‑1 (MQ=255)                                                                  tCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAcc                              <  2:2342695/100‑1 (MQ=255)                                                                  tCCGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAAc                               >  2:1504485/1‑99 (MQ=255)                                                                   ccGCACTGGTGATACCGTCTCTGGTAAGATTTTCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGCTAACTTCGACAACCCt                             >  2:8355/1‑100 (MQ=255)                                                                   ccGCACTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCt                             >  1:2043032/1‑100 (MQ=255)                                                                        cTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaaa                        >  2:825306/1‑100 (MQ=255)                                                                        cTGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGaaaa                        <  1:1833653/100‑1 (MQ=255)                                                                         tGGTGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTCACTTCGACAAACGTGaaa                         >  2:2291930/1‑98 (MQ=255)                                                                          ggTGATACCATCTCTGGTTAGATTTGCACTCCGTATTTAGTTGAACGCTATTTTGAGCTGCTGATACTTAACGAAGTTACCTTTGACAAATCt                             >  2:2194872/1‑93 (MQ=255)                                                                            tGATACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGcc                    >  1:1230338/1‑100 (MQ=255)                                                                              aTACCATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCg                  <  1:689706/100‑1 (MQ=255)                                                                                 ccATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaa               >  1:1698519/1‑100 (MQ=255)                                                                                  cATCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaac              >  1:1355023/1‑100 (MQ=255)                                                                                   aTCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaaca             >  2:842102/1‑100 (MQ=255)                                                                                   aTCTCTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGcaaca             >  2:1542200/1‑100 (MQ=255)                                                                                    tctcTGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCTCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCTcaacaa            >  2:1031999/1‑100 (MQ=255)                                                                                       ctGGTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAAt         <  1:1988797/100‑1 (MQ=255)                                                                                        tGGTAAGATTTGCCCGCCGAGAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAGATc        >  1:264144/1‑100 (MQ=255)                                                                                         ggTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc       <  1:1275302/100‑1 (MQ=255)                                                                                          gTAAGATTTGCCCGCCGAAAGACGGTGGATGCTAAGTTGCGCTGCTGAAAGATGACGGAGTAGACGTGTACAAACCGGAAGACGCCCGaaaataaa          >  2:286141/1‑91 (MQ=255)                                                                                          gTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGGTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATc        >  2:146917/1‑98 (MQ=255)                                                                                          gTAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGCTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATcc       <  2:340205/99‑1 (MQ=255)                                                                                           tAAGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctc     >  2:1778559/1‑100 (MQ=255)                                                                                            aaGATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCctct    >  1:255866/1‑100 (MQ=255)                                                                                              gATTTGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCttt  >  2:1719274/1‑100 (MQ=255)                                                                                                  |                                                                                                CAGCTCCTACCTCGCCGGTCCTGATGACATCTACGTTTCCCCTAGCCAAATCCGCCGTTTCAACCTCCGCACTGGTGATACCATCTCTGGTAAGATTCGCCCGCCGAAAGAAGGTGAACGCTATTTTGCGCTGCTGAAAGTTAACGAAGTTAACTTCGACAAACCTGAAAACGCCCGCAACAAAATCCTCTTT  >  NC_000913/3966623‑3966815

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑