Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,216,996 Δ2 bp [ymgC] [ymgC]

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,216,9950A.96.2% 399.9 / 10.1 106coding (248/249 nt)ymgCBlue light, low temperature and stress induced protein
Reads supporting (aligned to +/- strand):  ref base A (2/2);  new base . (54/48);  total (56/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01
*NC_0009131,216,9960G.96.2% 416.5 / 8.8 106coding (249/249 nt)ymgCBlue light, low temperature and stress induced protein
Reads supporting (aligned to +/- strand):  ref base G (2/2);  new base . (54/48);  total (56/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

Only 100 of 111 total aligned reads displayed.

GCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTT  >  NC_000913/1216859‑1217116
                                                                                                                                        ||                                                                                                                        
gCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaac                                                                                                                         >  1:927382/1‑137 (MQ=255)
          tGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAgg                                                                                                             >  2:694167/1‑139 (MQ=255)
          tGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAgg                                                                                                             >  2:143672/1‑139 (MQ=255)
          tGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGAAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGATCTCTT‑‑AACTTGAGTAAgg                                                                                                             >  1:589487/1‑139 (MQ=255)
           gCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGt                                                                                                            >  1:1433990/1‑139 (MQ=255)
             gcttAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTATCTT‑‑AACTTGAGTAAGGttt                                                                                                          <  1:1055010/137‑1 (MQ=255)
                        gAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACa                                                                                               <  2:1433990/139‑1 (MQ=255)
                            ccAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                           <  1:694167/139‑1 (MQ=255)
                            ccAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAAACTACACTAACGGATATCATTATTAATGATAGGTAATTCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                           <  1:1750431/139‑1 (MQ=255)
                                tgtgATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATc                                                                                             >  1:1075972/1‑133 (MQ=255)
                                tgtgATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATc                                                                                             <  2:1075972/133‑1 (MQ=255)
                                tgtgATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTg                                                                                        <  2:982468/138‑1 (MQ=255)
                                tgtgATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTg                                                                                        >  1:982468/1‑138 (MQ=255)
                                tgtgATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGt                                                                                       >  2:1277469/1‑139 (MQ=255)
                                          tctGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                           <  2:762046/125‑1 (MQ=255)
                                          tctGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTg                                                                                        >  1:338109/1‑128 (MQ=255)
                                          tctGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTg                                                                                        <  2:338109/128‑1 (MQ=255)
                                              aTCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAAc                                                                                                <  1:1162831/116‑1 (MQ=255)
                                              aTCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAAc                                                                                                >  2:1162831/1‑116 (MQ=255)
                                                        gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                               <  2:1402492/139‑1 (MQ=255)
                                                        gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                               >  1:90939/1‑139 (MQ=255)
                                                        gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATAACATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                           >  2:443459/1‑111 (MQ=255)
                                                        gCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATAACATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                           <  1:443459/111‑1 (MQ=255)
                                                               ttGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                >  1:1537675/1‑83 (MQ=255)
                                                               ttGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTa                                                                                                                <  2:1537675/83‑1 (MQ=255)
                                                                        aTTGAAATGCATAAGACTACACTAACTGATATCATTATTAATGACAGTTAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTATAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                               <  2:589487/139‑1 (MQ=255)
                                                                          tGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTGTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCt                                             >  1:1699120/1‑139 (MQ=255)
                                                                            aaaTGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTg                                                                 <  1:470573/117‑1 (MQ=255)
                                                                            aaaTGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTg                                                                 >  2:470573/1‑117 (MQ=255)
                                                                            aaaTGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCt                                             <  2:373569/137‑1 (MQ=255)
                                                                            aaaTGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCt                                             >  1:373569/1‑137 (MQ=255)
                                                                               tGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGc                                                     <  1:1645371/126‑1 (MQ=255)
                                                                                 cATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCt                                                                            >  1:839191/1‑103 (MQ=255)
                                                                                 cATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCt                                                                            <  2:839191/103‑1 (MQ=255)
                                                                                 cATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                           <  1:1782150/86‑1 (MQ=255)
                                                                                 cATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAg                                                                                           >  2:1782150/1‑86 (MQ=255)
                                                                                  aTAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCttatt                                     >  2:1381583/1‑139 (MQ=255)
                                                                                     aGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATg                                  <  1:1277469/139‑1 (MQ=255)
                                                                                      gACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGt                                 >  2:1089318/1‑139 (MQ=255)
                                                                                      gACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGt                                 <  1:1381583/139‑1 (MQ=255)
                                                                                          acacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGc                             >  2:190369/1‑139 (MQ=255)
                                                                                           cacTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTaac                                                                                                                         <  2:1536550/48‑3 (MQ=255)
                                                                                            acTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCaa                           >  2:877110/1‑139 (MQ=255)
                                                                                              tAACGGATATCATTATTAATGACAGGGGATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAata                         >  1:422645/1‑139 (MQ=255)
                                                                                              tAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCa                                                                                            >  2:1580876/1‑72 (MQ=255)
                                                                                              tAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCa                                                                                            <  1:1580876/72‑1 (MQ=255)
                                                                                                aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCa                                                                                >  1:1369325/1‑82 (MQ=255)
                                                                                                aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCa                                                                                <  2:1369325/82‑1 (MQ=255)
                                                                                                aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                               <  2:1763157/99‑1 (MQ=255)
                                                                                                aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                               >  1:1763157/1‑99 (MQ=255)
                                                                                                aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                               <  1:1017001/115‑1 (MQ=255)
                                                                                                aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                               >  2:1017001/1‑115 (MQ=255)
                                                                                                aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCttat                                      >  1:1204768/1‑124 (MQ=255)
                                                                                                aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCttat                                      <  2:1204768/124‑1 (MQ=255)
                                                                                                aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCa                            >  1:1281772/1‑134 (MQ=255)
                                                                                                aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATc                       >  1:427765/1‑139 (MQ=255)
                                                                                                aCGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATc                       >  1:387778/1‑139 (MQ=255)
                                                                                                    atatCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGgt                   <  1:387952/139‑1 (MQ=255)
                                                                                                     tatCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCt                                             <  2:1118040/114‑1 (MQ=255)
                                                                                                     tatCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCt                                             >  1:1118040/1‑114 (MQ=255)
                                                                                                      atCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTcc                               <  2:1016805/125‑1 (MQ=255)
                                                                                                      atCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTcc                               >  1:1016805/1‑125 (MQ=255)
                                                                                                      atCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGgtgt                 <  1:143672/139‑1 (MQ=255)
                                                                                                      atCATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGgtgt                 <  2:90939/139‑1 (MQ=255)
                                                                                                        cATTATTAATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCa               >  1:140448/1‑139 (MQ=255)
                                                                                                          ttattaATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTg                                                                                        >  2:179811/1‑64 (MQ=255)
                                                                                                          ttattaATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTg                                                                                        <  1:179811/64‑1 (MQ=255)
                                                                                                          ttattaATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCaa              >  1:395014/1‑138 (MQ=255)
                                                                                                          ttattaATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCaa              <  2:395014/138‑1 (MQ=255)
                                                                                                           tattaATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGTTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGa            >  2:46293/1‑139 (MQ=255)
                                                                                                             ttaATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTcc                               <  2:602577/118‑1 (MQ=255)
                                                                                                             ttaATGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTcc                               >  1:602577/1‑118 (MQ=255)
                                                                                                               aaTGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTtattat                                   <  2:601539/112‑1 (MQ=255)
                                                                                                               aaTGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTtattat                                   >  1:601539/1‑112 (MQ=255)
                                                                                                                aTGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAg                     >  2:632997/1‑125 (MQ=255)
                                                                                                                aTGACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAg                     <  1:632997/125‑1 (MQ=255)
                                                                                                                  gACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTa                                                 >  2:921629/1‑95 (MQ=255)
                                                                                                                  gACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTa                                                 <  1:921629/95‑1 (MQ=255)
                                                                                                                  gACAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGgaga     >  2:698352/1‑139 (MQ=255)
                                                                                                                   aCAGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGagag    >  1:922395/1‑139 (MQ=255)
                                                                                                                     aGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                               <  1:425932/78‑1 (MQ=255)
                                                                                                                     aGGGAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGtt  >  2:211860/1‑139 (MQ=255)
                                                                                                                       ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCa                                                                                            >  2:455912/1‑47 (MQ=255)
                                                                                                                       ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCAttt                                                          <  1:475837/81‑1 (MQ=255)
                                                                                                                       ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCAttt                                                          >  2:475837/1‑81 (MQ=255)
                                                                                                                       ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                               >  2:1115051/1‑92 (MQ=255)
                                                                                                                       ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATg                                               <  1:1115051/92‑1 (MQ=255)
                                                                                                                       ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCt                                             >  1:61486/1‑94 (MQ=255)
                                                                                                                       ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGga       >  1:459383/1‑132 (MQ=255)
                                                                                                                       ggAATCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGga       <  2:459383/132‑1 (MQ=255)
                                                                                                                         aaTCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCAc                                                                                  >  2:1153168/1‑55 (MQ=255)
                                                                                                                         aaTCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCAc                                                                                  <  1:1153168/55‑1 (MQ=255)
                                                                                                                         aaTCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCttatt                                     <  1:1009277/100‑1 (MQ=255)
                                                                                                                         aaTCCGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCttatt                                     >  2:1009277/1‑100 (MQ=255)
                                                                                                                             cGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTAt                                                                                     <  2:1195227/48‑1 (MQ=255)
                                                                                                                             cGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTcc                               >  1:1772738/1‑102 (MQ=255)
                                                                                                                             cGTGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTcc                               <  2:1772738/102‑1 (MQ=255)
                                                                                                                               tGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCt                                             >  1:1440039/1‑86 (MQ=255)
                                                                                                                               tGCTCTCTT‑‑AACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCt                                             <  2:1440039/86‑1 (MQ=255)
                                                                                                                                        tgAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAg                                                               >  1:1079962/2‑61 (MQ=255)
                                                                                                                                        ||                                                                                                                        
GCTTTACGCTTGCGTTTAAAAGATGAAACCAGTGTGATTCTTTCTGATCTCTATCTGCAAATATTGCAGTACATTGAAATGCATAAGACTACACTAACGGATATCATTATTAATGACAGGGAATCCGTGCTCTCTTAGAACTTGAGTAAGGTTTTAGAATAAACATCAGCTGTATCACCATGCTGATGCAAAGTGAGCATTTCAGGCGTTATGCTTTCTTATTATGTCCGCAATATCAGGTGTCAAGAATGGAGAGTT  >  NC_000913/1216859‑1217116

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