Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,431,170 A→C 6.1% R82R (CGT→CGG tsaC ← threonylcarbamoyl‑AMP synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,431,1700AC 6.1% 287.6 / 8.2 95R82R (CGT→CGGtsaCthreonylcarbamoyl‑AMP synthase
Reads supporting (aligned to +/- strand):  ref base A (49/40);  new base C (2/4);  total (51/44)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

CATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCAGGATCGCACCC  >  NC_000913/3431022‑3431316
                                                                                                                                                    |                                                                                                                                                  
cATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCAc                                                                                                                                                   <  2:204652/150‑1 (MQ=255)
cATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCAc                                                                                                                                                   >  1:1268264/1‑150 (MQ=255)
   aaGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCt                                                                                                                                                >  2:1496166/1‑150 (MQ=255)
    aGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGTGAAACCATTTTTTCCCGCTg                                                                                                                                               <  2:882466/150‑1 (MQ=255)
    aGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTg                                                                                                                                               <  2:364294/150‑1 (MQ=255)
    aGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTg                                                                                                                                               >  2:1132567/1‑150 (MQ=255)
    aGCCTGGCACAAAACAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTg                                                                                                                                               <  1:1393573/150‑1 (MQ=255)
     gCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCTAATCAAAGCGGCCCGTCAACCAGCGCTTTGTTTTCGCAGTCGCGTGAAATACAAAGGTGACAGGACCTGTCCAGCGTGAAACCATTTTTTCCCGCTgg                                                                                                                                              <  2:1063816/150‑2 (MQ=255)
      ccTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGTGAAACCATTTTTTCCCGCTggc                                                                                                                                             <  1:708478/150‑3 (MQ=255)
       cTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCATGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGTGAAACCATTTTTTCCCGCTggcc                                                                                                                                            <  2:129993/150‑4 (MQ=255)
       cTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGCCAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCAc                                                                                                                                            >  2:1007180/1‑150 (MQ=255)
        tGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACg                                                                                                                                           <  1:1228884/150‑1 (MQ=255)
           cacaAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCGCgtca                                                                                                                                        >  1:1044359/1‑150 (MQ=255)
           cacaAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACgtca                                                                                                                                        >  1:1402392/1‑150 (MQ=255)
           cacaAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACgtca                                                                                                                                        >  2:978966/1‑150 (MQ=255)
                 gCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAAc                                                                                                                                  <  1:910776/150‑1 (MQ=255)
                 gCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAAc                                                                                                                                  <  2:977976/150‑1 (MQ=255)
                     ccaccaACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATgg                                                                                                                              <  1:195891/150‑1 (MQ=255)
                        ccaACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGgtgt                                                                                                                           >  1:307869/1‑150 (MQ=255)
                           aCGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACg                                                                                                                                           >  1:244809/1‑131 (MQ=255)
                           aCGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACg                                                                                                                                           <  2:244809/131‑1 (MQ=255)
                           aCGGATGGGCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCTCAGGCGCTGGAAAGACAAAGGTGACAGGTCCTGGCCAGCGGGAAAATATGGTTTCACGCTGCACGTCAGTCAAGATGGTGtcat                                                                                                                        <  2:316869/150‑1 (MQ=255)
                             ggATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGtcatca                                                                                                                      <  1:978966/150‑1 (MQ=255)
                              gATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCaa                                                                                                                     >  1:1172821/1‑150 (MQ=255)
                                    cGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAgg                                                                                                               >  2:93749/1‑150 (MQ=255)
                                        gACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGttt                                                                                                           >  2:284246/1‑150 (MQ=255)
                                            cGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGTGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgc                                                                                                       <  1:139449/150‑1 (MQ=255)
                                            cGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgc                                                                                                       <  1:97748/150‑1 (MQ=255)
                                            cGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgc                                                                                                       <  1:92142/150‑1 (MQ=255)
                                              tACAGCAAGCGAATCAAAGCGGCCCGTCAAGCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgctg                                                                                                     <  2:1055760/150‑1 (MQ=255)
                                              tACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgctg                                                                                                     >  1:629459/1‑150 (MQ=255)
                                               aCAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGTAAATATGGTTTCACGCTGCACGTCAGTCAACATGGGGTCCTCAATATAGGGTTTAAgctgc                                                                                                    >  2:1389295/1‑150 (MQ=255)
                                               aCAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgctgc                                                                                                    >  2:250683/1‑150 (MQ=255)
                                                cAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAATACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAgctgct                                                                                                   <  2:307869/150‑1 (MQ=255)
                                                  gCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCg                                                                                                 >  2:680434/1‑150 (MQ=255)
                                                  gCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATCGTGTCATCAATATAGGGTTTAAGCTGCACg                                                                                                 >  1:608166/1‑150 (MQ=255)
                                                     aGCGAATCAAAGCGGCCCGTCAACCAGCGCGCTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGCCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCCTaa                                                                                              >  2:1289820/1‑150 (MQ=255)
                                                      gCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCTGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTGAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAAt                                                                                             <  2:284311/150‑1 (MQ=255)
                                                       cGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAAtt                                                                                            >  1:912456/1‑150 (MQ=255)
                                                          aTCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATAtgc                                                                                         <  2:808702/150‑1 (MQ=255)
                                                             aaaGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCTGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTtgctgc                                                                                      <  2:216688/150‑1 (MQ=255)
                                                             aaaGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTtgctgc                                                                                      >  2:1497875/1‑150 (MQ=255)
                                                              aaGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCg                                                                                     <  1:131421/150‑1 (MQ=255)
                                                                  ggCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTa                                                                                 >  2:671336/1‑150 (MQ=255)
                                                                  ggCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTa                                                                                 <  1:1282653/150‑1 (MQ=255)
                                                                     ccGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTaaaa                                                                              >  1:416338/1‑150 (MQ=255)
                                                                            ccAGCGCGCTGTTGTCCCAGGCGCGGGAAATACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCACTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGccc                                                                       >  2:1129491/1‑150 (MQ=255)
                                                                              aGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGAGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCt                                                                     >  1:7992/1‑150 (MQ=255)
                                                                                      gttgtCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAAcc                                                             >  2:1488482/1‑150 (MQ=255)
                                                                                       ttgtCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCg                                                            <  1:680434/150‑1 (MQ=255)
                                                                                       ttgtCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCg                                                            <  1:1497875/150‑1 (MQ=255)
                                                                                         gtCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGa                                                          >  2:518727/1‑150 (MQ=255)
                                                                                           cGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCg                                                            >  1:1136321/1‑146 (MQ=255)
                                                                                           cGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCg                                                            <  2:1136321/146‑1 (MQ=255)
                                                                                            gCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGc                                                       >  1:375965/1‑150 (MQ=255)
                                                                                            gCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGc                                                       <  1:749130/150‑1 (MQ=255)
                                                                                               ggCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGt                                                    <  2:1269925/150‑1 (MQ=255)
                                                                                                 cgcgGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTAt                                                                  >  2:485419/1‑134 (MQ=255)
                                                                                                 cgcgGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTAt                                                                  <  1:485419/134‑1 (MQ=255)
                                                                                                  gcgGGAAAGACAAAGGGTACAGGACCTGGCCAGCTGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGtttt                                                 <  2:973601/150‑1 (MQ=255)
                                                                                                   cgGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTTCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTa                                                <  1:1213572/150‑1 (MQ=255)
                                                                                                     ggAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAAc                                              <  2:1402392/150‑1 (MQ=255)
                                                                                                     ggAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAAc                                              <  1:93749/150‑1 (MQ=255)
                                                                                                      gAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACt                                             >  1:39972/1‑150 (MQ=255)
                                                                                                       aaaGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAATATCAGCCCCTTATCAACCGGACGCTGTTTTAACTc                                            <  1:1047662/150‑1 (MQ=255)
                                                                                                       aaaGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTc                                            <  2:935372/150‑1 (MQ=255)
                                                                                                        aaGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTca                                           >  2:489049/1‑149 (MQ=255)
                                                                                                          gACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCaa                                         >  2:1027445/1‑150 (MQ=255)
                                                                                                           aCAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAAc                                        <  2:994437/150‑1 (MQ=255)
                                                                                                            cAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACa                                       >  1:1116925/1‑150 (MQ=255)
                                                                                                            cAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACa                                       >  2:733478/1‑150 (MQ=255)
                                                                                                                  tGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCa                                 <  2:1005871/150‑1 (MQ=255)
                                                                                                                  tGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCa                                 <  1:875738/150‑1 (MQ=255)
                                                                                                                   gACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCAt                                <  1:613589/150‑1 (MQ=255)
                                                                                                                   gACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCAt                                >  2:805146/1‑150 (MQ=255)
                                                                                                                      aGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCAc                             >  1:1134486/1‑150 (MQ=255)
                                                                                                                       ggACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCGCCTTATCAACCGGACGCTGTTTTCACTCCGACAGTCGCATCagt                            >  2:1053930/1‑148 (MQ=255)
                                                                                                                       ggACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCAct                            >  2:313982/1‑150 (MQ=255)
                                                                                                                          ccTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCActgct                         >  2:1526230/1‑150 (MQ=255)
                                                                                                                           cTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCActgctg                        >  1:584741/1‑150 (MQ=255)
                                                                                                                            tGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATTGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGt                       >  2:1350968/1‑150 (MQ=255)
                                                                                                                            cgtgcAGCGTGAAACCATTTTTTCCCGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGt                       >  1:775977/5‑150 (MQ=255)
                                                                                                                             ggCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGtt                      >  1:516254/1‑150 (MQ=255)
                                                                                                                              tgcAGCGTGAAACCATTTTTTCCCGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGttt                     >  1:549029/3‑150 (MQ=255)
                                                                                                                               ccAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTc                    >  2:944783/1‑150 (MQ=255)
                                                                                                                                cAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCg                   <  2:908710/150‑1 (MQ=255)
                                                                                                                                cAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCg                   <  2:327755/150‑1 (MQ=255)
                                                                                                                                     ggAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAAc                                              <  2:1312290/118‑1 (MQ=255)
                                                                                                                                     ggAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAAc                                              >  1:1312290/1‑118 (MQ=255)
                                                                                                                                     ggAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATc              >  2:311791/1‑150 (MQ=255)
                                                                                                                                        aaaaaTGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCGGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCAgg           >  2:553810/1‑150 (MQ=255)
                                                                                                                                         aaaaTGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCAGGa          >  2:354057/1‑150 (MQ=255)
                                                                                                                                           aaTGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAATAGTCGCATCACTGCTGTTTCGCTATCAGGATc        >  2:137729/1‑150 (MQ=255)
                                                                                                                                                tttCACGCTGCACGTAAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATAAAAATCAGCCCCTTATCAACCCGACGCTCTTTTAACTCCAACAGTCGCATCACTGCTGATTCCCTATCAGGATCGCAcc   >  2:357391/1‑150 (MQ=255)
                                                                                                                                                 ttCACGCTGCGCGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCAGGATCGCAccc  <  2:1044359/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CATAAGCCTGGCACAAAGCAACCACCAACGGATGGTCGGTGACTCGTACAGCAAGCGAATCAAAGCGGCCCGTCAACCAGCGCGGTGTTGTCGCAGGCGCGGGAAAGACAAAGGTGACAGGACCTGGCCAGCGGGAAAAAATGGTTTCACGCTGCACGTCAGTCAACATGGTGTCATCAATATAGGGTTTAAGCTGCTCGTAATTTGCTGCGATTAAAATCAGCCCCTTATCAACCGGACGCTGTTTTAACTCCAACAGTCGCATCACTGCTGTTTCGCTATCAGGATCGCACCC  >  NC_000913/3431022‑3431316

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