Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,811,472 A→T 5.7% L68Q (CTG→CAG)  rpmB ← 50S ribosomal subunit protein L28

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,811,4720AT 5.7% 382.3 / 7.4 118L68Q (CTG→CAG) rpmB50S ribosomal subunit protein L28
Reads supporting (aligned to +/- strand):  ref base A (61/50);  new base T (3/4);  total (64/54)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 118 total aligned reads displayed.

TAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTGGG  >  NC_000913/3811330‑3811619
                                                                                                                                              |                                                                                                                                                   
tAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGt                                                                                                                                              >  1:175810/1‑150 (MQ=255)
 aGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTa                                                                                                                                             >  2:538821/1‑150 (MQ=255)
 aGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTa                                                                                                                                             >  2:1509513/1‑150 (MQ=255)
  gTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCAAGAATACCTTTAGCCAAGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCCGCCAGAACTGTAt                                                                                                                                            >  1:745825/1‑150 (MQ=255)
   tACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCTGGCTGAACTgcgtg                                                                                                                                           <  2:715753/150‑5 (MQ=255)
   tACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCTGGCTGAACTgcgtg                                                                                                                                           <  1:320342/150‑5 (MQ=255)
   tACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATc                                                                                                                                           <  2:647070/150‑1 (MQ=255)
     cGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGa                                                                                                                                         >  2:230893/1‑150 (MQ=255)
     cGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGa                                                                                                                                         <  1:538821/150‑1 (MQ=255)
       tttGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATg                                                                                                                                       >  1:78274/1‑150 (MQ=255)
          gTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCt                                                                                                                                    >  1:1393413/1‑150 (MQ=255)
             cttcGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCTGGCTGAACTgcgtg                                                                                                                                           <  2:1251866/140‑5 (MQ=255)
             cttcGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCtttt                                                                                                                                 >  1:1421344/1‑150 (MQ=255)
             cttcGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCtttt                                                                                                                                 <  2:599654/150‑1 (MQ=255)
             cttcGAACTGCTATAGAAGTAACAAGTAAAATCAGAAGAAACCACCTTGATTTTCTCAACAATAAATTTAGACATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCtttt                                                                                                                                 >  2:719136/1‑150 (MQ=255)
              ttcGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCttttt                                                                                                                                <  1:1080488/150‑1 (MQ=255)
              ttcGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCAAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCttttt                                                                                                                                <  2:172247/150‑1 (MQ=255)
               tcGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCtttttt                                                                                                                               <  2:505155/150‑1 (MQ=255)
                cGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTa                                                                                                                              >  2:617121/1‑150 (MQ=255)
                cGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTa                                                                                                                              >  1:364368/1‑150 (MQ=255)
                    gTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGa                                                                                                                          <  2:1226251/150‑1 (MQ=255)
                      ggTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGAtt                                                                                                                        >  2:289467/1‑150 (MQ=255)
                         ataGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCTGGCTGAACTgcgtg                                                                                                                                           >  2:353577/1‑124 (MQ=255)
                         ataGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCTGGCTGAACTgcgtg                                                                                                                                           <  1:353577/128‑5 (MQ=255)
                          taGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTacac                                                                                                                    <  2:179054/150‑1 (MQ=255)
                           aGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACg                                                                                                                   >  1:1152796/1‑150 (MQ=255)
                           aGAAGTGACCAGTACCAGCAGAAGAAACAAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTACTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACg                                                                                                                   >  2:7411/1‑150 (MQ=255)
                            gAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGc                                                                                                                  <  2:364368/150‑1 (MQ=255)
                            gAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGc                                                                                                                  >  2:622636/1‑150 (MQ=255)
                            gAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGc                                                                                                                  >  1:1267001/1‑150 (MQ=255)
                             aaGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTCCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCa                                                                                                                 >  1:996779/1‑150 (MQ=255)
                              aGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCAt                                                                                                                >  1:1463924/1‑150 (MQ=255)
                                 gACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATAcc                                                                                                             >  1:1275225/1‑150 (MQ=255)
                                    cAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCttt                                                                                                          <  1:1392201/150‑1 (MQ=255)
                                        aCCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCa                                                                                                      <  1:26319/150‑1 (MQ=255)
                                        aCCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCa                                                                                                      >  1:791933/1‑150 (MQ=255)
                                          cagcagAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGa                                                                                                    >  1:76523/1‑150 (MQ=255)
                                           agcagAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGAt                                                                                                   >  2:13274/1‑150 (MQ=255)
                                               gaagaaACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATAcgc                                                                                               >  1:380123/1‑150 (MQ=255)
                                               gaagaaACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATAcgc                                                                                               >  2:147375/1‑150 (MQ=255)
                                               gaagaaACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATAcgc                                                                                               >  1:26892/1‑150 (MQ=255)
                                                 agaaACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACgcgc                                                                                             >  1:730495/1‑150 (MQ=255)
                                                  gaaACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCa                                                                                            >  1:853937/1‑150 (MQ=255)
                                                  gaaACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCa                                                                                            >  2:1107768/1‑150 (MQ=255)
                                                      ccAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGt                                                                                        >  2:783176/1‑150 (MQ=255)
                                                       cAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTg                                                                                       <  2:291827/150‑1 (MQ=255)
                                                         gCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGAc                                                                                     >  1:290067/1‑150 (MQ=255)
                                                          cTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCttt                                                                                                          <  2:690686/128‑1 (MQ=255)
                                                          cTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCttt                                                                                                          >  1:690686/1‑128 (MQ=255)
                                                                 ttCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCt                                                                             <  1:13274/150‑1 (MQ=255)
                                                                     cACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTtctc                                                                         >  2:1305100/1‑150 (MQ=255)
                                                                     cACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTtctc                                                                         >  2:1312791/1‑150 (MQ=255)
                                                                      aCGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCg                                                                        >  1:208433/1‑150 (MQ=255)
                                                                              cTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAAc                                                                <  1:1312791/150‑1 (MQ=255)
                                                                                ttAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAAccc                                                              <  2:485184/150‑1 (MQ=255)
                                                                                 tAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCa                                                             >  1:285740/1‑150 (MQ=255)
                                                                                  aGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAg                                                            <  1:972917/150‑1 (MQ=255)
                                                                                   gCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGa                                                           >  2:1343784/1‑150 (MQ=255)
                                                                                   gCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGa                                                           >  1:1089791/1‑150 (MQ=255)
                                                                                       tGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACg                                                       >  2:279447/1‑150 (MQ=255)
                                                                                          tttatttCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTg                                                    <  1:1430915/150‑1 (MQ=255)
                                                                                           ttatttCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTga                                                   <  2:94458/150‑1 (MQ=255)
                                                                                            tatttCCTCTAAGTACTTAGTACATTTCGCCACGCGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCACATACGCGCAGGGTCACAAAACCCTTCTCGCTCTCAACCCAGAAACGGTgac                                                  >  2:1423670/1‑149 (MQ=255)
                                                                                             atttCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTAAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTgaga                                                 >  1:601248/1‑150 (MQ=255)
                                                                                              tttCCTCTAAGTACTTAGTACTTTTCGCCTCGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATAGGCGCAGGGTGACATAAGGCTTCTCGCTCTCAACCCAGAATGTGTGagag                                                <  2:1204551/150‑1 (MQ=255)
                                                                                              tttCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGagag                                                <  2:853937/150‑1 (MQ=255)
                                                                                              tttCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGagag                                                <  1:445435/150‑1 (MQ=255)
                                                                                                  ctctAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCa                                            >  2:798737/1‑150 (MQ=255)
                                                                                                  ctctAAGTACTTACTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATCCGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCa                                            >  2:329855/1‑150 (MQ=255)
                                                                                                        gtacttagtacttTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCACCCCAGAAACGGTGAGAGTGCAGGTTCg                                      >  2:1308886/1‑150 (MQ=255)
                                                                                                         tacttagtacttTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCgg                                     >  1:314233/1‑150 (MQ=255)
                                                                                                          acttagtacttTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGc                                    >  1:355578/1‑150 (MQ=255)
                                                                                                          acttagtacttTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGc                                    <  2:78274/150‑1 (MQ=255)
                                                                                                           cttagtacttTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCa                                   <  1:280557/150‑1 (MQ=255)
                                                                                                           cttagtacttTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCa                                   <  2:208433/150‑1 (MQ=255)
                                                                                                                  cttTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTtctc                                                                         <  1:1107735/105‑1 (MQ=255)
                                                                                                                      tCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGt                        >  2:930941/1‑150 (MQ=255)
                                                                                                                        gCCACGGGCACGCAGTTCAGCCAGAACTTTATCTATGCATTTTTTATCGATTACCCGCATACATATAGCAGATANCCGCAGGGTGACAAAACGCTTCTCGCTCTAAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACggacttt                      >  2:575511/1‑145 (MQ=255)
                                                                                                                        gCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGttt                      <  2:1089791/150‑1 (MQ=255)
                                                                                                                        gCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGttt                      <  2:314233/150‑1 (MQ=255)
                                                                                                                         ccACGGGCACGCAGTTCAGCCAGACCTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCATCTCGCACTCAACCCAGAACCGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTa                     >  2:620258/1‑150 (MQ=255)
                                                                                                                         ccACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAg                    >  1:1219246/1‑150 (MQ=255)
                                                                                                                          cACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAg                    <  1:1252300/150‑1 (MQ=255)
                                                                                                                           aCGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGt                   <  1:617121/150‑1 (MQ=255)
                                                                                                                            cGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTc                  <  1:1407536/150‑1 (MQ=255)
                                                                                                                             gggCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTcg                 <  1:368236/150‑1 (MQ=255)
                                                                                                                              ggCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTcgc                >  2:661235/1‑150 (MQ=255)
                                                                                                                               gCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTcgcg               <  2:355578/150‑1 (MQ=255)
                                                                                                                                cACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGt              >  2:581546/1‑150 (MQ=255)
                                                                                                                                  cGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTc            <  2:92815/150‑1 (MQ=255)
                                                                                                                                  cGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTc            <  2:1116915/150‑1 (MQ=255)
                                                                                                                                  cGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTc            >  2:1081832/1‑150 (MQ=255)
                                                                                                                                      gTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTg        <  2:1113118/150‑1 (MQ=255)
                                                                                                                                       ttCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGc       <  2:1275225/150‑1 (MQ=255)
                                                                                                                                         cAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGt     >  2:175446/1‑150 (MQ=255)
                                                                                                                                          aGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTg    <  1:1107768/150‑1 (MQ=255)
                                                                                                                                          aGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTg    >  1:1279601/1‑150 (MQ=255)
                                                                                                                                          aGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTg    >  1:208655/1‑150 (MQ=255)
                                                                                                                                            ccAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTggg  <  1:930941/150‑1 (MQ=255)
                                                                                                                                            ccAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTggg  <  1:604348/150‑1 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
TAGTACGTTTGTTCTTCGTAGTGGTATAGAAGTGACCAGTACCAGCAGAAGAAACCAGCTTGATTTTCTCACGAATACCTTTAGCCATGATTTATTTCCTCTAAGTACTTAGTACTTTTCGCCACGGGCACGCAGTTCAGCCAGAACTGTATCGATGCCTTTTTTATCGATTACACGCATACCTTTAGCAGATACGCGCAGGGTGACAAAACGCTTCTCGCTCTCAACCCAGAAACGGTGAGAGTGCAGGTTCGGCAGGAAACGGCGTTTAGTCGCGTTCAGTGCGTGGG  >  NC_000913/3811330‑3811619

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: