Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,107,876 A→C 7.2% I427S (ATT→AGT) ‡ speC ← ornithine decarboxylase, biosynthetic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,107,8760AC 7.2% 249.3 / 10.7 84I427S (ATT→AGT) ‡speCornithine decarboxylase, biosynthetic
Reads supporting (aligned to +/- strand):  ref base A (34/43);  new base C (4/2);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.59e-01

ATATCCTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAGGCGTTATT  >  NC_000913/3107728‑3108024
                                                                                                                                                    |                                                                                                                                                    
atatCCTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAAt                                                                                                                                                     <  2:427689/150‑1 (MQ=255)
     cTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCt                                                                                                                                                <  1:778534/150‑1 (MQ=255)
      ttCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGCCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAGCAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCtt                                                                                                                                               >  1:410878/1‑150 (MQ=255)
        cAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTgc                                                                                                                                             >  2:1589041/1‑150 (MQ=255)
          aaGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTgcgc                                                                                                                                           >  1:37709/1‑150 (MQ=255)
          aaGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTgcgc                                                                                                                                           <  2:1079038/150‑1 (MQ=255)
           aGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAAGACTGATGTCGGATAATCCTGCCTCAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTgcgcg                                                                                                                                          <  2:1277994/150‑1 (MQ=255)
           aGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGGCTATTCTTgcgcg                                                                                                                                          <  1:538411/150‑1 (MQ=255)
               gTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTc                                                                                                                                      <  1:819549/150‑1 (MQ=255)
                tGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCa                                                                                                                                     <  1:1303515/150‑1 (MQ=255)
                  ccACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGGCTATTCTTGCGCGCtgtaag                                                                                                                                   <  1:354713/150‑6 (MQ=255)
                  ccACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTAACAGCGCGCAAGAATAGCCTTGCGCGCTTCAAt                                                                                                                                   <  2:470292/150‑1 (MQ=255)
                   cACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATc                                                                                                                                  >  2:646400/1‑150 (MQ=255)
                           ccccGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTAtctc                                                                                                                          >  2:944868/1‑150 (MQ=255)
                            cccGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCa                                                                                                                         >  2:1017549/1‑150 (MQ=255)
                             ccGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCaa                                                                                                                        >  1:1012434/1‑150 (MQ=255)
                             acGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCaa                                                                                                                        <  1:1062805/149‑1 (MQ=255)
                              cGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAac                                                                                                                       <  1:537187/150‑1 (MQ=255)
                              cGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAac                                                                                                                       <  2:655392/150‑1 (MQ=255)
                               ggCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAaca                                                                                                                      <  1:1589041/150‑1 (MQ=255)
                                 cTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAacaca                                                                                                                    >  1:549472/1‑150 (MQ=255)
                                  gcAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGCGCGGGATAAACTGTAGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCAAATCTCAacacac                                                                                                                   <  2:1282960/149‑1 (MQ=255)
                                   cAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAACCGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACt                                                                                                                  >  2:1343661/1‑150 (MQ=255)
                                     aaaCTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCa                                                                                                                >  2:547378/1‑150 (MQ=255)
                                         tGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGccc                                                                                                            >  1:1361810/1‑150 (MQ=255)
                                                 cggcggTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGa                                                                                                    <  1:1374794/150‑1 (MQ=255)
                                                      gTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGccc                                                                                               <  2:184017/150‑1 (MQ=255)
                                                        cGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCAc                                                                                             <  1:1431918/150‑1 (MQ=255)
                                                        cGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCAc                                                                                             >  1:1312513/1‑150 (MQ=255)
                                                              cTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTccc                                                                                       >  2:429563/1‑150 (MQ=255)
                                                                 acacTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCtt                                                                                    <  1:109354/150‑1 (MQ=255)
                                                                   acTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAAGGGGCGGGATATACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCa                                                                                  <  2:1600527/150‑1 (MQ=255)
                                                                    cTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCAt                                                                                 >  1:1525644/1‑150 (MQ=255)
                                                                      gATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGa                                                                               >  1:1460146/1‑150 (MQ=255)
                                                                        tGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAAt                                                                             >  1:319224/1‑150 (MQ=255)
                                                                             gATAATCCTGCCACAATTTGCCTTCAAGAACGTGCGGGATATACGGGCGTAACAGCTTAGAGGGGGCAAGAATAGCCTTGGGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTgg                                                                        <  2:37709/150‑1 (MQ=255)
                                                                              aTAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGc                                                                       <  2:140384/150‑1 (MQ=255)
                                                                                 aTCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACCCTCAGCCCACAGCCGACCCCCACTCTCCCCTTCATGAATTTTGGCGtt                                                                    >  2:1547456/1‑150 (MQ=255)
                                                                                  tCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTa                                                                   <  1:944868/150‑1 (MQ=255)
                                                                                   ccTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTaa                                                                  <  1:916505/150‑1 (MQ=255)
                                                                                    cTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAAc                                                                 <  2:1370630/150‑1 (MQ=255)
                                                                                      gCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACAt                                                               >  2:730351/1‑150 (MQ=255)
                                                                                       ccACAATTTACCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATc                                                              <  1:1260622/150‑1 (MQ=255)
                                                                                        cacaATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGGGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATcc                                                             <  2:1095882/150‑1 (MQ=255)
                                                                                        cacaATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATcc                                                             <  2:1361810/150‑1 (MQ=255)
                                                                                        cacaATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATcc                                                             >  1:111693/1‑150 (MQ=255)
                                                                                        cacaATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATcc                                                             <  1:1364081/150‑1 (MQ=255)
                                                                                             tttGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTgc                                                        <  1:1533146/150‑1 (MQ=255)
                                                                                              ttGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTgca                                                       <  2:549472/150‑1 (MQ=255)
                                                                                                 ccATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTgcagca                                                    >  2:335991/1‑150 (MQ=255)
                                                                                                  cATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCaa                                                   >  2:288893/1‑150 (MQ=255)
                                                                                                    tCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAAc                                                 <  2:1312513/150‑1 (MQ=255)
                                                                                                    tCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAAc                                                 >  1:1312232/1‑150 (MQ=255)
                                                                                                    tCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAAc                                                 <  1:1573013/150‑1 (MQ=255)
                                                                                                     caacaaCGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACa                                                >  1:447680/1‑150 (MQ=255)
                                                                                                      aacaacGGGCGGGATAAACGGGCGGAACCGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAg                                               >  2:1450982/1‑150 (MQ=255)
                                                                                                      aacaacGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGtt                                                                    <  1:1045333/129‑1 (MQ=255)
                                                                                                      aacaacGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGtt                                                                    >  2:1045333/1‑129 (MQ=255)
                                                                                                      aacaacGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACa                                                                >  2:1063882/1‑133 (MQ=255)
                                                                                                      aacaacGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACa                                                                <  1:1063882/133‑1 (MQ=255)
                                                                                                       acaacGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGc                                              >  1:1336953/1‑150 (MQ=255)
                                                                                                         aacGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGTgg                                            <  2:447680/150‑1 (MQ=255)
                                                                                                         aacGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCgg                                            <  1:291252/150‑1 (MQ=255)
                                                                                                             ggcgggATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAg                                        >  1:1270108/1‑150 (MQ=255)
                                                                                                                 ggATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAg                                    >  2:969420/1‑150 (MQ=255)
                                                                                                                    tAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGTGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGCCGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGc                                 >  1:1625495/1‑150 (MQ=255)
                                                                                                                        cGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCACCACACTCAGGCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGt                             <  2:1041373/150‑1 (MQ=255)
                                                                                                                         gggCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTa                            <  1:1218574/150‑1 (MQ=255)
                                                                                                                          ggCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAg                           <  2:167489/150‑1 (MQ=255)
                                                                                                                                  aGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATgg                   <  2:1012434/150‑1 (MQ=255)
                                                                                                                                   gCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGa                  >  2:778455/1‑150 (MQ=255)
                                                                                                                                     ttACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGc                <  1:730351/150‑1 (MQ=255)
                                                                                                                                      tgaaGCGCGCAAGGCTATTCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCa               >  1:300566/4‑150 (MQ=255)
                                                                                                                                       gaaGCGCGCAAGGCTATTCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCAt              >  1:227227/3‑150 (MQ=255)
                                                                                                                                       aCAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTGTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCAt              >  2:872593/1‑150 (MQ=255)
                                                                                                                                        aaGCGCGCAAGGCTATTCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATa             >  1:1159312/2‑150 (MQ=255)
                                                                                                                                        aaGCGCGCAAGGCTATTCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAACAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATa             >  1:376437/2‑150 (MQ=255)
                                                                                                                                          gcgcgcAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATaaa           <  1:1273513/150‑1 (MQ=255)
                                                                                                                                          gcgcgcAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATaaa           >  1:663942/1‑150 (MQ=255)
                                                                                                                                           cgcgcATGATTAGCCTTGCGCGCTTAAATCCGAATCTCAACACACTCAGCCCACAGCCGACGCGCACTCTCCCCTACATGAATTTTGGCGTTAACATCGAGTGCAGCTTACAGGGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAg          <  2:1491839/150‑1 (MQ=255)
                                                                                                                                             cgcAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAGGc        <  2:410878/150‑1 (MQ=255)
                                                                                                                                              gcAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACCCTCAGCCCACAGACGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATACAGTGCAGCAAACAGCTGATACAACGGGCTGGTAGAAGCATGGAGCAGaaatcct       >  2:1298150/1‑146 (MQ=255)
                                                                                                                                              gcAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAGGCg       <  1:547378/150‑1 (MQ=255)
                                                                                                                                                  gAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTTCAGCAAACAGCGGATAGAAAGGGCTTGTAGAAGCATGGAGCATAAAGGCGTTAt   <  1:877850/150‑1 (MQ=255)
                                                                                                                                                   taTAGGCTTGAGCGCTTCAATCCCTATCTCAACACACTCAGCCGACAGCCGACGCCCACTCTCCGCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGATAGGGCTGGTAGAAGCATGGAGCATAAAGGCGTTAtt  <  2:599540/149‑1 (MQ=255)
                                                                                                                                                   aaTAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAGGCGTTAtt  <  1:969420/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATATCCTTCAAAGCCGTGCCACTTCGCCCCCGGCTCAAAACTGAAAAAACGGCGGTCGCTGGCTAACACTGATGTCGGATAATCCTGCCACAATTTGCCATCAACAACGGGCGGGATAAACGGGCGGAACAGCTTACAGCGCGCAAGAATAGCCTTGCGCGCTTCAATCCCTATCTCAACACACTCAGCCCACAGCCGACGCCCACTCTCCCCTTCATGAATTTTGGCGTTAACATCCAGTGCAGCAAACAGCGGATAGAAAGGGCTGGTAGAAGCATGGAGCATAAAGGCGTTATT  >  NC_000913/3107728‑3108024

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