Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,663,395 G→A 6.0% G136R (GGG→AGG)  ynfG → oxidoreductase, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,663,3950GA 6.0% 207.5 / 3.1 66G136R (GGG→AGG) ynfGoxidoreductase, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base G (33/29);  new base A (2/2);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

CGTGAAGATGGCTTTGTGGTGGTTGATGAAGATGTTTGTATCGGCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATCGTTATCAAACC  >  NC_000913/1663248‑1663543
                                                                                                                                                   |                                                                                                                                                    
cGTGAAGATGGCTTTGTGGTGGTTGATGAAGATGTTTGTATCGGCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAggg                                                                                                                                                     <  1:400566/149‑1 (MQ=255)
   gAAGATGGCTTTGTGGTGGTTGATGAAGATGTTTGTATCGGCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGa                                                                                                                                                   >  5:157357/1‑148 (MQ=255)
    aaGATGGCTTTGTGGTGGTTGATGAAGATGTTTGTATCGGCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGaaa                                                                                                                                                 <  1:161943/149‑1 (MQ=255)
           cattGTGGTGGTTGATGAAGATGTTTGTATCGGCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACCAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAAccc                                                                                                                                           >  8:397008/3‑148 (MQ=255)
            tttGTGGTGGTTGATGAAGATGTTTGTATCGGCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTCGTTATTCGCGCGTCGCCGAggg                                                                                                                                                     >  1:149346/1‑137 (MQ=255)
            tttGTGGTGGTTGATGAAGATGTTTGTATCGGCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTCGTTATTCGCGCGTCGCCGAggg                                                                                                                                                     <  2:149346/137‑1 (MQ=255)
               gtggtgGTTGATGAAGATGTTTGTATCGGCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCatat                                                                                                                                       <  7:397008/148‑1 (MQ=255)
               ggggtggtTGATGACGATGTGTGTATCGGCTGCCGGTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAAGCCatat                                                                                                                                       <  6:331415/146‑1 (MQ=255)
                   tggtTGATGAAGATGTTTGTATCGGCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTc                                                                                                                                  <  7:126599/149‑1 (MQ=255)
                         atgaAGATGTTTGTATCGGCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATcc                                                                                                                            >  3:151491/1‑149 (MQ=255)
                               aTGTTTGTATCGGCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCg                                                                                                                      >  8:312885/1‑149 (MQ=255)
                                      tATCGGCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTg                                                                                                               <  6:347243/149‑1 (MQ=255)
                                       aTCGGCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCatat                                                                                                                                       >  5:32843/1‑124 (MQ=255)
                                       aTCGGCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCatat                                                                                                                                       <  6:32843/124‑1 (MQ=255)
                                           gCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACtt                                                                                                          >  3:226085/1‑149 (MQ=255)
                                           gCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGAGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACtt                                                                                                          >  1:109677/1‑149 (MQ=255)
                                           gCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGAGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACtt                                                                                                          >  1:104243/1‑149 (MQ=255)
                                               ccGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGt                                                                                                      >  6:223986/1‑149 (MQ=255)
                                                     cTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGTGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTc                                                                                                >  4:356557/1‑149 (MQ=255)
                                                        ccACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATACTGCCCGCTGCGTGCACTTGAGTTCGGTCCaa                                                                                             >  5:90724/1‑149 (MQ=255)
                                                          acaTGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAAtt                                                                                           >  2:106600/1‑149 (MQ=255)
                                                          acaTGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAAtt                                                                                           >  5:207665/1‑149 (MQ=255)
                                                               gCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGaa                                                                                        >  1:67267/1‑147 (MQ=255)
                                                                cTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGaa                                                                                        >  8:364942/1‑146 (MQ=255)
                                                                 ttGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGc                                                                                    >  1:163014/1‑149 (MQ=255)
                                                                 ttGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTAATCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGc                                                                                    >  7:309851/1‑149 (MQ=255)
                                                                      cGTACGGAGCGCCACAGTACAATGCTGAAATAGGGCTCATGACGATGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTgcgc                                                                               <  6:157357/149‑1 (MQ=255)
                                                                           ggCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGaa                                                                          <  5:223986/149‑1 (MQ=255)
                                                                           ggCGCGCCACAGTACAATGCGGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTTTTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGaa                                                                          <  8:309851/149‑1 (MQ=255)
                                                                            gcgcgcCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGt                                                                                                                <  7:257191/110‑1 (MQ=255)
                                                                            gcgcgcCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGt                                                                                                                >  8:257191/1‑110 (MQ=255)
                                                                                 ccACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACgg                                                                    <  2:163014/149‑1 (MQ=255)
                                                                                  cacaGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATCGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACgg                                                                    >  6:118288/1‑148 (MQ=255)
                                                                                        aCAATGCGGAAAACGGGCACCTGACGACGTGCGATGGTCGTTATCCGCGCGCCGCCGACGCGAACCAACCCATATGGGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCtt                                                             <  4:251813/149‑1 (MQ=255)
                                                                                         cAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCtt                                                             <  2:67267/148‑1 (MQ=255)
                                                                                         cAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTtg                                                            <  6:235537/149‑1 (MQ=255)
                                                                                          aaTGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTtgc                                                           <  1:106600/149‑1 (MQ=255)
                                                                                              cTGAAAAATGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTGGCCGAGGCGAAACAAACGATATGTGTCGACTCCTGCCCCCTGCGTTCACTTGAGTTCGGTCCAATTGAATAGCTGCGCCAGAAGCACGGCACTCTtgctgc                                                        <  1:201197/148‑1 (MQ=255)
                                                                                               tGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCAAATGTGTCGAATCCTGCCAGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTtgctgcag                                                      >  5:57569/1‑147 (MQ=255)
                                                                                                gAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGt                                                     >  3:119586/1‑149 (MQ=255)
                                                                                                   aaaGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTcgc                                                  <  4:265367/149‑1 (MQ=255)
                                                                                                          acaTGACGAAGTGCGATGGTTGTTATACGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTg                                           >  2:373243/1‑149 (MQ=255)
                                                                                                             tGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAGGCACGGCACTCTTGCTGCCGTCGCGCCACTGCcg                                        >  4:299015/1‑149 (MQ=255)
                                                                                                              gACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCcgc                                       >  7:360177/1‑149 (MQ=255)
                                                                                                               aCGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCcgcg                                      >  2:20792/1‑149 (MQ=255)
                                                                                                                 gAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCcgcgcg                                    >  2:197762/1‑149 (MQ=255)
                                                                                                                 gAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCcgcgcg                                    >  7:371893/1‑149 (MQ=255)
                                                                                                                   aGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCcgcgcgcg                                  <  8:374959/149‑1 (MQ=255)
                                                                                                                   aGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCcgcgcgcg                                  >  7:374959/1‑149 (MQ=255)
                                                                                                                        gATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCAttt                             >  8:280514/1‑149 (MQ=255)
                                                                                                                         aTGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCAttt                             >  4:236627/1‑148 (MQ=255)
                                                                                                                          tGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAGGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTca                           <  3:299015/149‑1 (MQ=255)
                                                                                                                             ttgttATTCGCGCGTCGCCGAGGGGAAACAACCCatat                                                                                                                                       <  2:299041/38‑1 (MQ=255)
                                                                                                                             ttgttATTCGCGCGTCGCCGAGGGGAAACAACCCatat                                                                                                                                       >  1:299041/1‑38 (MQ=255)
                                                                                                                                ttATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAAc                     >  4:19698/1‑149 (MQ=255)
                                                                                                                                ttATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAAc                     >  1:257095/1‑149 (MQ=255)
                                                                                                                                     cgcgcgTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCaa                 <  7:280514/148‑1 (MQ=255)
                                                                                                                                     cgcgcgTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCaa                 <  3:19698/148‑1 (MQ=255)
                                                                                                                                      gcgcgTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAata               <  7:312885/149‑1 (MQ=255)
                                                                                                                                        gcgTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATc             <  8:360177/149‑1 (MQ=255)
                                                                                                                                        gcgTCGCCGAGAGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATc             <  2:109677/149‑1 (MQ=255)
                                                                                                                                        gcgTCGCCGAGAGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATc             <  2:104243/149‑1 (MQ=255)
                                                                                                                                         cgTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCAAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATCg            <  7:364942/149‑1 (MQ=255)
                                                                                                                                          gTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATCGt           >  7:321362/1‑149 (MQ=255)
                                                                                                                                           tCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATCGtt          <  1:20792/149‑1 (MQ=255)
                                                                                                                                           tCGCCGAGGGGAAACAACCCATATGTGGCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGCAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATCGtt          <  4:226085/149‑1 (MQ=255)
                                                                                                                                             gCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATCGTTAt        <  5:118288/149‑1 (MQ=255)
                                                                                                                                             gCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATCGTTAt        >  2:210879/1‑149 (MQ=255)
                                                                                                                                               cGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATCGTTATCa      <  1:197762/149‑1 (MQ=255)
                                                                                                                                                  ggggAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATCGTTATCAAAc   <  6:207665/149‑1 (MQ=255)
                                                                                                                                                   gggAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATCGTTATCAAAcc  >  2:219920/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CGTGAAGATGGCTTTGTGGTGGTTGATGAAGATGTTTGTATCGGCTGCCGCTACTGCCACATGGCTTGTCCGTACGGCGCGCCACAGTACAATGCTGAAAAAGGGCACATGACGAAGTGCGATGGTTGTTATTCGCGCGTCGCCGAGGGGAAACAACCCATATGTGTCGAATCCTGCCCGCTGCGTGCACTTGAGTTCGGTCCAATTGAAGAGCTGCGCCAGAAGCACGGCACTCTTGCTGCCGTCGCGCCACTGCCGCGCGCGCATTTCACAAAACCCAATATCGTTATCAAACC  >  NC_000913/1663248‑1663543

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