Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,208,210 T→C 5.4% S63P (TCC→CCC)  ttdT → L‑tartrate/succinate antiporte

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,208,2100TC 5.4% 269.8 / 5.1 92S63P (TCC→CCC) ttdTL‑tartrate/succinate antiporte
Reads supporting (aligned to +/- strand):  ref base T (40/46);  new base C (3/2);  total (44/48)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.32e-01

GGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCTTTATGT  >  NC_000913/3208065‑3208357
                                                                                                                                                 |                                                                                                                                                   
ggCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTGTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACGGATCCTCGCACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCcaa                                                                                                                                                 >  1:185095/1‑149 (MQ=255)
 gCGGTCATCGCGATTATTGCTCTACTTCCGGTTCCAGCAGGGCTCGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGAACGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCcatc                                                                                                                                                >  2:303603/1‑150 (MQ=255)
 gCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCcatc                                                                                                                                                <  1:1010144/150‑1 (MQ=255)
    gTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCatcatc                                                                                                                                             >  2:1287755/1‑150 (MQ=255)
     tCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTTGTGGCGATGATGGATATACCCACCATCg                                                                                                                                            <  1:1389150/150‑1 (MQ=255)
       aTCGCCATTATTGCTCTACTTCAGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGAGCGCTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGATGGATATACCCACCATCGcc                                                                                                                                          <  1:828755/150‑1 (MQ=255)
        tCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCa                                                                                                                                         <  1:337867/150‑1 (MQ=255)
         cGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCAt                                                                                                                                        <  2:1395064/150‑1 (MQ=255)
         cGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCACCAt                                                                                                                                        <  2:476350/150‑1 (MQ=255)
          gCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATa                                                                                                                                       >  2:642804/1‑150 (MQ=255)
          gCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGATTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGAGGGTATATCCATCATCGCCATa                                                                                                                                       >  2:647392/1‑150 (MQ=255)
            cATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATAct                                                                                                                                     <  2:75869/150‑1 (MQ=255)
             attattGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCGCTCCTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATActc                                                                                                                                    <  1:309692/150‑1 (MQ=255)
                 ttGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACtctctc                                                                                                                                <  2:928450/150‑1 (MQ=255)
                    ctctACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCt                                                                                                                             >  1:672716/1‑150 (MQ=255)
                    ctctACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCt                                                                                                                             <  2:826072/150‑1 (MQ=255)
                              ggTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGtt                                                                                                                   >  1:152227/1‑150 (MQ=255)
                               gTTCCCGCAGGGCTGGAGAACCATTCCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGtta                                                                                                                  >  1:496791/1‑149 (MQ=255)
                                 tCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAg                                                                                                                >  2:791684/1‑150 (MQ=255)
                                  cccGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTAATGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCCATGGTGGGTATATACATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGc                                                                                                               >  1:747222/1‑150 (MQ=255)
                                   ccGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGcc                                                                                                              >  2:947022/1‑150 (MQ=255)
                                   ccGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGATGTTCAGcc                                                                                                              >  2:1308344/1‑150 (MQ=255)
                                   ccGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGATGGACTGATCCTCTATCCCGTGCCGGGTGCCGTGTTGGCGATGGTGGGTATATACATAATCTCCATACTCTCTCACTGGATGCTGTTAAGcc                                                                                                              >  2:68410/1‑150 (MQ=255)
                                    cGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGCACCTTGGACTGATCCTAGAACCCCTGCCGCGTGCCGTGGTGGCGATGGTGGGTATAGCCAAAATCGCCAAACTCACACCCTGGCAGCGGTTCAGccc                                                                                                             >  1:1558786/1‑150 (MQ=255)
                                    cGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGGGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGccc                                                                                                             >  1:807734/1‑150 (MQ=255)
                                      cAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCgg                                                                                                           <  1:570036/150‑1 (MQ=255)
                                      cAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGGGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCgg                                                                                                           <  2:931357/150‑1 (MQ=255)
                                         ggCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGagc                                                                                                        <  1:1345622/150‑1 (MQ=255)
                                          gCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTCGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCCGCCCGGagct                                                                                                       >  1:221091/1‑149 (MQ=255)
                                            tGGAGAACCATACCTGGCTCTACTTTGCAGTTATTACTGGCGTGATCGTTGGACTGATAATCGAACCCGTGCCCGGTGCCGTGGTGGCGATGGGGGGTATATCCATCATAGCAATTCTCTCACCCTGGCTGATTATCAGAACCGagcagc                                                                                                     >  2:639175/1‑150 (MQ=255)
                                             ggAGAACCATACCTGGCTCTACTTTGCCGTTTTTCCTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAgct                                                                                                    <  1:290533/150‑1 (MQ=255)
                                             ggAGAACCATACCTGGCGCTACTTAGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAgct                                                                                                    <  2:1131100/150‑1 (MQ=255)
                                               agaACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAgctcg                                                                                                  <  1:1708133/150‑1 (MQ=255)
                                                gaACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAgctcgc                                                                                                 >  2:1715089/1‑150 (MQ=255)
                                                     aTACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGc                                                                                            >  2:1221778/1‑150 (MQ=255)
                                                     aTACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGc                                                                                            >  2:1519927/1‑150 (MQ=255)
                                                      tACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGcc                                                                                           >  1:28898/1‑150 (MQ=255)
                                                         cTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCCTACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAgg                                                                                        >  2:929874/1‑150 (MQ=255)
                                                            gCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCtt                                                                                     >  1:1040885/1‑150 (MQ=255)
                                                            gCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCAGTGGTGGCGATGGTGGGTATATACATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCtt                                                                                     >  1:1055087/1‑150 (MQ=255)
                                                             ctctACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCttt                                                                                    >  2:640875/1‑150 (MQ=255)
                                                             ctctACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCttt                                                                                    >  1:1180787/1‑150 (MQ=255)
                                                                tACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTTCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTaaa                                                                                 <  1:870303/150‑1 (MQ=255)
                                                                 aCTTTGCCGTTTTTACTGGCGTGACCGCTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAAt                                                                                <  2:1465559/150‑1 (MQ=255)
                                                                  cTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAAtt                                                                               >  1:1088625/1‑150 (MQ=255)
                                                                  cTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAAtt                                                                               >  1:1566858/1‑150 (MQ=255)
                                                                    ttGCCGGTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCa                                                                             <  1:1646923/150‑1 (MQ=255)
                                                                     tGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCAc                                                                            <  1:947022/150‑1 (MQ=255)
                                                                       ccGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTg                                                                          <  2:56318/150‑1 (MQ=255)
                                                                       ccGTTTTTACTGGCGGGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTg                                                                          <  1:243537/150‑1 (MQ=255)
                                                                              tACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTTTATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATc                                                                   <  2:218005/150‑1 (MQ=255)
                                                                              tACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATc                                                                   >  2:173140/1‑150 (MQ=255)
                                                                                cTGGCGTGATCGTTTGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATccc                                                                 <  2:398492/150‑1 (MQ=255)
                                                                                    cGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCtctc                                                             <  2:450301/150‑1 (MQ=255)
                                                                                     gTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCg                                                            >  1:69470/1‑150 (MQ=255)
                                                                                     gTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGATCAGCCCGGAGCAGCTCGCTCATCCAGGCTTTAAATTCACTGCAAAATCCCTCTCg                                                            <  1:303603/150‑1 (MQ=255)
                                                                                       gATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTg                                                          <  2:1476603/150‑1 (MQ=255)
                                                                                         tCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTggg                                                        <  2:1088625/150‑1 (MQ=255)
                                                                                            ttGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCg                                                     <  2:672716/150‑1 (MQ=255)
                                                                                                 cTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTcc                                                <  2:425158/150‑1 (MQ=255)
                                                                                                 cTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCAGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTcc                                                >  1:853676/1‑150 (MQ=255)
                                                                                                   gATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCgg                                              >  1:1272985/1‑150 (MQ=255)
                                                                                                   gATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCgg                                              >  1:783131/1‑150 (MQ=255)
                                                                                                       cTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGtttt                                          <  1:1287755/150‑1 (MQ=255)
                                                                                                         cGAACCCGTGCCGTGTGCCTTGGTGGCGATGGTGGGTATATCAATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGTCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGtttttg                                        <  1:68410/150‑2 (MQ=255)
                                                                                                         cGAACCCATGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTc                                        <  1:944928/150‑1 (MQ=255)
                                                                                                          gAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAATATCCCTCTCGTGGGCCGTTTCCGGTTTTTCt                                       <  1:647392/150‑1 (MQ=255)
                                                                                                          gAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCt                                       <  1:173140/150‑1 (MQ=255)
                                                                                                           aaCCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTa                                      <  2:358611/150‑1 (MQ=255)
                                                                                                           aaCCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTa                                      <  2:13112/150‑1 (MQ=255)
                                                                                                           aaCCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTa                                      <  2:1180787/150‑1 (MQ=255)
                                                                                                           aaCCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTa                                      >  2:474137/1‑150 (MQ=255)
                                                                                                            aCCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTaa                                     <  1:1191545/150‑1 (MQ=255)
                                                                                                               cGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTc                                  <  2:982317/150‑1 (MQ=255)
                                                                                                                 tGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCgg                                <  2:1566858/150‑1 (MQ=255)
                                                                                                                         gCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTACTTCGGTAATCTgg                        >  2:1142645/1‑150 (MQ=255)
                                                                                                                          ccGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGc                       >  1:992913/1‑150 (MQ=255)
                                                                                                                            gtggtgGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTg                     >  2:1407455/1‑150 (MQ=255)
                                                                                                                             gtatggCGATGATGGATATACCCACCATCGCCATACTCTCTCCCTGGATGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGGTTTAAATTCACTGCAAAATCCCTCACGTGGGCCGATTCCGGTGTTTCTAATTCGGTTATCTGGCTGa                    >  1:1107589/4‑150 (MQ=255)
                                                                                                                                tggCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGAttt                 <  1:1715089/150‑1 (MQ=255)
                                                                                                                                tggCGATGATGGATATACCCACCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCCGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGAttt                 >  2:177635/1‑150 (MQ=255)
                                                                                                                                tggCGATGATGGATATACCCACCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGAttt                 >  2:928860/1‑150 (MQ=255)
                                                                                                                                   cGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTcg              >  2:57763/1‑150 (MQ=255)
                                                                                                                                   cGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTcg              <  1:1221778/150‑1 (MQ=255)
                                                                                                                                   cGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTcg              >  1:1321846/1‑150 (MQ=255)
                                                                                                                                   cGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTcg              <  1:1519927/150‑1 (MQ=255)
                                                                                                                                     aTGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCATCCAGGCTTTAAATTCACTGCAAAATCCCTGTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGGTATCTGGCTGATTTTcgcc            <  2:152227/150‑1 (MQ=255)
                                                                                                                                       ggtggGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGTGCAGCTCGCTCAGCCAGGCTTTATATTCACTTCAAAATCCCTCGCGTGGGCGGTGTCCGGTTTGTCTAATTCGGTTATCTGTCTGATTTTcgccgc          >  1:311865/1‑150 (MQ=255)
                                                                                                                                        gtggGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCgccgcc         <  1:1131293/150‑1 (MQ=255)
                                                                                                                                           ggTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAATTCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCttt      <  1:219968/150‑1 (MQ=255)
                                                                                                                                           ggTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCttt      >  1:1365182/1‑150 (MQ=255)
                                                                                                                                            gTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCTTTa     <  2:783131/150‑1 (MQ=255)
                                                                                                                                               tatCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCTTTATGt  >  1:346615/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GGCGGTCATCGCCATTATTGCTCTACTTCCGGTTCCCGCAGGGCTGGAGAACCATACCTGGCTCTACTTTGCCGTTTTTACTGGCGTGATCGTTGGACTGATCCTCGAACCCGTGCCGGGTGCCGTGGTGGCGATGGTGGGTATATCCATCATCGCCATACTCTCTCCCTGGCTGCTGTTCAGCCCGGAGCAGCTCGCTCAGCCAGGCTTTAAATTCACTGCAAAATCCCTCTCGTGGGCCGTTTCCGGTTTTTCTAATTCGGTTATCTGGCTGATTTTCGCCGCCTTTATGT  >  NC_000913/3208065‑3208357

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