Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 993,258 G→T 5.0% intergenic (+24/+19) pepN → / ← ssuB aminopeptidase N/aliphatic sulfonate ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913993,2580GT 5.0% 500.9 / 10.9 139intergenic (+24/+19)pepN/ssuBaminopeptidase N/aliphatic sulfonate ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (69/63);  new base T (3/4);  total (72/67)
Fisher's exact test for biased strand distribution p-value = 7.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 140 total aligned reads displayed.

GATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCAAA  >  NC_000913/993111‑993403
                                                                                                                                                   |                                                                                                                                                 
gATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAgc                                                                                                                                                  >  1:577887/1‑149 (MQ=255)
        tGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCAt                                                                                                                                          <  2:1259947/149‑1 (MQ=255)
         gAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGcgccgc                                                                                                                                             >  4:351245/1‑145 (MQ=255)
         gAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGAACTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCAtt                                                                                                                                         <  1:384240/149‑1 (MQ=255)
           aaCGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCGCTATTGCCGCCATTCg                                                                                                                                       <  4:1091101/149‑1 (MQ=255)
            aCGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCgg                                                                                                                                      >  1:146617/1‑149 (MQ=255)
            aCGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCgg                                                                                                                                      <  1:1269351/149‑1 (MQ=255)
              gTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCgggg                                                                                                                                    <  2:1012604/149‑1 (MQ=255)
               ttACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGa                                                                                                                                   >  4:1405436/1‑149 (MQ=255)
               ttACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGa                                                                                                                                   <  4:1011499/149‑1 (MQ=255)
                     tGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAc                                                                                                                             <  4:61335/149‑1 (MQ=255)
                     tGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAc                                                                                                                             <  1:62716/149‑1 (MQ=255)
                         aaaCGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAGTAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCt                                                                                                                         >  1:468508/1‑149 (MQ=255)
                         aaaCGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCt                                                                                                                         >  3:341663/1‑149 (MQ=255)
                         aaaCGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCt                                                                                                                         >  4:1417596/1‑149 (MQ=255)
                         aaaCGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCt                                                                                                                         >  1:318258/1‑149 (MQ=255)
                            cGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGtt                                                                                                                      <  3:448122/149‑1 (MQ=255)
                              tCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCg                                                                                                                                       <  3:807814/130‑1 (MQ=255)
                              tCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGtttt                                                                                                                    <  1:59402/149‑1 (MQ=255)
                               cAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTtc                                                                                                                   <  4:192656/149‑1 (MQ=255)
                                aGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAc                                                                                                                             <  2:1402827/138‑1 (MQ=255)
                                aGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTAc                                                                                                                             >  1:1402827/1‑138 (MQ=255)
                                 ggAGAAAATGTGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTtctc                                                                                                                 >  1:679346/1‑149 (MQ=255)
                                  gagaAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTtctcc                                                                                                                <  4:1358327/149‑2 (MQ=255)
                                  gagaAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCa                                                                                                                >  4:627887/1‑149 (MQ=255)
                                  gagaAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCa                                                                                                                >  2:881911/1‑149 (MQ=255)
                                            gcgcgGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTc                                                                                                      <  3:452282/149‑1 (MQ=255)
                                            gcgcgGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTc                                                                                                      <  2:1450806/149‑1 (MQ=255)
                                             cgcgGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCa                                                                                                     <  4:624980/149‑1 (MQ=255)
                                               cggcgCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGa                                                                                                   <  2:882184/149‑1 (MQ=255)
                                               cggcgCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGa                                                                                                   <  4:1213544/149‑1 (MQ=255)
                                               cggcgCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGa                                                                                                   <  1:1448965/149‑1 (MQ=255)
                                                ggcgCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGAt                                                                                                  >  2:828508/1‑149 (MQ=255)
                                                    cTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcac                                                                                              <  4:223812/149‑1 (MQ=255)
                                                     tGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcacc                                                                                             >  2:1198220/1‑149 (MQ=255)
                                                      ggAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcacca                                                                                            >  3:1409523/1‑149 (MQ=255)
                                                       gAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcaccac                                                                                           >  3:1228602/1‑149 (MQ=255)
                                                         aCAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGt                                                                                         <  2:1058528/149‑1 (MQ=255)
                                                         aCAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGt                                                                                         <  2:1423118/149‑1 (MQ=255)
                                                             ttGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCa                                                                                     >  3:686336/1‑149 (MQ=255)
                                                                aaaGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCa                                                                                     >  4:1212820/1‑146 (MQ=255)
                                                                aaaGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCa                                                                                     <  3:1212820/146‑1 (MQ=255)
                                                                     gCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGt                                                                             <  4:640359/149‑1 (MQ=255)
                                                                       tGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTg                                                                           >  2:1383961/1‑149 (MQ=255)
                                                                        ggAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGc                                                                          <  4:1414605/149‑1 (MQ=255)
                                                                        ggAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGc                                                                          <  4:938702/149‑1 (MQ=255)
                                                                          aaaaTCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCaa                                                                        <  1:164444/149‑1 (MQ=255)
                                                                          aaaaTCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCAGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCaa                                                                        <  4:175888/149‑1 (MQ=255)
                                                                               ctctctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCt                                                                   >  2:621219/1‑149 (MQ=255)
                                                                                tctctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTc                                                                  >  3:1111186/1‑149 (MQ=255)
                                                                                tctctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAAGAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTAACCACGTTGCATCACCCGTTGCAACACCTc                                                                  >  1:1031454/1‑149 (MQ=255)
                                                                                 ctctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTcc                                                                 <  2:703369/149‑1 (MQ=255)
                                                                                 ctctGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTcc                                                                 >  2:562041/1‑149 (MQ=255)
                                                                                    tGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCa                                                                                     <  4:1114001/126‑1 (MQ=255)
                                                                                    tGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCt                                                              >  3:50435/1‑149 (MQ=255)
                                                                                    tGGCGAACTGTACGAGAAGAGAAGTAAAGCACTGCCTTGATAAAGAACAGAATGGCGGCAATAGCGCCGCCATTCGGGGAAATTACCCCTGTTTTCGCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCt                                                              >  2:615198/1‑149 (MQ=255)
                                                                                      gCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTc                                                            <  4:1289109/149‑1 (MQ=255)
                                                                                       cGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTcc                                                           <  1:1383961/149‑1 (MQ=255)
                                                                                        gATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCa                                                          >  3:928455/1‑149 (MQ=255)
                                                                                             gTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTcc                                                     >  4:840964/1‑149 (MQ=255)
                                                                                               aCGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGc                                                   <  1:1198220/149‑1 (MQ=255)
                                                                                                cGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCt                                                  <  2:220919/149‑1 (MQ=255)
                                                                                                   gaagaTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCgg                                                                                                                                      <  2:560640/62‑1 (MQ=255)
                                                                                                   gaagaTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAAt                                               >  4:710397/1‑149 (MQ=255)
                                                                                                       aTAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCa                                           >  4:809484/1‑149 (MQ=255)
                                                                                                         aaCTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTcc                                                                 >  4:888274/1‑125 (MQ=255)
                                                                                                         aaCTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCAcc                                         <  4:341663/149‑1 (MQ=255)
                                                                                                           cTAAAGCACTGGCTTGATAAATAACCGAATGGTGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGa                                       >  4:1448235/1‑149 (MQ=255)
                                                                                                            tAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCt                                                              >  3:971019/1‑125 (MQ=255)
                                                                                                            tAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATc                                              <  2:1153250/141‑1 (MQ=255)
                                                                                                            tAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATc                                              >  1:1153250/1‑141 (MQ=255)
                                                                                                            tAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGAc                                      >  3:238499/1‑149 (MQ=255)
                                                                                                                  aCTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAAt                                <  4:956639/149‑1 (MQ=255)
                                                                                                                   cTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATc                               >  2:408084/1‑149 (MQ=255)
                                                                                                                    tGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTaa                                 <  1:926291/146‑1 (MQ=255)
                                                                                                                      gCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTtct                                                                                                                  >  1:1573693/1‑63 (MQ=255)
                                                                                                                      gCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTtct                                                                                                                  <  2:1573693/63‑1 (MQ=255)
                                                                                                                      gCTTGATAAATAACCGAATGGAGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAg                                                                                                               >  1:928797/1‑66 (MQ=255)
                                                                                                                      gCTTGATAAATAACCGAATGGAGGCAATAGCGCCACCATTCGGGGAATTTACCCCTGTTTTCTCAg                                                                                                               <  2:928797/66‑1 (MQ=255)
                                                                                                                        ttGATAAATAACCGAATGGTGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCa                                                                                     >  3:1297013/1‑90 (MQ=255)
                                                                                                                        ttGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGAcgc                          >  4:776502/1‑149 (MQ=255)
                                                                                                                          gATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCgg                        >  1:774244/1‑149 (MQ=255)
                                                                                                                          gATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCgg                        >  1:924202/1‑149 (MQ=255)
                                                                                                                              aataaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACgg                    >  1:67945/1‑149 (MQ=255)
                                                                                                                                 aaCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACggggg                 <  3:210787/149‑1 (MQ=255)
                                                                                                                                  aCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGa                >  4:1156866/1‑149 (MQ=255)
                                                                                                                                  aCCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAAAACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGa                <  1:1040643/149‑1 (MQ=255)
                                                                                                                                   ccGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGat               >  1:1509645/1‑149 (MQ=255)
                                                                                                                                    cGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCg                                                    <  2:1519059/111‑1 (MQ=255)
                                                                                                                                    cGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCg                                                    >  1:1519059/1‑111 (MQ=255)
                                                                                                                                     gAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcac                                                                                              >  2:1043164/1‑68 (MQ=255)
                                                                                                                                     gAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTcac                                                                                              <  1:1043164/68‑1 (MQ=255)
                                                                                                                                     gAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGc                                                                          >  3:1412874/1‑88 (MQ=255)
                                                                                                                                     gAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGatat             <  4:662042/149‑1 (MQ=255)
                                                                                                                                      aaTGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCg                                                    <  1:757446/109‑1 (MQ=39)
                                                                                                                                      aaTGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCg                                                    >  2:757446/1‑109 (MQ=39)
                                                                                                                                      aaTGGCGGCGCTATTGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATc            >  1:1431971/1‑149 (MQ=255)
                                                                                                                                          gcggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACt        >  4:1082799/1‑149 (MQ=255)
                                                                                                                                            ggcAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGt      <  2:1006337/149‑1 (MQ=255)
                                                                                                                                               aaTAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCaa   >  2:586912/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GATTCGCCTGAAACGTTACGATGCCAAACGTCAGGAGAAAATGCGCGCGGCGCTGGAACAGTTGAAAGGGCTGGAAAATCTCTCTGGCGATCTGTACGAGAAGATAACTAAAGCACTGGCTTGATAAATAACCGAATGGCGGCAATAGCGCCGCCATTCGGGGAATTTACCCCTGTTTTCTCAGGCGAGTTTCAGATTCACCACGTTGCATCACCCGTTGCAACACCTCCGCTTCCAGCTCCGCTAATCGCACCGACCCTAATCGACGCGGACGGGGGATATCCACTGTCAAA  >  NC_000913/993111‑993403

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: