Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,613,461 C→T 12.5% G415S (GGT→AGT)  sad ← succinate semialdehyde dehydrogenase, NAD(P)+‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,613,4610CT12.5% 264.2 / 25.5 96G415S (GGT→AGT) sadsuccinate semialdehyde dehydrogenase, NAD(P)+‑dependent
Reads supporting (aligned to +/- strand):  ref base C (44/40);  new base T (6/6);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.56e-01

GTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATT  >  NC_000913/1613314‑1613609
                                                                                                                                                   |                                                                                                                                                    
ttCAGATCCGGTCTGTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATAGATAAACACCCCACCg                                                                                                                                                     >  2:114985/2‑149 (MQ=255)
    gATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCAtt                                                                                                                                                 <  1:948148/149‑1 (MQ=255)
      tCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTc                                                                                                                                                >  3:633811/1‑148 (MQ=255)
      tCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTc                                                                                                                                                >  3:277235/1‑148 (MQ=255)
       ccGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCa                                                                                                                                              <  2:595998/149‑1 (MQ=255)
        cGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACTGCATTCCAg                                                                                                                                             >  2:1059022/1‑149 (MQ=255)
        cGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAg                                                                                                                                             >  1:752690/1‑149 (MQ=255)
         ggTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAg                                                                                                                                             <  4:310568/148‑1 (MQ=255)
          gTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAGTAACCATTGATAAACACCCCACTGCATTCCAGAc                                                                                                                                           <  2:1044536/149‑1 (MQ=255)
                  acacCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTgccgcc                                                                                                                                   >  1:192453/1‑149 (MQ=255)
                         cTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACTGCATTCCAGACGTGCCGCCAtctgtc                                                                                                                            <  2:730577/149‑1 (MQ=255)
                         cTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCAtctgt                                                                                                                             <  4:1005537/148‑1 (MQ=255)
                            gATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACTGCATTCCAGACGTGCCGCCAtctgtctg                                                                                                                          >  4:831634/1‑148 (MQ=255)
                             atatTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACTGCATTCCAGACGTGCCGCCATCTGTCTGGc                                                                                                                        >  2:730153/1‑149 (MQ=255)
                               atTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCt                                                                                                                      <  2:1014014/149‑1 (MQ=255)
                               agtACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGccc                                                                                                                       >  3:934043/3‑148 (MQ=255)
                                    aGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGtt                                                                                                                 >  1:681356/1‑149 (MQ=255)
                                      aaTTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTc                                                                                                               <  2:681356/149‑1 (MQ=255)
                                        ttCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCg                                                                                                              >  4:219261/1‑148 (MQ=255)
                                         tCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGt                                                                                                             >  4:929356/1‑148 (MQ=255)
                                             gtAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTg                                                                                                        <  1:330872/149‑1 (MQ=255)
                                                aGCCGAAGTGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGt                                                                                                      >  3:467700/1‑148 (MQ=255)
                                                   cGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTaaa                                                                                                   >  4:887428/1‑148 (MQ=255)
                                                     aaaTGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACTGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTaaaaa                                                                                                 >  3:678365/1‑148 (MQ=255)
                                                        tGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATgg                                                                                              <  4:934043/148‑1 (MQ=255)
                                                           gAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGc                                                                                          <  2:496312/149‑1 (MQ=255)
                                                              aGCTCACGACCAAAGCCACTCTTTTTCAAGCCACCAAAGGCCACTCGCGCGTCGGTGGCACAATAACAATTGATAAACACCCCATCGCTTTCCAGACGTGCCGCAATCTGTCTGGCCTATGAATCGTCAGTGGTAAAAATGGTCGCTg                                                                                        >  4:656973/1‑148 (MQ=255)
                                                               gCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTgg                                                                                                                         >  2:954278/1‑114 (MQ=255)
                                                               gCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTgg                                                                                                                         <  1:954278/114‑1 (MQ=255)
                                                                 tCACGCCCACAGCCACTCTTTTTCCCGCCACCCAAGCCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGCGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAACTGGTCGCTGaaa                                                                                     <  3:419901/148‑1 (MQ=255)
                                                                 tCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAg                                                                                    >  2:636504/1‑149 (MQ=255)
                                                                    cGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGc                                                                                  <  4:749026/148‑1 (MQ=255)
                                                                    cGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGc                                                                                  >  3:706665/1‑148 (MQ=255)
                                                                     gACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACTGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGcc                                                                                 <  3:831634/148‑1 (MQ=255)
                                                                      aCCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCg                                                                                <  4:372216/148‑1 (MQ=255)
                                                                       ccAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGaa                                                                              >  1:96392/1‑149 (MQ=255)
                                                                       ccAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGaa                                                                              <  1:707866/149‑1 (MQ=255)
                                                                            gCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCAc                                                                         <  2:417681/149‑1 (MQ=255)
                                                                             ccACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACt                                                                        <  1:1071514/149‑1 (MQ=255)
                                                                             ccACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACt                                                                        >  1:731386/1‑149 (MQ=255)
                                                                               aCTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTa                                                                       <  4:460338/148‑1 (MQ=255)
                                                                                ctctTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTc                                                                           <  4:225127/143‑1 (MQ=255)
                                                                                ctctTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTc                                                                           >  3:225127/1‑143 (MQ=255)
                                                                                ctctTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTAt                                                                      >  4:357459/1‑148 (MQ=255)
                                                                                ctctTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTAt                                                                      >  4:964830/1‑148 (MQ=255)
                                                                                ctctTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTAt                                                                      >  3:908729/1‑148 (MQ=255)
                                                                                 tctTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATc                                                                     <  3:401968/148‑1 (MQ=255)
                                                                                   tttttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACGATTGATAAACACCCCACTGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCAtt                                                                  >  2:140619/1‑149 (MQ=255)
                                                                                   ttgttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTGAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTTTCAtt                                                                  >  2:386771/4‑149 (MQ=255)
                                                                                    ttttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACTGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTa                                                                 <  1:730153/149‑1 (MQ=255)
                                                                                    ttttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCa                                                                    >  3:553145/1‑146 (MQ=255)
                                                                                    ttttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCAtt                                                                  >  3:129688/1‑148 (MQ=255)
                                                                                    ttttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTa                                                                 <  2:752690/149‑1 (MQ=255)
                                                                                    ttttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACGCCCTATCa                                                                    <  4:553145/146‑1 (MQ=255)
                                                                                    tattCACGCCACCAAAGGCCACTCGCGCGTCGCTAGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCAtt                                                                  <  3:634509/146‑1 (MQ=255)
                                                                                     tttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATAGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGT‑AAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGc                                                               <  1:114985/149‑1 (MQ=255)
                                                                                      ttCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGc                                                               <  2:365682/149‑1 (MQ=255)
                                                                                            ccaccaAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGtt                                                          >  3:731290/1‑148 (MQ=255)
                                                                                              accaAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTcc                                                        <  3:357459/148‑1 (MQ=255)
                                                                                                   aGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTcc                                                        >  2:270820/1‑143 (MQ=255)
                                                                                                   aGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTcc                                                        <  1:270820/143‑1 (MQ=255)
                                                                                                   aGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCa                                                  <  2:998439/149‑1 (MQ=255)
                                                                                                        aCTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGtt                                              <  4:706665/148‑1 (MQ=255)
                                                                                                        aCTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGtt                                              <  4:949697/148‑1 (MQ=255)
                                                                                                          tCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTg                                           >  1:435712/1‑149 (MQ=255)
                                                                                                            gcgcgTCGCTGGCACAATAACCATTGATAAACACCCCACTGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGc                                          <  4:678365/148‑1 (MQ=255)
                                                                                                                 tCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCttt                                    <  1:558791/149‑1 (MQ=255)
                                                                                                                  cGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTc                                   <  2:96392/149‑1 (MQ=255)
                                                                                                                    cTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCg                                  <  4:467700/148‑1 (MQ=255)
                                                                                                                      ggCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCaa                               >  2:850423/1‑149 (MQ=255)
                                                                                                                      ggCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCa                                >  4:747140/1‑148 (MQ=255)
                                                                                                                       gCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCaa                               <  3:373931/148‑1 (MQ=255)
                                                                                                                        cacaATAACCATTGATAAACACCCCACTGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATg                             <  1:1059022/149‑1 (MQ=255)
                                                                                                                        cacaATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAAt                              >  3:26346/1‑148 (MQ=255)
                                                                                                                          caATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCg                                                                                           <  4:1050377/85‑1 (MQ=255)
                                                                                                                          caATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCg                                                                                           >  3:1050377/1‑85 (MQ=255)
                                                                                                                           aaTAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGt                           >  4:285360/1‑148 (MQ=255)
                                                                                                                                ccATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGAATCTTTCGCAATGGTGATTg                      >  4:607625/1‑148 (MQ=255)
                                                                                                                                  aTTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGcc                    >  3:1038854/1‑148 (MQ=255)
                                                                                                                                  aTTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGcc                    >  4:624054/1‑148 (MQ=255)
                                                                                                                                   ttGATAAACACCCCACTGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCg                   >  3:776664/1‑148 (MQ=255)
                                                                                                                                   ttGATAAACACCCCACCGGATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCg                   >  3:646204/1‑148 (MQ=255)
                                                                                                                                   ttGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAgg                                                                                   >  4:288749/1‑84 (MQ=255)
                                                                                                                                   ttGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAgg                                                                                   <  3:288749/84‑1 (MQ=255)
                                                                                                                                   ttGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCg                   >  3:203720/1‑148 (MQ=255)
                                                                                                                                   ttGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGc                  >  2:197051/1‑149 (MQ=255)
                                                                                                                                      aTAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCa                                                       <  1:1001751/109‑1 (MQ=255)
                                                                                                                                      aTAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCa                                                       >  2:1001751/1‑109 (MQ=255)
                                                                                                                                        aaaCACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCaa                >  3:13327/1‑146 (MQ=255)
                                                                                                                                        aaaCACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCaa                <  4:13327/146‑1 (MQ=255)
                                                                                                                                            acCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGcc         <  2:435712/149‑1 (MQ=255)
                                                                                                                                            acCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGcc         >  2:210174/1‑149 (MQ=255)
                                                                                                                                                cACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCa                                                                          >  3:143833/1‑80 (MQ=255)
                                                                                                                                                cACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCa                                                                          <  4:143833/80‑1 (MQ=255)
                                                                                                                                                 aCCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTc                                                                           <  3:439283/78‑1 (MQ=255)
                                                                                                                                                 aCCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTc                                                                           >  4:439283/1‑78 (MQ=255)
                                                                                                                                                  ccGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACa    <  4:129688/148‑1 (MQ=255)
                                                                                                                                                   cGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGATCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACAtt  >  1:78045/1‑149 (MQ=255)
                                                                                                                                                   cGCATTCCAGACGTGCCGCCATCTGTCTGGCCTCTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACAt   >  3:474965/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTCAGATCCGGTCTTTCCACACCGTCTGGATATTACAGAATTCGTGTAAGCCGAAATGGGAAAGCTCACGACCAAAGCCACTCTTTTTCACGCCACCAAAGGCCACTCGCGCGTCGCTGGCACAATAACCATTGATAAACACCCCACCGCATTCCAGACGTGCCGCCATCTGTCTGGCCTGTGTTTCGTCAGTGGTAAAAATGGTCGCTGAAAGGCCGAACTCACTATCATTAGCCAGTTCCAGTGCATGTTCTGCATCTTTCGCAATGGTGATTGCCGCAACGGGGCCAAACATT  >  NC_000913/1613314‑1613609

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