Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,873,917 G→C 5.9% H26D (CAT→GAT) ‡ cysC ← adenosine 5'‑phosphosulfate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,873,9170GC 5.9% 534.5 / 17.7 150H26D (CAT→GAT) ‡cysCadenosine 5'‑phosphosulfate kinase
Reads supporting (aligned to +/- strand):  ref base G (70/71);  new base C (3/6);  total (73/77)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.79e-01

Only 100 of 152 total aligned reads displayed.

GATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTT  >  NC_000913/2873770‑2874065
                                                                                                                                                   |                                                                                                                                                    
  tCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGAcc                                                                                                                                                   <  4:607850/149‑1 (MQ=255)
     cTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtg                                                                                                                                                <  4:1180103/149‑1 (MQ=255)
     cTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtg                                                                                                                                                <  2:1014880/149‑1 (MQ=255)
        cATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtggtg                                                                                                                                             >  3:1183369/1‑149 (MQ=255)
         aTAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGgtgt                                                                                                                                            >  1:565481/1‑149 (MQ=255)
          tAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTa                                                                                                                                           <  1:15043/149‑1 (MQ=255)
          tAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTa                                                                                                                                           >  1:607248/1‑149 (MQ=255)
           aaTCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTAg                                                                                                                                          <  4:1162575/149‑1 (MQ=255)
            aTCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTAgt                                                                                                                                         <  1:594223/149‑2 (MQ=255)
            aTCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGc                                                                                                                                         <  4:1060968/149‑1 (MQ=255)
               cGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTcg                                                                                                                                      >  4:80357/1‑149 (MQ=255)
                     gAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAGCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGTACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgt                                                                                                                                >  2:378420/1‑149 (MQ=255)
                     gAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgt                                                                                                                                <  2:238900/149‑1 (MQ=255)
                      aaCATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                               >  1:1255660/1‑149 (MQ=255)
                           tGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACGCCACGATGACCGTGGTGTAGCTCGCGTTGTTGcaca                                                                                                                          <  3:1343962/149‑1 (MQ=255)
                           tGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGcaca                                                                                                                          <  1:1434489/149‑1 (MQ=255)
                           tGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGcaca                                                                                                                          <  2:518241/149‑1 (MQ=255)
                              cGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTg                                                                                                                       <  3:725401/149‑1 (MQ=255)
                               gCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGa                                                                                                                      >  4:1234951/1‑149 (MQ=255)
                                 cATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtggt                                                                                                                                              <  4:605703/123‑1 (MQ=255)
                                 cATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtggt                                                                                                                                              >  3:605703/1‑123 (MQ=255)
                                 cATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGAcc                                                                                                                    <  1:892820/149‑1 (MQ=255)
                                     cagcagATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGAt                                                                                                                <  2:1466687/149‑1 (MQ=255)
                                     cagcagATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGAt                                                                                                                <  4:439675/149‑1 (MQ=255)
                                       gcagATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATgg                                                                                                              >  3:390046/1‑149 (MQ=255)
                                           aTACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                               <  4:1360241/128‑1 (MQ=255)
                                           aTACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                               >  3:1360241/1‑128 (MQ=255)
                                             aCGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGc                                                                                                        <  3:80357/149‑1 (MQ=255)
                                             aCGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGc                                                                                                        >  1:1138882/1‑149 (MQ=255)
                                                tACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAg                                                                                                     <  2:1179063/149‑1 (MQ=255)
                                                   tGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtg                                                                                                                                                >  2:977449/1‑103 (MQ=255)
                                                   tGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtg                                                                                                                                                <  1:977449/103‑1 (MQ=255)
                                                   tGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAgacg                                                                                                  >  4:830099/1‑149 (MQ=255)
                                                   tGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAgacg                                                                                                  >  2:108317/1‑149 (MQ=255)
                                                     aCGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAgacgac                                                                                                >  4:1529991/1‑149 (MQ=255)
                                                     aCGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAgacgac                                                                                                <  1:18228/149‑1 (MQ=255)
                                                       gCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGt                                                                                              <  2:898533/149‑1 (MQ=255)
                                                             tttATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGt                                                                                        <  1:570736/149‑1 (MQ=255)
                                                             tttATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGt                                                                                        <  2:1515515/149‑1 (MQ=255)
                                                             tttATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGCGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGt                                                                                        <  2:1299766/149‑1 (MQ=255)
                                                               tATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCa                                                                                      <  1:644611/149‑1 (MQ=255)
                                                                aTGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCTCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCAt                                                                                     <  2:866082/149‑1 (MQ=255)
                                                                aTGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCAt                                                                                     >  3:1294952/1‑149 (MQ=255)
                                                                aTGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCAt                                                                                     <  4:292736/149‑1 (MQ=255)
                                                                   tAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTAt                                                                                                          <  4:1150685/125‑1 (MQ=255)
                                                                   tAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCa                                                                                  <  4:684714/149‑1 (MQ=255)
                                                                    aaCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTAt                                                                                                          >  3:1465927/1‑124 (MQ=255)
                                                                    aaCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTAt                                                                                                          <  4:1465927/124‑1 (MQ=255)
                                                                     aCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGtttt                                                                                           >  1:1409715/1‑138 (MQ=255)
                                                                     aCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAgc                                                                                >  3:41972/1‑149 (MQ=255)
                                                                       gCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAgcgc                                                                              >  2:761974/1‑149 (MQ=255)
                                                                             tccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCattat                                                                        <  4:41972/149‑1 (MQ=255)
                                                                             tccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCattat                                                                        >  3:1373270/1‑149 (MQ=255)
                                                                             tccAGCGCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCattatt                                                                       >  1:1341212/1‑149 (MQ=255)
                                                                               cAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgt                                                                                                                                <  4:1317973/91‑1 (MQ=255)
                                                                               cAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAttattt                                                                      >  1:1469807/1‑149 (MQ=255)
                                                                                 gcgcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAttatttat                                                                    >  1:296986/1‑149 (MQ=255)
                                                                                 gcgcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAttatttat                                                                    >  2:1486757/1‑149 (MQ=255)
                                                                                  cgcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGt                                                                                              <  3:701894/122‑1 (MQ=255)
                                                                                  cgcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGt                                                                                              >  4:701894/1‑122 (MQ=255)
                                                                                    ccccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCCCGATGACCGTGGTGTAGCTCTCGTTGTTGCACAGCGACCGGATGGCGATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATccc                                                                 <  1:1083629/149‑1 (MQ=255)
                                                                                    ccccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATccc                                                                 <  3:323939/149‑1 (MQ=255)
                                                                                    ccccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATccc                                                                 >  4:1056113/1‑149 (MQ=255)
                                                                                    ccccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATccc                                                                 >  2:242235/1‑149 (MQ=255)
                                                                                     cccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATcccc                                                                <  3:1529991/149‑1 (MQ=255)
                                                                                      ccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgtt                                                                                                                                  <  4:1121908/82‑1 (MQ=255)
                                                                                      ccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgtt                                                                                                                                  >  3:1121908/1‑82 (MQ=255)
                                                                                      ccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATccccc                                                               <  2:565481/149‑1 (MQ=255)
                                                                                            aCCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCaaa                                                         <  1:367256/149‑1 (MQ=255)
                                                                                              cGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTTGTGTAGCTCGCCTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATc                                                       <  2:1134991/149‑1 (MQ=255)
                                                                                              cGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATc                                                       <  4:788051/149‑1 (MQ=255)
                                                                                                    tttACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcgcgcgc                                                 >  4:1195477/1‑149 (MQ=255)
                                                                                                         cGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTc                                                                                          <  4:290233/103‑1 (MQ=255)
                                                                                                          ggACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCa                                           >  3:391940/1‑149 (MQ=255)
                                                                                                              ccGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGc                                       >  3:563695/1‑149 (MQ=255)
                                                                                                              ccGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGc                                       >  3:864532/1‑149 (MQ=255)
                                                                                                                ggAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCgg                                     >  2:1407451/1‑149 (MQ=255)
                                                                                                                   gagGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGaaa                                  <  1:242235/149‑1 (MQ=255)
                                                                                                                     ggCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGt                                <  1:1322507/149‑1 (MQ=255)
                                                                                                                          gTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGcga                           >  1:1532236/1‑149 (MQ=255)
                                                                                                                           tAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGcgac                          <  1:1013952/149‑1 (MQ=255)
                                                                                                                           tAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGcgac                          >  1:830922/1‑149 (MQ=255)
                                                                                                                               ccACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAAc                      >  3:1369297/1‑149 (MQ=255)
                                                                                                                               ccACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAAc                      >  1:1198652/1‑149 (MQ=255)
                                                                                                                                 acaGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCa                    >  2:1113634/1‑149 (MQ=255)
                                                                                                                                 acaGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCa                    >  1:201805/1‑149 (MQ=255)
                                                                                                                                  caGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCag                   <  2:1421068/149‑1 (MQ=255)
                                                                                                                                   aGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCaga                  >  2:1110044/1‑149 (MQ=255)
                                                                                                                                   aGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCaga                  <  3:350021/149‑1 (MQ=255)
                                                                                                                                        tACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCAtt             >  4:952855/1‑149 (MQ=255)
                                                                                                                                        tACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCAtt             <  2:593030/149‑1 (MQ=255)
                                                                                                                                           accacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGc                                       >  4:836749/1‑120 (MQ=255)
                                                                                                                                           accacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGc                                       <  3:836749/120‑1 (MQ=255)
                                                                                                                                           accacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCaa          >  3:499003/1‑149 (MQ=255)
                                                                                                                                             cacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAtt        >  1:250946/1‑149 (MQ=255)
                                                                                                                                              acGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTc       <  2:1138882/149‑1 (MQ=255)
                                                                                                                                              acGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTc       >  1:1170192/1‑149 (MQ=255)
                                                                                                                                                gATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATcccc                                                                <  1:1116098/90‑1 (MQ=255)
                                                                                                                                                gATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCa     <  2:1394653/149‑1 (MQ=255)
                                                                                                                                                   gACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGtt  >  1:809592/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCAGTT  >  NC_000913/2873770‑2874065

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: