Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,669,521 A→C 5.1% A24A (GCT→GCG 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,5210AC 5.1% 247.5 / 4.8 79A24A (GCT→GCGynfLputative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (41/33);  new base C (2/2);  total (43/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 = 6.06e-01

GATTTGCCGACTATCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTCCTACCCTATGTATAAGCCTGA  >  NC_000913/1669374‑1669669
                                                                                                                                                   |                                                                                                                                                    
gATTTGCCGACTATCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAgc                                                                                                                                                    >  1:52572/1‑150 (MQ=255)
          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                                                                                                        <  2:1348411/150‑1 (MQ=255)
                                            gCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAAc                                                                                                        <  1:543123/150‑1 (MQ=255)
                                               aTACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCa                                                                                                     <  1:18230/150‑1 (MQ=255)
                                               aTACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCa                                                                                                     >  1:1022364/1‑150 (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                                                                                                                                   >  1:29620/1‑106 (MQ=255)
                                                             gTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCa                                                                                                                                   <  2:29620/106‑1 (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                                                          <  1:1496519/150‑1 (MQ=255)
                                                                                          ttgCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTgag                                                          <  2:1418719/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                                                     >  2:417898/1‑150 (MQ=255)
                                                                                               ccAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGattat                                                     >  1:1401017/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                         >  1:981111/1‑129 (MQ=255)
                                                                                                                                                ggCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGAttt                         <  2:981111/129‑1 (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)
                                                                                                                                                   |                                                                                                                                                    
GATTTGCCGACTATCTGCAAGAAACTGTTTTCCTGCTGCCGTCAGCAATACACTGCGATTGGTTCGTGCCAGCAGTCGGGCACCAATTTGTTGCTCCAGCGCCTGAATCTGCTGACTTAGCGGCGGTTGCGAAATATTCAGGCGGGCAGCGGCGCGCCCGAAATGCAGCTCTTCCGCAACAGCAACAAAGTAACGCAGATGACGAAGTTCAATATTCATATTTAAAACATCTTATTTGAGATTATTAATATATTAGACAGAACAATTCGATTTTCCTACCCTATGTATAAGCCTGA  >  NC_000913/1669374‑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: