Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,669,524 G→T 5.1% A23A (GCC→GCA ynfL ← putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,669,5240GT 5.1% 262.3 / 3.1 78A23A (GCC→GCAynfLputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (40/34);  new base T (2/2);  total (42/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

CTATCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTCCTACCCTATGTATAAGCCTGA  >  NC_000913/1669384‑1669669
                                                                                                                                            |                                                                                                                                                 
cTATCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGCGCCGCTGCCCGcctg                                                                                                                                          <  2:997510/150‑3 (MQ=255)
cTATCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGCGCCGCTGCCCGcctg                                                                                                                                          <  1:183435/150‑3 (MQ=255)
 tATCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGa                                                                                                                                         <  1:379799/150‑1 (MQ=255)
  aTCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACAGTGCGATTGGCTCGGGCCAGCAGTCGCGCACCAATTTGGTGCTCCAGCGCCGGAATCTGCTGACTTAGCGGCGGTTACGAAATATTCAGGCGGGCGCGCCGCGgcccgcc                                                                                                                                       <  1:1090976/150‑5 (MQ=255)
    cTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAAt                                                                                                                                      >  1:1418719/1‑150 (MQ=255)
     tGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATg                                                                                                                                     >  1:1494823/1‑150 (MQ=255)
      gCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGc                                                                                                                                    <  2:596793/150‑1 (MQ=255)
       cAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCa                                                                                                                                   >  1:979314/1‑150 (MQ=255)
         aGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGc                                                                                                                                 >  1:53943/1‑150 (MQ=255)
                     ccTGCTGCCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCgcaaca                                                                                                                     >  2:1230523/1‑150 (MQ=255)
                       tgctgcCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCgcaacagc                                                                                                                   >  1:385312/1‑150 (MQ=255)
                        gctgcCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCgcaacagca                                                                                                                  <  2:963919/150‑1 (MQ=255)
                           gcCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCgcaacagcaaca                                                                                                               <  1:1457889/150‑1 (MQ=255)
                             cGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACaaa                                                                                                             >  2:120185/1‑150 (MQ=255)
                              gTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAg                                                                                                            <  1:823260/150‑1 (MQ=255)
                                  gCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAAc                                                                                                        <  1:543123/150‑1 (MQ=255)
                                  gCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAAc                                                                                                        <  2:1348411/150‑1 (MQ=255)
                                     aTACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCa                                                                                                     >  1:1022364/1‑150 (MQ=255)
                                     aTACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCa                                                                                                     <  1:18230/150‑1 (MQ=255)
                                        cacTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGAt                                                                                                  >  1:584121/1‑150 (MQ=255)
                                        cacTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGAt                                                                                                  >  1:28641/1‑150 (MQ=255)
                                          cTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGa                                                                                                >  1:883623/1‑150 (MQ=255)
                                           tGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGAc                                                                                               >  1:734118/1‑150 (MQ=255)
                                               aTTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAg                                                                                           <  1:172006/150‑1 (MQ=255)
                                                   gTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCa                                                                                                                                   <  2:29620/106‑1 (MQ=255)
                                                   gTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCa                                                                                                                                   >  1:29620/1‑106 (MQ=255)
                                                   gTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCa                                                                                       >  1:1455233/1‑150 (MQ=255)
                                                    ttCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCaa                                                                                      >  1:162000/1‑150 (MQ=255)
                                                    gtCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCaa                                                                                      <  1:1436628/149‑1 (MQ=255)
                                                     tCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAat                                                                                     >  1:671088/1‑150 (MQ=255)
                                                      cGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAata                                                                                    >  1:1305785/1‑150 (MQ=255)
                                                         gCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTc                                                                                 >  2:276005/1‑150 (MQ=255)
                                                              cagTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATAtt                                                                            >  1:17530/1‑150 (MQ=255)
                                                                gTCGGTCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTa                                                                          <  1:1245316/150‑1 (MQ=255)
                                                                  cGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTaaa                                                                        >  2:1169976/1‑150 (MQ=255)
                                                                    ggCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAAc                                                                      >  2:830104/1‑150 (MQ=255)
                                                                      cACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACAt                                                                    >  2:983033/1‑150 (MQ=255)
                                                                      cACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATGAATATTGaacttcg                                                                    <  1:138574/150‑6 (MQ=255)
                                                                              tgttgCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTg                                                            <  2:327288/150‑1 (MQ=255)
                                                                              tgttgCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTg                                                            <  1:1044015/150‑1 (MQ=255)
                                                                              tgttgCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTg                                                            <  1:1123289/150‑1 (MQ=255)
                                                                                ttgCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTgag                                                          <  2:1418719/150‑1 (MQ=255)
                                                                                ttgCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTgag                                                          <  1:1496519/150‑1 (MQ=255)
                                                                                  gCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGat                                                        <  2:734118/150‑1 (MQ=255)
                                                                                  gCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGat                                                        >  1:1023753/1‑150 (MQ=255)
                                                                                    tCCAGGGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGatta                                                      <  1:1058496/150‑1 (MQ=255)
                                                                                     ccAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGattat                                                     >  1:1401017/1‑150 (MQ=255)
                                                                                     ccAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGattat                                                     >  2:417898/1‑150 (MQ=255)
                                                                                         cgcCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAat                                                 <  2:1494823/150‑1 (MQ=255)
                                                                                             tGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAatatat                                             >  1:1151975/1‑150 (MQ=255)
                                                                                                  ctgctgACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGAc                                        <  2:1022364/150‑1 (MQ=255)
                                                                                                  ctgctgACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGAc                                        <  1:120185/150‑1 (MQ=255)
                                                                                                   tgctgACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACa                                       <  1:1519066/150‑1 (MQ=255)
                                                                                                   tgctgACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACa                                       >  1:104660/1‑150 (MQ=255)
                                                                                                   tgctgACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACa                                       >  2:1477195/1‑150 (MQ=255)
                                                                                                      tgACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGaa                                    <  1:1230523/150‑1 (MQ=255)
                                                                                                      tgACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGaa                                    >  2:508563/1‑150 (MQ=255)
                                                                                                       gACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAAc                                   <  2:1151975/150‑1 (MQ=255)
                                                                                                            aGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTc                              >  2:803665/1‑150 (MQ=255)
                                                                                                               ggcggTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGAt                           >  2:1503838/1‑150 (MQ=255)
                                                                                                                 cggTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGAttt                         <  2:446271/150‑1 (MQ=255)
                                                                                                                    ttGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTcc                      <  1:831775/150‑1 (MQ=255)
                                                                                                                       gctgcatttCGGGCGCGCCGCTGCCCGCCTGAAATGCAGCTCTTCCGCCACAGCAACAAAGTAACGCAGATGACGAAGGTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTCCTAc                   >  1:249692/8‑150 (MQ=255)
                                                                                                                       gctgcatttCGGGCGCGCCGCTGCCCGCCTGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTCCTAc                   >  1:63694/8‑150 (MQ=255)
                                                                                                                       cGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATg                                                                                                 <  2:941677/72‑1 (MQ=255)
                                                                                                                       cGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATg                                                                                                 >  1:941677/1‑72 (MQ=255)
                                                                                                                        gAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTCCTAcc                  >  2:1365244/1‑150 (MQ=255)
                                                                                                                        gAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTCCTAcc                  >  2:631437/1‑150 (MQ=255)
                                                                                                                             atTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTCCTACCCTATg             >  2:1148147/1‑150 (MQ=255)
                                                                                                                              ttCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTCCTACCCTATGt            <  1:830104/150‑1 (MQ=255)
                                                                                                                                     gggCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTCCTACCCTATGTATAAGcc     <  1:276005/150‑1 (MQ=255)
                                                                                                                                      ggCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGAttt                         <  2:981111/129‑1 (MQ=255)
                                                                                                                                      ggCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGAttt                         >  1:981111/1‑129 (MQ=255)
                                                                                                                                      ggCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTCCTACCCTATGTATAAGCCt    >  2:393257/1‑150 (MQ=255)
                                                                                                                                      ggCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTCCTACCCTATGTATAAGCCt    >  1:926021/1‑150 (MQ=255)
                                                                                                                                       gCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTCCTACCCTATGTATAAGCCTg   >  2:943737/1‑150 (MQ=255)
                                                                                                                                       gCAGCGGCGCGCCCGAAATGCAGCCCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTCCTACCCTATGTATAAGCCTg   >  1:500078/1‑150 (MQ=255)
                                                                                                                                        cAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTCCTACCCTATGTATAAGCCTGa  <  1:85711/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
CTATCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTCCTACCCTATGTATAAGCCTGA  >  NC_000913/1669384‑1669669

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: