Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,679,641 A→T 8.4% N29Y (AAT→TAT)  ydgI → putative arginine/ornithine antiporter transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,679,6410AT 8.4% 274.7 / 8.8 93N29Y (AAT→TAT) ydgIputative arginine/ornithine antiporter transporter
Reads supporting (aligned to +/- strand):  ref base A (42/43);  new base T (6/2);  total (48/45)
Fisher's exact test for biased strand distribution p-value = 2.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGC  >  NC_000913/1679492‑1679780
                                                                                                                                                     |                                                                                                                                           
cTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGCCTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCaaa                                                                                                                                             >  2:1047318/1‑150 (MQ=255)
 tCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCaaaa                                                                                                                                            >  2:1586854/1‑150 (MQ=255)
     tttttATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGGTGATATTAAGCTCAATGATGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATg                                                                                                                                        >  2:36578/1‑150 (MQ=255)
       tttATGAGTTTTGCCCTAATATTTATTATTATTATCACGTTTTAATCACTGGATATCGAAGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAAAGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATggc                                                                                                                                      <  2:1018804/150‑1 (MQ=255)
         tATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATggcgg                                                                                                                                    >  1:536254/1‑150 (MQ=255)
          ttGCTTTTTGCCCTTATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCTATGGAAAAGAAACTGGGACTGAGCGGACTTACGGCTCTGGTATTAAGCTCTATGCTGGGCGCGTTTTTTTTCAGTCTGCCGCAAAATATggcggc                                                                                                                                   <  1:1687398/149‑1 (MQ=255)
              gTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCCATATTTTGCGGCAgac                                                                                                                               <  1:1796677/150‑3 (MQ=255)
               ttttGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTg                                                                                                                              <  1:13482/150‑1 (MQ=255)
                tttGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCCATATTTTGCGGCAgactg                                                                                                                             <  2:1591606/150‑5 (MQ=255)
                 ttGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGcc                                                                                                                            <  2:1280606/150‑1 (MQ=255)
                   gCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAg                                                                                                                          >  1:695198/1‑150 (MQ=255)
                   gCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAg                                                                                                                          <  1:1381336/150‑1 (MQ=255)
                    cccTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGc                                                                                                                         <  1:703812/150‑1 (MQ=255)
                          tatttattCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGgca                                                                                                                   >  1:850292/1‑150 (MQ=255)
                          tatttattCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGAGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAATATATGGCGGCAGTTGCCAGCCCGgcc                                                                                                                   >  1:1256955/1‑149 (MQ=255)
                           atttattCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGgcag                                                                                                                  >  2:446595/1‑150 (MQ=255)
                           atttattCTTCTTATCACGTTTAAATCACTGGATATCGATCGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGgcag                                                                                                                  <  1:378933/150‑1 (MQ=255)
                            tttattCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGgcagc                                                                                                                 <  1:695327/150‑1 (MQ=255)
                                 tcttctTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTGAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGc                                                                                                            >  2:1110714/1‑150 (MQ=255)
                                 tcttctTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGc                                                                                                            >  1:856129/1‑150 (MQ=255)
                                  cttcttATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCt                                                                                                           >  2:342923/1‑150 (MQ=255)
                                   ttcttATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTc                                                                                                          <  2:1104093/150‑1 (MQ=255)
                                        aTCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCgg                                                                                                     >  2:1652207/1‑150 (MQ=255)
                                        aTCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCgg                                                                                                     >  2:894774/1‑150 (MQ=255)
                                        aTCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCgg                                                                                                     >  1:219728/1‑150 (MQ=255)
                                          cACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCt                                                                                                   >  2:1104820/1‑150 (MQ=255)
                                           aCGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTg                                                                                                  <  2:1351812/150‑1 (MQ=255)
                                            cGTTTTAATCACTGGATATCGATGGATAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTgg                                                                                                 <  2:1512037/150‑1 (MQ=255)
                                              ttttAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTgggg                                                                                               <  1:217548/150‑1 (MQ=255)
                                                ttAATCACTGGATATCGATGGATAAGAACCTGGGACTGCGCGCACTCACCGCGCTGGGATTAAGCTCAATGGTGGGCGCGGGTGTTTTCACTCTGCCGCAAAATATGGCGGCAGTTTCCAGCCCGGCCGCAGGGCTAAGCGGCGGGGGTa                                                                                             >  2:1151713/1‑150 (MQ=255)
                                                  aaTCACTGGATATCGATGGAAAAGAAACTGGGACTGCGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTAtt                                                                                           >  2:1787224/1‑150 (MQ=255)
                                                   aTCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTa                                                                                          >  1:1544782/1‑150 (MQ=255)
                                                       cTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTgg                                                                                      <  1:1635096/150‑1 (MQ=255)
                                                         ggATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGgcg                                                                                    <  2:446676/150‑1 (MQ=255)
                                                         ggATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGgcg                                                                                    >  1:1546562/1‑150 (MQ=255)
                                                             atCGATGGACAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTgg                                                                                <  1:104742/150‑1 (MQ=255)
                                                             atCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTgg                                                                                <  1:530222/150‑1 (MQ=255)
                                                                gATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCAGCGGCTGGGGTATTACTGGCGCTGGCAt                                                                             >  2:251764/1‑150 (MQ=255)
                                                                gATGGAAAAGAAACTGGGACTGAGCGAACACACCGCGCTGCTATAAAGCTCAAAGCTGGGCGCGGGTGTTTTCAGTCTCCCGCAAAATAAGGCGGCACTTGCCAGCCCGGCAGCACTGCTCATCGCCTGGGGTATTACTGGCGCTGTCaa                                                                             >  2:1804404/1‑149 (MQ=255)
                                                                 aTGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCAtt                                                                            >  1:1703235/1‑150 (MQ=255)
                                                                    gAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATtata                                                                         >  1:541997/1‑150 (MQ=255)
                                                                    gAAAAGAAACTGGGACTGAGCGCACTCACCGCGCAGGTATTAAGCTCTATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCATCCCGGCAGCCCTGCTCATCGTCTGGGGTATTACTGGCGCTGGCAGTTTa                                                                         >  2:1640121/1‑150 (MQ=255)
                                                                     aacaGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTAt                                                                        <  1:383259/147‑1 (MQ=255)
                                                                      aaaGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGACAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTAtt                                                                       <  2:568622/150‑1 (MQ=255)
                                                                         gAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCt                                                                    >  1:780268/1‑150 (MQ=255)
                                                                           aaCTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTgg                                                                  >  1:1436646/1‑150 (MQ=255)
                                                                             cTGGGACTGAGCGCACTCACCTCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGcc                                                                <  2:1446369/150‑1 (MQ=255)
                                                                             cTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGcc                                                                <  2:490548/150‑1 (MQ=255)
                                                                             cTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGcc                                                                <  2:1533528/150‑1 (MQ=255)
                                                                              tGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCTCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGACt                                                               >  1:1510713/1‑150 (MQ=255)
                                                                                 gACTGCGCGCACTCACCGCGCTGGGATTAAGCTCAATGCTGGCCGCGGGTGGTTTCAGGCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGAATTACTGGCGCTGGCATTTTATTGCTGGCCTTTg                                                            <  2:277100/150‑1 (MQ=255)
                                                                                 gACTGAGCGCACTCACCGCGCTTGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTg                                                            <  1:573331/150‑1 (MQ=255)
                                                                                 gACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTg                                                            >  1:1587389/1‑150 (MQ=255)
                                                                                   cTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGcc                                                          <  1:1367652/150‑1 (MQ=255)
                                                                                   cTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGcc                                                          <  1:1257721/150‑1 (MQ=255)
                                                                                     gAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCAt                                                        >  2:307076/1‑150 (MQ=255)
                                                                                     gAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGGGTTTTCAGTCTGACGCAAAATATGTCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCAGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCAt                                                        >  2:193336/1‑150 (MQ=255)
                                                                                        cgcACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGGGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCt                                                     <  2:1423654/150‑1 (MQ=255)
                                                                                           aCTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTACCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGAt                                                  <  2:563695/150‑1 (MQ=255)
                                                                                            cTCACCTCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGTGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATc                                                 <  1:1529387/150‑1 (MQ=255)
                                                                                            cTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCCCTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATc                                                 >  1:1145031/1‑150 (MQ=255)
                                                                                              cACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCt                                               >  2:539201/1‑150 (MQ=255)
                                                                                                ccGCGCTGGTATAAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCa                                             <  2:1659946/150‑1 (MQ=255)
                                                                                                     cTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACgcgc                                        >  1:518384/1‑150 (MQ=255)
                                                                                                       ggTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATACGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCAt                                      <  2:246527/150‑1 (MQ=255)
                                                                                                         tATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCCCGCGCATTc                                    >  1:848155/1‑150 (MQ=255)
                                                                                                         tATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTc                                    <  2:1205503/150‑1 (MQ=255)
                                                                                                         tATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTc                                    >  1:413702/1‑150 (MQ=255)
                                                                                                          aTTAAGCTCAAGGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCg                                   <  2:1677108/150‑1 (MQ=255)
                                                                                                              aGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTccc                               >  2:1488280/1‑150 (MQ=255)
                                                                                                              aGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTccc                               >  1:1661835/1‑150 (MQ=255)
                                                                                                               gCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCg                              >  2:680620/1‑150 (MQ=255)
                                                                                                               gCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCg                              >  1:1577635/1‑150 (MQ=255)
                                                                                                               gCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGAATTTTATTGCTGGCCATTGCCATGCTGATCCTCACGCGCATTCGTcccc                              >  2:1513914/1‑149 (MQ=255)
                                                                                                                 tCACTGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGaa                            <  2:1122762/150‑1 (MQ=255)
                                                                                                                 tCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGGTCCTCACGCGCATTCGTCCCGaa                            >  1:1318966/1‑150 (MQ=255)
                                                                                                                  cAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAAc                           <  1:770219/150‑1 (MQ=255)
                                                                                                                    aTGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACtt                         >  1:654583/1‑150 (MQ=255)
                                                                                                                      gCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGa                       <  2:1861986/150‑1 (MQ=255)
                                                                                                                       cTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAc                      <  1:1277293/150‑1 (MQ=255)
                                                                                                                       cTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAc                      <  1:219435/150‑1 (MQ=255)
                                                                                                                       cTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAc                      <  1:1844211/150‑1 (MQ=255)
                                                                                                                            cgcgGGTGTTTTCAGTCAGCCGCAAAATATGGCGGCTGTTGCCAGCCCGGCAGCACTGCTGATCGGCTGTGGTATTACTGGCGCTGGCATTTTTTTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGAcggcgg                 <  1:1565276/150‑1 (MQ=255)
                                                                                                                             gcgGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGt                <  2:1770252/150‑1 (MQ=255)
                                                                                                                                  tgtTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCtt           >  1:510133/1‑150 (MQ=255)
                                                                                                                                     tttCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTAc        <  1:1652207/150‑1 (MQ=255)
                                                                                                                                      ttCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTAcc       <  2:850292/150‑1 (MQ=255)
                                                                                                                                      ttCAGGCTGCCGAAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTAcc       <  1:1682199/150‑1 (MQ=255)
                                                                                                                                       tCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCCCTGCTCATCGGCGGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTCCCATGCTGATCCTCACGCGCCTTCGTCCCGAACTTGACGGTGGTATCTTTACCt      >  2:949565/1‑150 (MQ=255)
                                                                                                                                       ggcaacTGCCGCCATATTTTGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCt      >  1:999022/6‑150 (MQ=255)
                                                                                                                                         caacTGCCGCCATATTTTGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTAt    >  2:1002691/4‑150 (MQ=255)
                                                                                                                                         aGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGTTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTAt    >  1:874642/1‑150 (MQ=255)
                                                                                                                                          aacTGCCGCCATATTTTGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATg   >  1:497816/3‑150 (MQ=255)
                                                                                                                                           acTGCCGCCATATTTTGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGc  >  1:715053/2‑150 (MQ=255)
                                                                                                                                           acTGCCGCCATATTTTGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGAACCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGc  >  2:1366852/2‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                           
CTCCATTTTTATGCGTTTTGCCCTAATATTTATTCTTCTTATCACGTTTTAATCACTGGATATCGATGGAAAAGAAACTGGGACTGAGCGCACTCACCGCGCTGGTATTAAGCTCAATGCTGGGCGCGGGTGTTTTCAGTCTGCCGCAAAATATGGCGGCAGTTGCCAGCCCGGCAGCACTGCTCATCGGCTGGGGTATTACTGGCGCTGGCATTTTATTGCTGGCCTTTGCCATGCTGATCCTCACGCGCATTCGTCCCGAACTTGACGGCGGTATCTTTACCTATGC  >  NC_000913/1679492‑1679780

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: