Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,873,915 A→G 7.5% H26H (CAT→CAC) ‡ cysC ← adenosine 5'‑phosphosulfate kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,873,9150AG 7.5% 544.2 / 24.7 172H26H (CAT→CAC) ‡cysCadenosine 5'‑phosphosulfate kinase
Reads supporting (aligned to +/- strand):  ref base A (81/78);  new base G (5/8);  total (86/86)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.84e-01

Only 100 of 176 total aligned reads displayed.

CGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCA  >  NC_000913/2873767‑2874062
                                                                                                                                                    |                                                                                                                                                   
cGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGa                                                                                                                                                     >  4:785707/1‑149 (MQ=255)
cGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGa                                                                                                                                                     >  3:1222050/1‑149 (MQ=255)
 gagaTCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGAt                                                                                                                                                    <  4:1289911/149‑1 (MQ=255)
 gagaTCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGAt                                                                                                                                                    >  3:757772/1‑149 (MQ=255)
 gagaTCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGAt                                                                                                                                                    <  2:96254/149‑1 (MQ=255)
         tccATAACCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGGGTTTATGTAACGCCTCCTCCAGCGCCCCCGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAACACTACACCACGATAACCctgg                                                                                                                                            >  2:636525/3‑149 (MQ=255)
             tAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTa                                                                                                                                        <  1:611360/149‑1 (MQ=255)
             tAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTa                                                                                                                                        <  2:1424209/149‑1 (MQ=255)
             tAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTa                                                                                                                                        <  2:1218900/149‑1 (MQ=255)
                 ccGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTc                                                                                                                                    >  3:1047960/1‑149 (MQ=255)
                       cGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttg                                                                                                                              <  1:1007114/149‑1 (MQ=255)
                        gAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGAcc                                                                                                                                                >  1:992398/1‑130 (MQ=255)
                        gAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgt                                                                                                                             <  2:814831/149‑1 (MQ=255)
                        gAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgt                                                                                                                             >  2:741887/1‑149 (MQ=255)
                         aaCATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                            <  3:457191/149‑1 (MQ=255)
                         aaCATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                            <  2:670984/149‑1 (MQ=255)
                         aaCATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                            <  1:1018178/149‑1 (MQ=255)
                         aaCATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCgttgtt                                                                                                                            <  2:100981/149‑1 (MQ=255)
                           cATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGc                                                                                                                          >  4:703453/1‑149 (MQ=255)
                              tGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGcaca                                                                                                                       <  2:1104480/149‑1 (MQ=255)
                               gTCGCCATCCAGCAGATACGTACTGACGCCGAGTTGATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGCCCGTGGTGTAGCTCGCGTTGTTGCACAg                                                                                                                      >  4:1119047/1‑149 (MQ=255)
                                tCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtgg                                                                                                                                            >  3:782327/1‑126 (MQ=255)
                                   ccATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGAc                                                                                                                  >  2:51760/1‑149 (MQ=255)
                                    cATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGAcc                                                                                                                 >  2:43545/1‑149 (MQ=255)
                                      tCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCgg                                                                                                               >  1:244010/1‑149 (MQ=255)
                                        cagcagATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACggtggt                                                                                                                                           >  4:950945/1‑119 (MQ=255)
                                        cagcagATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGAt                                                                                                             <  3:1087495/149‑1 (MQ=255)
                                        cagcagATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGAt                                                                                                             >  4:158471/1‑149 (MQ=255)
                                        cagcagATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGAt                                                                                                             >  4:601121/1‑149 (MQ=255)
                                        cagcagATACGTACTGACGCCGAATTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCgtggt                                                                                                                                           <  3:950945/119‑1 (MQ=255)
                                              aTACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCg                                                                                                                >  1:1473531/1‑140 (MQ=255)
                                              aTACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTNt                                                                                                       <  3:340928/149‑1 (MQ=255)
                                                  gTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCa                                                                                                   >  4:677512/1‑149 (MQ=255)
                                                          gCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGt                                                                                           >  1:756042/1‑149 (MQ=255)
                                                           ccGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGca                                                                                                                         <  2:872685/118‑1 (MQ=255)
                                                           ccGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGtt                                                                                          >  3:161574/1‑149 (MQ=255)
                                                            cGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGttt                                                                                         <  3:1324656/149‑1 (MQ=255)
                                                                tttATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTc                                                                                                                                    <  4:387259/102‑1 (MQ=255)
                                                                tttATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTc                                                                                                                                    >  3:387259/1‑102 (MQ=255)
                                                                   aTGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTa                                                                                                                                        <  2:173929/95‑1 (MQ=255)
                                                                    tGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATg                                                                                 <  1:91434/149‑1 (MQ=255)
                                                                      tAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCa                                                                               >  1:1296970/1‑149 (MQ=255)
                                                                       aaCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTAgt                                                                                                                                      >  3:13053/1‑92 (MQ=39)
                                                                       aaCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTAgt                                                                                                                                      <  4:13053/93‑2 (MQ=39)
                                                                       aaCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGGTCATCGTGGTGTAgt                                                                                                                                      <  4:1160992/93‑2 (MQ=38)
                                                                        aCGCCTCCTCCGGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAgc                                                                             >  3:477177/1‑149 (MQ=255)
                                                                        aCGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCa                                                                                                   <  1:255353/127‑1 (MQ=255)
                                                                          gCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAgcgc                                                                           >  4:394751/1‑149 (MQ=255)
                                                                            ctcctcCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCa                                                                         >  2:1482476/1‑149 (MQ=255)
                                                                               ctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCatta                                                                      >  1:1233222/1‑149 (MQ=255)
                                                                               ctccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGCCATGCAGCGCCatta                                                                      >  2:530468/1‑149 (MQ=255)
                                                                                tccAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCattat                                                                     <  1:185789/149‑1 (MQ=255)
                                                                                    gcgcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAttatttat                                                                 >  3:1184702/1‑149 (MQ=255)
                                                                                    gcgcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAttatttat                                                                 <  1:793128/149‑1 (MQ=255)
                                                                                    gcgcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCAttatttat                                                                 <  2:231810/149‑1 (MQ=255)
                                                                                      gcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATcc                                                               <  4:1415237/149‑1 (MQ=255)
                                                                                      gcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATcc                                                               <  2:1296970/149‑1 (MQ=255)
                                                                                      gcCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATcc                                                               <  4:732266/149‑1 (MQ=255)
                                                                                        cccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATcccc                                                             >  1:503824/1‑149 (MQ=255)
                                                                                         ccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATccccc                                                            <  1:159539/149‑1 (MQ=255)
                                                                                         ccGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATccccc                                                            <  1:881621/149‑1 (MQ=255)
                                                                                                 cGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATc                                                    >  1:609348/1‑149 (MQ=255)
                                                                                                  gTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcg                                                   >  4:1111749/1‑149 (MQ=255)
                                                                                                  gTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcg                                                   >  3:879309/1‑149 (MQ=255)
                                                                                                   ttGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATACCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcgc                                                  >  1:439207/1‑149 (MQ=255)
                                                                                                     gATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcgcgc                                                >  1:807352/1‑149 (MQ=255)
                                                                                                      aTTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATcgcgcg                                               >  1:1035824/1‑149 (MQ=255)
                                                                                                            cGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATcccc                                                             >  3:764361/1‑129 (MQ=255)
                                                                                                            cGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATcccc                                                             <  4:764361/129‑1 (MQ=255)
                                                                                                              gACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAg                                       >  1:868587/1‑149 (MQ=255)
                                                                                                                 ccGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGc                                    <  1:1401763/149‑1 (MQ=255)
                                                                                                                 ccGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGc                                    <  3:1350002/149‑1 (MQ=255)
                                                                                                                      gagGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGaaa                               <  1:599431/149‑1 (MQ=255)
                                                                                                                         gCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTg                            <  1:687794/149‑1 (MQ=255)
                                                                                                                         gCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTg                            <  4:879309/149‑1 (MQ=255)
                                                                                                                             gTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGcga                        >  1:827528/1‑149 (MQ=255)
                                                                                                                             gTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGcga                        <  3:394751/149‑1 (MQ=255)
                                                                                                                              tAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGcgac                       >  4:287226/1‑149 (MQ=255)
                                                                                                                                 aCCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGaa                    >  2:439606/1‑149 (MQ=255)
                                                                                                                                   cacaGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCatta                                                                      <  3:156138/97‑1 (MQ=255)
                                                                                                                                    acaGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCa                 >  3:796417/1‑149 (MQ=255)
                                                                                                                                      aGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCaga               <  3:965916/149‑1 (MQ=255)
                                                                                                                                       gCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCagag              <  3:1108779/149‑1 (MQ=255)
                                                                                                                                         aCTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCa            <  4:1184702/149‑1 (MQ=255)
                                                                                                                                         aCTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCa            >  1:130965/1‑149 (MQ=255)
                                                                                                                                          cTACACCACGGTCATCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCAt           >  2:289717/1‑149 (MQ=255)
                                                                                                                                           tACACCACGGTCATCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCAtt          >  1:1452248/1‑149 (MQ=255)
                                                                                                                                           tACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATc                                                                <  1:321259/95‑1 (MQ=255)
                                                                                                                                             caccacGGTCATCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCa        >  3:12792/1‑149 (MQ=255)
                                                                                                                                              accacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCaa       <  2:756042/149‑1 (MQ=255)
                                                                                                                                               ccacGATGACCGTGGTGTAGCTCTCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAt      >  1:390543/1‑149 (MQ=255)
                                                                                                                                               ccacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTg                            <  4:1333185/127‑1 (MQ=255)
                                                                                                                                               ccacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAt      >  4:465047/1‑149 (MQ=255)
                                                                                                                                               ccacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAt      <  2:366894/149‑1 (MQ=255)
                                                                                                                                               ccacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAt      >  4:105262/1‑149 (MQ=255)
                                                                                                                                                cacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAtt     >  4:269226/1‑149 (MQ=255)
                                                                                                                                                cacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAtt     >  2:883644/1‑149 (MQ=255)
                                                                                                                                                cacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAtt     >  1:1355897/1‑149 (MQ=255)
                                                                                                                                                cacGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAAtt     >  3:394350/1‑149 (MQ=255)
                                                                                                                                                   gATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCa  <  1:43545/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CGAGATCGCTGCATAATCCGTGGCGAACATTGTCGCCATCCAGCAGATACGTACTGACGCCGAGTTTATGTAACGCCTCCTCCAGCGCCCCGGCGACCGTTGATTTACCGGACCCGGAGAGGCCGGTAAACCACAGCACTACACCACGATGACCGTGGTGTAGCTCGCGTTGTTGCACAGTGACCGGATGGCTATGCCAGACGACGTTTTCGTCATGCAGCGCCATTATTTATCCCCCAGCAAATCGCGCGCGCCCCAGTGCGGAAAGTGGCGACGAACCAGAGCATTCAATTCCA  >  NC_000913/2873767‑2874062

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: